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La Gravifique Einsteinienne 


Six conférences données à l'Institut Henri Poincaré 


PAR 
Tx. DE DONDER 


Professeur à l’Université Libre de Bruxelles 


Ces pages renferment un exposé succinct des six Conférences que 
j'ai eu l’honneur de donner à l’Institut Henri Poincaré en novembre 
et décembre 1920. Elles avaient respectivement pour objet : 

I. Le champ gravifique massique. 

II. Le champ gravifique électromagnétique. 

III. Application à la mécanique ondulatoire. 

IV. Electrodynamique des corps en mouvement. 

V. Electromagnétostriction et thermodynamique relativistes. 

VI. Synthèse. 

Je me suis efforcé de montrer que la relativité générale fournit 
un instrument parfaitement adapté à l'étude de ces problèmes. Con- 
trairement à une opinion assez répandue, la gravifique einsteinienne 
ne doit pas être modifiée : comprise dans un sens suffisamment large, 
elle demeure en harmonie avec les théories les plus modernes de la 
mécanique ondulatoire. 


La Gravifique einsteinienne 


1. LES ÉQUATIONS FONDAMENTALES DU CHAMP GRAVIFIQUE. — 
Considérons une fonction A? ne dépendant que des g2ê, gxf, .…. Sa 
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variance par rapport à un changement quelconque de variables x,, 
Ya, Xa, Xa est celle d’un multiplicateur ou facteur de densité. Nous 
appellons 9 la fonction caractéristique gravifique. 

Considérons, en outre, une fonction ‘t, de même variance, mais 
qui peut dépendre, indépendamment des g2Ë, g2f#, d’autres fonctions 
telles w%, À,, etc. Cette fonction sera explicitée plus loin; elle sera 
désignée sous le nom de fonchion caractéristique phénoménale. 

Le principe variationnel fondamental de la gravifique consiste à 
égaler à zéro, les dix dérivées variationnelles par rapport aux gzf. 
On obtient ainsi les dix équations fondamentales de la gravifique ; 
à savoir 

2 
@) PE = 0. 


# e Æ La ® ù LA - 
La dérivée variationnelle PT s'écrit : 
dpt 


() Du, Rat D Js d à d? CR Vies 
5g%8 CE age dx; ag") Le dxjdx, ag*Pir 
Posons 
QUE NU? es 9 re 
(3) Cas Deré ORNE 
d’où 


Nous appellerons Ge le fenseur gravifique covariant symétrique, 
Gag le tenseur phénoménal covariant symétrique ou simplement /e 
lenseur symétrique. 

Désignons par C l’invariant de courbure et par a et b deux constantes 
universelles. Prenons pour 9 la valeur 
(5) AU = (a + DC)V— &. 


En effectuant les opérations indiquées en (9 II), on obtient 


(6) —; (a + bC)gss + be = Ts 
où k 
O8 
(7 To 
AV er 
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et Cag sont les composantes du tenseur de Riemann bien connu. 


Remarque. — Te principe variationnel tel que nous venons de 
l'exposer revient évidemment à généraliser le principe de Hamilton, 
c'est-à-diie à annuler la variation 


(8) Ô 1. (UE + ON) dx, dx,dx;dx, = 0 
ts 


O étant une portion d’espace-temps sur les frontières de laquelle 
les variations doivent s’annuler. 


2. LES IDENTITÉS GRAVIFIQUES. — En appliquant à la fonction 99 
les identités de la gravifique (1), on obtient 


do SAUCE 1. ONU 
et A es le nt AS — 


ou encore, en vertu de (3), 


0 RE 
Æ 11 AE. — 
(x0) noie C0 0 
où 
(x1) GE — go. 
8. THÉORÈME DE L'IMPULSION ET DE L'ÉNERGIE. — Les dix équa- 


tions (4) jointes aux quatre identités (10), donnent immédiatement 
les quatre équations 
2, 7! 


GÉNEE Ü RARE = 17, a 2 
(12) Bt y +28" %y—0 


où nous avons posé, comme en (24 Il), 
(13) ü, = LR 


Nous dirons que les quatre équations (12) expriment le théorème 
de l'impulsion et de l'énergie. On peut aussi exprimer ce théorème 
en disant que la force généralisée F, est nulle. 


(x) Th. DE DonNDER. Théorie invariantive du Calcul des varialions. Bulletin Ac. R. Belg. CI. des 
Sciences. 5° série, t. XV, 1929. Voir $ 12. 
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4. ONDES ET RAYONS GRAVIFIQUES. — Choisissons de nouvelles 
variables x, %, X:, %, telles que les nouveaux g,8 satisfassent aux quatre 
équations complémentaires (1) 


(14) ET ]= 0 : 
ou, ce qui revient au même, 


(15) g 26 — 0. 


En utilisant ces nouvelles variables, les équations fondamentales 
du champ gravifique deviennent 


(16) VAE = (5, ©) 


où les seconds membres (5,7) ne contiennent pas de dérivées secondes 
des potentiels einsteiniens. Remarquons de plus que chacun des 
premiers membres de (16) ne contient que des dérivées secondes 
d’un seul potentiel einsteinien. 

Nous dirons, d’après HADAMARD et VESSIOT, que les solutions 
He JU 2 7, vide 


(x7) 


fournissent des variétés caractéristiques ou ondes gravifiques (élémen- 
taires) : 
(x8) f=0, 

Considérons maintenant les caractéristiques (de Caucxv) de G = 0 ; 
on aura ainsi les bicaractéristiques ou les rayons gravifiques. VESSIOT 


a montré (!) que ce sont des géodésiques de longueur nulle du champ 
gravifique. 


5. CHAMP GRAVIFIQUE MASSIQUE. — Dans le cas d’un champ gra- 
vifique dû à des masses, nous posons 


(19) I = — g'É(TOu, ue + &,a) 


(x) Th. DE DonDEr. La Gravifique einsteinienne. Annales de l'Observatoire R. de Belgique, r921 
(ou chez Gauthier-Villars, Paris), Voir $ 20. . 


ce 
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où 1 est la densité massique généralisée, #, les composantes cova- 
riantes de la vitesse et 42 les tensions massiques. 
En utilisant (3) et (13) on obtient le tenseur 


8 le 2 
(20) GE = Vou,u' + LS. 


Le théorème de l'impulsion et de l'énergie (12) devient alors 


(21) PF, =, +4, =0 
où 
(22) M, = A, + 0, (Teuf) 
et 

a EL T7 = 
(23) ion D Beraté 


En multipliant (21) par w* et en sommant on obtient l'équation 
de continuité : 


(24) — (avu*) + CNET 


x 


6. CHAMP GRAVIFIQUE ÉLECTROMAGNÉTIQUE ET MASSIQUE. — 
Considérons maintenant le cas où le champ gravifique est dû à des 
charges électriques et à des masses. 

A cet effet, nous introduisons la fonction caractéristique : 


(25) M = 8 | LÉ nt en 4 S'HH;; |. 


Les fonctions H4x sont les composantes de la force électromagné- 
tique ; elles sont antisymétriques, c’est-à-dire que 


(26) Hi = — H;,. 
Grâce à (3) et (13), on aura le tenseur 
en Gén +el+ eV EH, HV + VE HE}. 


I1 résulte immédiatement de là, le {héorème de l'impulsion et de 
l'énergie 
(28) F, =, +2, +3 —0 


NT —— 
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où %, et 2, ont été définis précédemment (22 et 23) et où 


Ha. ti) 
ai = Nestes ip eee? 
(29) F, es [v— FLE E7A H,; dX; ; 

La notation H; signifie qu’il faut affubler H de deux indices for- 
mant avec ?j une permutation paire ijij des indices 1234. En mul- 
tipliant (28) par w2 et en sommant par rapport à >, on obtient 
l'équation de continuité (1) : 

à NE 
(30) (ou?) ee (a, AS J# — 0. 


CEA 


En multipliant (28) par A% et en sommant par rapport à +, on 
obtient la relation (1) 


a 


— (sf +e,) 
A A° 


(31) = 


7. ÉQUATIONS MAXWELLIENNES. — Introduisons maintenant les 
équations électromagnétiques maxwelliennes : 


FRET RE 
(32) \ Re 
236% 
(33) le 


où nous avons posé 


(34 HR =H; 
De (33) et de (34), il résulte immédiatement que 
o®P, Ps 
(35) Homes 


où ® est le potentiel électromagnétique. 
Il est important de noter que les équations de MAXWELL (32) 
peuvent être tirées d’une fonction fondamentale 


(36) De) = ob, + — £ gÉeŸH, He; 


(x) Th. DE DONDER. Théorie des Champs gravifiques. Mémorial des Sciences mathématiques. 
Fasc. XIV: Paris 1926 (Voir chapitres II et Vi). 


L 


O2, — 
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en prenant la dérivée variationnelle par rapport à d et en supposant 
que les H,8 qui y figurent ont la forme (35). Le symbole 5 est un 
facteur de densité électrique. 

Il résulte immédiatement de (32) que l’on a 


(37) à (ou%) = 0 


qui exprime la conservation de la charge électrique en mouvement. 
En tenant compte des équations (32) et (33), dans la valeur (20) 


de gl), on trouve 


(38) 4 = + SU, 
En multipliant ge) par u2 et en sommant, on obtient 
(39) gl = 0. 


Introduisons la valeur de gl?) dans (31) ; nous obtenons ainsi 
une relation entre les facteurs de densité 1% et & ; à savoir (1). 


de — 
(40) . 


On appelle, par définition, équation complémentaire de MAXWELL, 
la relation suivante entre les potentiels électromagnétiques 


PR be) 
où nous avons écrit x — gas. Les équations électromagnétiques 
de MAxwELL (32) peuvent alors être simplifiées. Après quelques 
calculs on obtient 
cu* " 2°®, 


où K, ne contient pas de dérivées secondes de ®,, ®,, D, ®D,. 


8. FORME LAGRANGIENNE ET FORME CANONIQUE DU THÉORÈME 
DE L'IMPULSION ET DE L'ÉNERGIE. — Considérons le cas d’une masse 


G) Loc. cit, éq. (341). 
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incohérente, c’est-à-dire où les &ag sont nuls. Le théorème de l'im- 
plusion et de l'énergie (28), si l’on tient compte de (38), peut s'écrire 


(43) F,—=% [ECET) — ls :[4(2) -()] Le 


où l’on a posé 


(44) WE = g,guuf 
(45) U= ue, 
Posons (1) 
(46) dx) = Mix x OX ÔX, ; drle) = x,0%,0430%. 


11 résulte alors des équations de continuité (30) et (37), en tenant 
compte de (39) et de ce que £:3 = 0, que 


d (mm) — . d a — 
(47) sf" l=0; Fe art 0 


En multipliant les deux derniers membres de (43) par un élément 
de volume et en intégrant, on obtient 


en Jo LÉCET) ET fee é() 5e 


La mécanique ondulatoire de de Broglie-Schrodinger 


9. MÉCANIQUE RELATIVISTE DES CHARGES PONCTUELLES. — Pour 
étudier la dynamique des charges ponctuelles dans l’espace-temps, 
nous introduisons grâce à (46), les deux constantes 7") et 7 qui 
caractérisent la particule au point de vue masse et charge. 

Nous écrivons alors l'équation (48) sous la forme 


d (= —(#)= 
ds \ou* PRE | 


(50) L = Y4W° + Us 


(49) 


où nous avons posé 


(x) Loc. cit, équations (188) et (184). 


» 
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et 


<(m) 


à — 
(51) = 


Introduisons les variables canoniques 


(52) ‘1 — 2 = or cp, 
ou 

et la fonction hamiltonienne 

(53) H=—L + pu. 


Il est facile de calculer la valeur de cette fonction, grâce à (52). 
On trouve 


(54) H = + 
Les équations (49) jointes à (52) sont alors équivalentes au système 


dE, ion 
\55 ds — op,’ C7 NANTES 


La fonction hamiltonienne (53) peut être exprimée en fonction 
des variables canoniques D, grâce à (52) ; on obtient 


(56) Hz gps — .) (pe — cb). 


Proposons-nous maintenant de trouver l'équation de JaACcoBr du 
système (55) ; posons, à cet effet, 


(57) Pa —= 2x. 
D'où, en vertu de (56), l’équahion jacobrenne : 
- CS Te 2S l 

(58) P+ret(e—e.) (Es) = 0 


où S est la fonction de JACOBI. La théorie classique de JACOBI nous 
apprend d'autre part que 


(59) donc ee DA 
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En substituant dans (58), on obtient 


3 [dS 2$ 
(60) g*e le — #,) re — #,) — I = ©, 
; 8 ; 


Remarquons enfin que la fonction de JACOBI peut s’écrire, grâce 


à (59), 


(61) S = — : Sie SX; La, %3, Hyde 


10. ÉQUATION RELATIVISTE DE LA MÉCANIQUE ONDULATOIRE. — 
Posons 


(62) ÊSE log y 


où k est une constante universelle dont la valeur sera donnée plus 
loin. On aura donc 


Se: 2 oW 

(63) F—— ess 
et 

ni: 

2S I dX, 

dS 


Reprenons alors l’équation jacobienne (60). Celle-ci peut donc 
s’écrire 


2 CR Oo oW\ /r oW ow ow'\? 
— 00 ae ep : LR et EE — 
(CSST Er Ë x, nee a Ç 2% Ron ) ei % 


Pour obtenir l'équation relativiste de la mécanique ondulatoire, 
proposons-nous d’extrémer l'expression (65) ; écrivons donc 


(66) AVIS ee, 


NI 
\ [0] 5 s * PAL . . 
où le symbole ;;; signifie qu'il faut prendre la dérivée variationnelle 


de Jy— g par rapport à W, à savoir 


4 
(67) Rasa; \ dif vo 
7 WoW cd ax, à oÙ | TEEN. Vu, Vo, X3, Xe 
= 0) LE En 


[ 


ER M 
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L'équation (66) explicitée donne l'équation de DE BROGLIE-SCHRo- 
DINGER généralisée (1) 


y ; 
CO Y — 2k: De. nt + (ke)? (F—5)v — 0 | 
(4 


où l’on a posé 
ee PONT 
GO Es (v eÙe sx) 
(70) le De qu DEA 
V— £ = 2x, +) 


(71) F = Diese 


(68) 


MEL N2, 314 


La constante universelle À est donnée par 


MC? 


(72) es k == 21T bc 


où m, est la masse de l’électron au repos. On aura, en vertu de (72), 


(73) e = — 2,42.10 ° À. 

Dans le second membre de (73) figure la longueur d'onde correspon- 
dant à la transformation du contenu énergétique mc? d’un électron 
en un quantum lumineux. 

Champ de Minkowski. — Supposons que le champ gravifique 
soit de MINKOWSKI et que les composantes ®,. P,, ®, soient nulles. 


Posons, en outre 
(74) dev Cid'on "ph 


Le système considéré étant stationnaire, posons avec O. KTEIN, 


(75) væ=ee ! 


où ® ne dépend plus de # et où E est une constante. 


(x) Th. DE DoNDER. Bull. Ac. Roy. Belgique. CI. des Sc. 5° série, t. XIIT. Séance du 5 février 1927. 
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Substituons dans l’équation (68) ; d’où, en vertu de (74), 


2 
(76) AY + À UE — eV) + 2(6 — eV)m cr = 0 
où l’on a posé 
(77) 8—E — mc 
Approximativement, l'équation précédente pourra s’écrire 
2 
(78) (EE ee Moc (6 — EV)Y — 0 
où 
w—N d°Y 
(79) À _ dX 2" 
t 


Cette dernière équation est l'équation ordinaire de DE BROGLIE- 
SCHRODINGER. 


Electrodynamique des corps en mouvement 


11. EQUATIONS GRAVIFIQUES. — Nous avons vu que les 10 équations 
fondamentales du champ gravifique pouvaient se déduire du principe 
variationnel (1). Nous poserons ici pour la fonction 9 : 


(80) IC — — VW? + am l”) + am) RS 


où T est le facteur tensoriel massique, où 
ve 8 
(81) W? == ag” ui ml 


et où ru”), anis}, Dire ‘el représentent respectivement :les fonctions 
caractéristiques des phénomènes massiques, massiques-électromagné- 
hiques et purement électromagnétiques. Nous dirons aussi que rlmse) 
caractérise les phénomènes d’électromagnétostriction ou plus sim- 


plement la sérichion. 


12. HQUATIONS ÉLECTROMAGNÉTIQUES. — Nous écrivons les équa- 
tions du champ maxwellien généralisées sous la forme 
dn*® 
82 “y — — 
( ) axe ——*\ 
dx"? 
8 NES ox 
(83) €, 
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où 
(84) CE AL + ee 
(85) Cu = oqujé + Lu) 


où l'indice (e) signifie électrique, et l'indice (u), magnétique. Les 
expressions o#% et 4 représentent respectivement les composantes 
des courants de convection et de conduction généralisés. 

Dans le cas le plus général, nous définissons la force électromagné- 
ligue par le tenseur antisymétrique 


(86) RT = HT —e(, 


et la force électromagnétique adjointe par le tenseur antisymétrique 


TE 
(87) Han 
où 

8 
(88) HR — 3£ 


Ecrivons maintenant avec FINSTEIN, la force de polarisation définie 


par les six composantes ei sous la forme 


CAC) (6 GRACE 


où {e) sont les quatre composantes contravariantes tensorielles de 


dx 
l'intensité de polarisation électrique. Rappelons que #4 — Fr - 


Exprimons de même la force de polarisation magnétique au moyen 
de l'intensité de polarisation magnétique par 


(90) 2ù = Lu — L{uy# : 
De (82) et (83) on tire immédiatement 


DL Mas 
(91) GOT dr 


Afin de pouvoir donner facilement l'interprétation physique des 
termes ci-dessus, supposons que %,, %, x représentent des coordonnées 
rectangulaires dextrogyres, et que x, représente le temps {. Au lieu 
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de x,, %, %, nous emploierons aussi la notation x, y, :. Il est alors 
commode d'employer les notations habituelles de l'électromagnètisme, 
en posant : 


: à _ 
KO — cK, \* =, 
st * 
(Ke — ch, x, = &, 
RE = x a 
(92) u 

Ré ——8. Ex 

KR ——B, LE = 

Rt=—8. =R 


Dans ces tableaux, les symboles ont la signification physique 
suivante : 


(x, Fy, 33) — composantes de la force magnétique 
(Rs, By, Ps) — composantes de l'induction magnétique, 
D'autre part, on a posé 
Où H=H-H; M=H-H 1 
où 


. (fe, Hy, Hi) = composantes de la force électrique 
(H°, H°, H°) — composantes de la force électrique appliquée (einge- 
prägt). 
Enfin 
(Bs, By, B:) — composantes de l'induction électrique. 


Il est presque inutile d'ajouter que toutes ces expressions doivent 
être entendues dans un sens généralisé. 

Grâce aux notations (92), les équations électromagnètiques (82) 
et (83) des corps en mouvement, conservent la forme maxwellienne : 


(94) div B— Ste + %) rot à = | 2 + set + ta) 
(05) div 8 — ju +4) rot H° — 1 (—Ÿ Ta ui) 
Les composantes rectangulaires du vecteur ju + % sont 
— qù — 
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CAPE Lin À = 1, 2, 3). er 
Te + An). 


13. Rerour AUx ÉQUArIONS GRAVIFIQUES. — Prenons la dérivée 
De nn de hs fonction, (60) par rapport à g,: et posons 


u 2) 2H) 291 (9) 
—_. (22 D = + 
(#6) ap — FA 27 = = APE ge 
Nous aurons donc 
(97) Can, + 67) + CS + C0, 


Dans le cas d’un ‘champ SE quelconque, nous 
posons (1) 


(98) 
| Il en résulte que 

Oo) = VE PRAT +R RT) TE nt 
Posons ensuite 


(100) Ps + 17; 
d’où 
(x0x) Le = Nu, + Pis + TS 


_ Il est facile de déduire de (99), (100) et (101) la valeur des compo- 
santes mixtes GA), PÉ, TÉ Le théorème de l'impulsion et de l'énergie 
(12) devient ici 
(102) | $, =6,;=R, +2, + 50) 
où M, et 2, ont été définis respectivement en (22) et (23) en tenant 
compte de (100). On a, en outre, 
mono =if=— Ale (RAT +K x) 


+= VE eg Æ. (k,xP)]: 


@) Th. De Donnez. The mathematical Theory oj Rélatvity (Massachusetts Institute of Techno- 
ogg: Carnbridge, Mass, US. A., 1927). Voir page 77- 


— 91 — 
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14. HYSTÉRÈSE ÉLECTROMAGNÉTIQUE. — Par définition, l’hystérèse 
électromagnétique est un quadri-vecteur #{°) (x — 1, 2,3, 4) donné 


par () 


(104) F6, —=4,—3 
où l’on a posé 
(105) RER ES à For e 


ou, en vertu de (82) et (83), 


Fa me y en di 
“2 as 6 as er 
(106) 4, = RTE TR LE 


On voit que l'expression (105) est identique à (29) dans les systèmes 
dépourvus de polarisations électrique et magnétique de sorte que dans 
ce cas 3e{e) = 0. Il en résulte que l’hystérèse, telle que nous l’avons 
définie, est essentiellement due à ces polarisations. 

En introduisant les notations (92) dans (106), il est facile de voir 
que les trois premières composantes de (106) peuvent se mettre sous 
forme vectorielle ; à savoir 


4 
(og) Sa (4 = [6-8] + met + [6-8] + ce”. 
7e) 


Cette expression généralise l'expression classique de la force de T,ORENTZ 
(multipliée par c). 

La quatrième composante peut s’écrire 
(108) L, = — c(H*.€) + (6.0, 


ce qui est l'expression générale de l'effet JOULE. 


15. HQUATION DE CONTINUITÉ. — Cette équation s'obtient en 
multipliant les composantes ge) par x et en sommant. On trouve 
( ES 
(109) sue ) De (2, + ge) }u* = 0 


x 


(x) Th. DE Donper. C. K. Ac. des Sc. de Paris, 2 juillet 1928. 
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ou encore 
(x10) . [AVS(x + :x6)] — — (2, Fe F0) jué*(a x), 
Il est facile de transformer (109) et de le mettre sous la forme : 


(III) d\reu”) 
dx” 


+ (ot) + cite H Ge) }u” #0 
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Les expressions entre parenthèses peuvent s'intégrer par parties : 
on aura par exemple 


: DA gétre),,* 
(112) (AVS — AE) — le) 
où l’on a posé 
(me) 2 I felme) | gaime)\ (4u* | duf 
en ct sed (c: CR ) dx; & CEA 


DORE Ke rs gare) du” _. du MS de: _ | 
he x dr, dx,] — 2 ds °me) 
On aura des expressions analogues pour H(*) et je), 
Grâce aux formules précédentes, l'équation de continuité peut donc 
enfin s’écrire 


(I14) ss À — Er + Be . A) 
où 
(115) 1 NTI ER 


16. PRINCIPE FONDAMENTAL DE L'ÉLECTROMAGNÉTOSTRICTION. — 
Nous dirons que Ale) défini par (113) est la puissance (par unité 
de volume) de la striction. En généralisant une hypothèse de la théorie 
classique de l’électromagnétostriction, nous admettrons (!) qu’on à 


(116) ge) = lu 

Autrement dit : La puissance de striction est égale à la punssance 
d'hystérèse. 

(x) Loc. cil., ci-dessus. 


ANNALES DE L'INSTITUT H. POINCARÉ. 7 
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Le principe (116) peut encore s’écrire 
(117) ge) — (self "È el) ur, = ee (PES, 2). 


; or us Pt 
En assimulant w* à =, on aura, en première approximation, 


en Eee) etre) 
(118) Kim) = (Hit + H }. 


D'autre part, grâce aux notations (92), on peut écrire Ho défini 
par (104), sous la forme : 


eo aff 9) (54) (nt) (2% ] 
+ © div. ([H*-%] — [B-R)). 


Dans le cas de déformations oscillantes du corps considéré, nous 
admettrons qu’on ait le bilan : 


(120) Hi*/o + > div. ([H*-%] — [B-#]) — o. 


Nous pourrons alors écrire finalement le principe fondamental (116) 
sous la forme : 


fe een x, 0B) _ {R,.H* ( 2 @ . 2% 
a) K°* = Z|(H ) (8 )+ KT —(s.)] 


17. CALCUL DE K{#%) EN FONCTION DES DÉFORMATIONS. — Posons 
(122) Hi 4% + (af, 44, «5, 7) 1=1,4,3 


où À représente le déplacement infiniment petit à partir du point 

ARS RER k sut 

initial x; (? — 1, 2, 3) ; pour { — O0, on a xi — x°. Le corps considéré 

dans son ensemble est donc au repos par rapport au trièdre (x,, x, #3) 
Il en résulte que 


: d}; 
I2 i— +, 
(123) =; 


Posons 


(128) 
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Adoptons aussi les notations classiques 


2 hu DAS : PRE dh3 . 
bar RMC 7e Carr 
PR ON En ons MR Di on Dhs  dh 
(34) sie LEE oy° ox? un RAT Qi Th En dy° 
o?. oh, ds 
X, — FES AC — 2 - Z,— —Ÿ 
dois ES RT 
et 
(125) RP LR t pet re SRE NE 
REC TRE DAT ur 
Ii est bon de se rappeler qu’à un infiniment petit près, on peut 
identifier x? avec x (1 — 1, 2, 3) ; nous pouvons donc, dans (8o), 


remplacer les x} par les x:. Il en résulte que les expressions qui figurent 
dans (113) deviennent 


/ dv! 2h, 2 dv? : dU° . 
Ze —= —= X;, y» ST = À 
ox dx d%2 dX3 
(126) | | 
| DOTE OU er p SD OUTRE. x DD 00e 
Ds LOT, A 7 MS AT EP ARRET ie 
EE 
: US GNT 2 dv! dv? ë DUR OT. 
Pr de LOS FREE RER 
(127) ds 2y oz ? y oz dx dX dy 


le point au-dessus des Xz, yy, … indiquant une dérivée partielle par 


rapport à £. 
La puissance de striction K(”*) peut donc s’écrire (113), dans cette 


I 
À : : CRT 
même approximation (x OC =) 


Ti %, + Témeÿ, + TS. 


ee er : : E ! 
he CU a nee ET m)E 
K (mc) —=—— 
— 5 (TROT En) (TE) TE On Jo, — (T4 Om — TI), 
JE Témeg, + Tiemeÿ, + TÉtm.02,. 


18. CALCUL DU SECOND MEMBRE DE (121). — Retournons au second 
membre de (121) et admettons qu’on ait 
H = DR F6; — > m'y 
E 


(729) 


ÿ 
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Posons 
(130) W= Tite), 


Donc W représente la densité de l'énergie électromagnétique l’oca- 
lisée dans le volume v 


(131) f Wèv = [(H-B) + (36.8) dv 
© U 
HS D'>euBB; Æ u/ 38 R ;) \d. 
2 y 
151 


Le second membre de (121) pourra donc s’écrire 


o 


me 


(132) 7 ea 


où les indices B,B servent à rappeler que la dérivée partielle par 
rapport à / se fait en laissant les Bs et les Bi (2 — I, 2, 3) constants ; 
on aura donc 


E à ‘oW NN I dE';; du'i; var 
(133) — El (Sr lee? = — — VN ES BB, + Den a, )d. 


Supposons que les coefficients £'i; et u' soient des fonctions des 
déformations linéaires et angulaires xz---2z, des rotations « et des 


vitesses v*. Autrement dit, la densité W de l'énergie électromagnétique 
dépend explicitement de 


(RER Ms Vrs 


O», Wy, W: Xe Ve Ze 


(134) 


Le second membre de (121) pourra donc s'écrire 


Fete (at nat 


s) 20] (+) e L FRE _. +3, g?* HS . < 

(13 _ dv = 
d /BR Jo Hé à + si a+ () À 
PE * ( LC d2 Je at 


ER 


da ex ds 


F 


4 L L = ñ == 
19, TEXSEUR DE L'ÉLECTROSTRICTION. 
: __{3W. 
| | «2j me F 
En - ram = (2 a 
Dr er : WW 
À TA — ME 
TAED + Tim) — 
$ orme (ET), 
on Jememo-(, 
| TM ot ls _ 
[2 (UYm0 —Tim0) =— Re 
" on F = NW 
a) D nr nr 
z Lio SM e) ee 
ci “ Ed 5e 
(627 | E) >” 
mis — (2W 
| | \2z/B2, 
mere 
| (140) Ti, Tim, Tim, Ti, 
__ Cependant, remarquons que l’on a d’une façon générale : 
(4) Ton = Tinp = DER = LE. 


Dans un champ de Mnsxowski, ces relations deviennent 


+ 7 ER jai 
a = Ts = — 77 — 


(x42) T9 = mA = & 
| ne 
Tn= a. 


point y#ms 


DO—X.2:3 


FAN 


(A 
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I1 résulte de (142) qu'ici les ji, sont symétriques (a, b = I, 2, 3); 
donc, en vertu de (138), la fonction .W ne renfermera pas 6%, y, 6) 


(me) ne diffèrent de 


explicitement. Il résulte de (142), que les TY 
r4{#) fournis par (139), que par le facteur (— c?). 

Dans un champ de MINKowSKI, nous avons donc calculé toutes 
les composantes du tenseur dat (x, BE %2) 407 seau pete ; on 
sait que er a les dimensions de l'énergie par unité de volume ; 
il y aura donc lieu de l’égaler à la densité de l'énergie élastique (de 
déformation). On pourrait aussi reprendre toutes ces théories dans 
l’espace-temps grâce aux corps parfaits définis dans notre Théorie 
des champs gravifiques (loc. cit. éq. 328). 


20. TENSIONS DE RADIATIONS. — ÆEnvisageons,. avec Léon 
BRILLOUIN (1), le cas d’une onde électromagnétique se propageant 
suivant l’axe des x ; soit Hy (ou Hy) le champ électrique, #, le champ 
magnétique. On aura, en vertu de (130) et de la définition d’une onde 
électromagnétique 


(143) — W—=B-H — 8.4. 


La vitesse de l’onde sera donc 
C 
(144) V = — 
1. 
où € et u sont respectivement le pouvoir inducteur spécifique et 


la perméabilité magnétique du système considéré. On a 


vI08D/ > “» 


(X45) sr Tale bib) 


où D représente la densité massique ; d’où, pour un corps isotrope 


oW 2W\ oW 
146 CE = ee ak 
(146) dXx Rs Cu lB 8 CE Fe 
(role = (ri) 
(Xe + Yy + 22) B,8 a n) log D B,8 


(x) Léon BRILLOUIN, Les tensions de Radiation ; leur interprétation en Mécanique classique et 
en relativité (J. de Ph., série 6, t. VI, n° 11, pp. 337-353). 


& 
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En vertu de (143) et (144), on obtient aisément 
(x47) A (Es) 


En formant grâce aux formules (136 à 139) et aux formules (142), 
le tenseur, à la fois électromagnétique et de striction 
(148) | mai(e) ae ra (me) | 


on obtient enfin le fenseur de 1. BRILLOUIN complété 


21. FORMULE DE HELMHOL1IZ-L1PPMAN (Champ de MINKOWSKI). — 
Considérons le cas de corps massiques parfaits, isotropes et au repos 
et supposons négligeables la variation de la densité massique dans le 
temps, ainsi que les fensions de MAxWELI. D'autre part, soit f(D), 
l'expression de la pression en fonction de la densité. En tenant compte 
des formules que nous avons trouvées pour le tenseur de l’électros- 
triction, on obtient très facilement 


2f(D) ù log € o log w 
Co CS), DE EN ne Vie D 
où 
(l I HP 
(57) Nec A se 


ce 
FE 
Nr log D é 5 He 2 x B,@ 
- o log V :>W 
: ie à log D ( 2: B @ 
. o log V pus 
(o) (0) — W.: | 
ologD 92: /B,8 
(v) 
1+(S 
M'A AU 2W un) 2W Tr Dette) 
vw] 2 an ae en B,8 ( oz, Es =LW ]. 
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Thermodynamique relativiste 


22. PREMIER PRINCIPE DE LA THÉRMODYNAMIQUE. — Nous écrivons 
le premier principe sous la forme (1) 
d Les 24 
(152) ds (auë(x, ue x) Se (x à Xi) Vu: Rè(x, #84) 


ou sous la forme équivalente 


DEUTREE" 
(153) _ pe 


Le symbole U est le facteur de densité d'énergie 2nterne du système ; 
2 est un facteur d’apports calorifiques à ce système. Un élément de 
l’espace-temps est représenté par d(x! - . . x). Le symbole X a été défini 
en (115). Voici les dimensions de ces symboles, en coordonnées car- 
tésiennes : 
95xdyozdt — énergie 
RÔxdyOzÔt — énergie 
UIxSyÈzd — énergie X longueur. 
Utilisons les coordonnées rectangulaires dextrogvres. La relation (x 


53) 
s’écrira 


a(UV-1*) » 
(154) mr Liu 2—R 
où 
dx? LES 
(155) DR "dE et V = di 
Multiplions les deux membres de (154) par dv défini par 
du — dxdydz; d’où en vertu de la théorie des invariants intégraux. 
(156) à [uv-è0 = f 250 — [ae 
di Tv Ca uv 


Posons 


(157) Vi fur= fuv-ce 
ral rAl 


(1) Th. DE DONDER. Comptes rendus de l’Ac. des Sc. de Paris, t. 186, pp. 1599-1601, 1928, 
et t. 187, pp. 28-30, 1928. 
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où V-1 représente la moyenne de V-1 prise dans tout le système consi- 
déré. 
Posons de même 


(158) Psp 
EIRE où HE 
dt n © € Jo 


en représentant par V la moyenne de V prise dans le système considéré. 
Alors l'équation (156) s’écrira 


(159) . [evU | = He face 


Nous dirons que U est l'énergie interne du corps à l'instant { et 
que dQ est la chaleur reçue par ce corps pendant d/f. 

Dans le cas de corps massiques parfaits, il est facile, en tenant 
compte de (115) d'écrire (159) sous la forme 


(160) DV = RE — HV ne + op 2 20 


où V, représente la composante suivant la demi-normale » extérieure 
aux corps en un point de l’élément do pris sur la surface fermée s 
limitant ce corps. 

Or, on a évidemment 


dv 
(161) pe, do — dt ° 


Nous retrouvons donc en première approximation, le premier 
principe de la thermodynamique classique : 


(162) | du — dQ — pôv. 


23. DEUXIÈME PRINCIPE DE LA THERMODYNAMIQUE. — En rela- 
tivité générale, nous écrirons ce second principe comme suit : 


(x64) AE) 


G(x'---x°) 
où % est le facteur de densité RS où 2* est le facteur (positif) 


TOI 
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d'apport de chaleur non compensée (ou de viscosité physico-chimique) ; 
où enfin T* est le scalaire (positif) /hermique. 
De (164), on déduit que ; 


o(gu) _2+2% 


(165) 0%, 1 gle 


Retournons à l’image euclidienne ; on aura, comme au (156), 


(166) 3 J sv-%0 — = Le + 9#)ôv 


où T* est la moyenne de T* dans le système considéré. 
Posons maintenant : 


(167) S — 4 SV— 130 
v 
et, comme au (158) 
(168) 7 9#3. 
Alors l'équation (166) deviendra 
dS _ c [dQ  dQ* 
(x69) di = ss (à A } 
Posons 
(170) T = T* ‘J 


C 


d’où l'équation différentielle exprimant le second principe de la ther- 
modynamique classique 

RS ar ON ET CNE 
(71) | RE Lee 0 


0 


24. THERMODYNAMIQUE DES SYSTÈMES  ÉLECTROMAGNÉTIQUES 
DOUÉS D'HYSTÉRÈSE ET ANIMÉS D'UN MOUVEMENT QUELCONQUE. — 
Utilisons le premier principe de la Thermodynamique en substituant 
(115) dans (152) ; d’où (1) 


d N / m AOC .(e) Ts 
(172) TS TUA(at et) 90 26) — [rl + gellut + xlO) out. à) 
(1) Addenda. Pour compléter cette théorie énergétique, il y aura lieu de poser 


YU — Ul”) + qe) + qe) 


et d'admettre que 
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Utilisons maintenant le second principe de la Thermodynamique en 
substituant (172) dans (164) ; d’où 


a7D RW 2) — TA À [Bee eat] = 
2Æ: al 7 Ru pe Rex. 2245) — DK O( xt. xt), 


On pourrait faire apparaître dans cette relation, la fonction 4 —T*y 
qui généraliserait l'énergie libre. 


Mécanique ondulatoire de Dirac généralisée 


25. HQUATIONS GRAVIFIQUES ET ÉLECTRONIQUES. — Présentons 
sous une forme tout à fait synthétique l'étude des systèmes électro- 
magnétiques les plus généraux que nous avons considéré au début 
du présent travail. À cet effet, rappelons d’abord les dix équations 
gravifiques ; à savoir 
PR ui ne 
où les fonctions caractéristiques 99, etc. écrites ci-dessus ont été 
définies précédemment. 

Passons maintenant aux quatre équations électroniques ou maxwel- 
liennes dont nous avons indiqué la généralisation dans (82) et (83), et 
supposons qu'il n’y ait pas de courant magnétique, c’est-à-dire que 
CZ — 0. Dans ce cas, on pourra écrire 


e af 


(175) ñ 


— — : 
. Pa,B — Pha 


Autrement dit, ces quantités dérivent d’un potentiel vecteur électro- 
magnétique D, - -. ®,. Nous allons aussi écrire ces équations maxwel- 
liennes sous forme de dérivées variationnelles ; à savoir : 


* * 


(or — aow2 + on f") + on (4) om) 
(176) ee TT NUE Lou 
@ 


—= ©, 


Précisons la manière dont la dérivée variationnelle par rapport à 
D doit être prise dans (176) pour qu’on obtienne les équations (82) 


nr 103 — 


TH. DE DONDER 


et (83) susmentionnées. Pour cela, reportons-nous à HU(*) dont la 
valeur est explicitement donnée en (98). Dans cette expression, nous 
considérons les K,; comme fonctions wniquement de x, x, (et non 
de ®,) ; les variations de ces fonctions par rapport aux ®, seront 
donc nulles. 

Retournons à (176). Nous savons que At?, W?, M (2) ne dépendent 


# 


pas des da. 
Posons 
AMC 
(177) i db, a Ce) 


où €, représente le courant électrique (total). 


Les dix équations gravifiques (174) peuvent alors s’écrire explicite- 
ment 


(178) _ : (a + Ogg + 0Gue = Nuits + Pis + ne 
et les quatre équations maxwelliennes s’écriront explicitement, grâce 
à (177), 


(179) ns Ce} 


26. FQUATIONS PHOTONIQUES. — Nous admettons que les équations 
du champ photonique sont de la même forme qjue les équations électro- 
niques (maxwelliennes) données par (82) et (83) ; posons donc pour 
ces équations photoniques (1) 


AUS. 
(180) dx: = Con) 
PF 
A Jo 
(x81) ax ia C(px) 


où UF représente la force photonique et su la force dualistique pho- 
tonique. Ces tenseurs antisymétriques sont donc les analogues de 


caf Ç 4 
h 


et 50e Les symboles C},,, et C3, Sont les courants fofaux dans 
le champ photonique considéré. 


(x) Th. DE DONDER. Le champ photonique Bull. Ac. Roy. Belg., 2 juin 1928, pp. 307-372. 
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Toutes les grandeurs qui figurent dans ce paragraphe seront consi- 
dérés comme étant complexes ; elles renferment donc chacune une 
partie réelle et une partie imaginaire pure. 

Pour établir un lien entre le champ photonique et le champ électro- 
nique (de Maxwell), posons (1) : 


(182) Co) = AP QU + D 
(189) on) = PEUT + 2. 


Donc les 4* et 1” définissent les courants photoniques diminués 
à P q 


des courants d'interaction entre le champ photonique et électrique. 
En substituant (182) et (183) dans (180) et (181), on aura une 
nouvelle forme des équations photoniques : 


| au°* 
(184) dx | d X° 
CAT 
où l’on a posé 
arts 4 : 
(186) [2 | —_— dx as «PQ. 
On a posé 
(187) Ps ci 


Par analogie avec la manière d'exprimer les forces électriques et 
magnétiques au moyen des potentiels ® électromagnétiques, posons 
ici 


(x88) Url] [01 
(189) 10, si [P;,] qi [Q...] 
où 


4P,  dE,T. dQ,  dQ, 
er) a, = 2 [Q]= Fe + al 


M 
où les crochets figurant dans les seconds membres de (190) doivent être 
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appliqués à chacun des termes qu'ils renferment ; ils ont la même signi- 
fication qu’en (186). Le symbole 0? signifie qu'il faut prendre la 
dualistique de Qu. D'autre part, les P, et Q, sont les potentiels photo- 


niques (1). 
On a posé 
___ 2ÈrMC 
(791) 1=— 7 


Toujours par analogie, avec le champ électronique maxwellien, 
écrivons les équations complémentaires photoniques ; à savoir 


ds" 
(192) [| = 
49° 
(193) [e. | = 48 


où # et # sont les porentiels de l’éther (?). 
Diminuons maintenant les courants photoniques 1* et 4° des gra- 


dients (contravariants) de ces potentiels $ et B de l’éther ; d’où 


DE — . ot as 

(194) PR VS a. | 
SR A Et 

(195) D = LV g gf A 


Nous dirons que 9” et 9° sont les courants photoniques proprement 
dits. 

Remarquons qu’en vertu de (180) et (181), on aura les théorèmes de 
conservation photonique. 


dC; 
(ph) 
(196) SON api 
FNAC 
x(ph) 
(197) es Ds 


() J. M. WuitrAkER. Proc. Roy. Soc., n° 788, p. 543 (3 décembre 192 }. 
(2) Th. DE DONDER, loc. cit, (2 juin 1928). 


» 
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27. PRINCIPE DE CORRESPONDANCE POUR LES SYSTÈMES ONDULA- 
TOIRES. — Introduisons maintenant la fonction caractéristique 


(198) ME + ue) on 
Si nous comparons cette fonction à celle qui figure dans (174), 
on voit qu’on a remplacé dans cette dernière tous les termes massiques 


[— AW? + AL é LE ar” d ] par le seul terme photonique 91 (24), Nous 


allons ane définir la fonction AK (1) au moyen des potentiels 
photoniques P, et Qu et des potentiels électromagnétiques d,. Nous 
écrivons, pour simplifier 


(199) Eu = *Pye 
Posons, avec J. M. Whittaker 
(200) NPA) = =» { DRUr, + 2(SS — BB) — 2(P Por — Q"Q,)} 


les tirets figurant dans (200) indiquant qu’il faut prendre l’imaginaire 
conjuguée de l'expression qui en est affectée. 

Prenons les dérivées variationnelles de I + Hu (°) + ot P#) par 
rapport à gx, à ®, et aux potentiels photoniques P:, Qu, PAU: 
Nous aurons alors les équations gravifiques de la mécanique ondulatoire : 


o(onf af ar (* e) Fa HUPI) 
se EE 


(201) 


les équations électroniques ou maxwelliennes de la mécanique ondulatoire : 


a(one + omlf) + one”) 
(202) Sr de 0 
d 


et enfin les équations photoniques de la mécanique ondulatoire : 


sm A ru(e) +R Ca gone + oncle) + on (#4) | 
Le LE — 0: PAPIER —9 
(203) [ES a £ èQ,, 
(om soie Pl) o(ons + on + me) 
(204) — 2% "—0; = So Pas = O: 
Ep u 
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En effectuant les calculs indiqués dans (201), nous obtenons expli- 
citement 


(205) — 2 (a + 0C)g + 0G = TES + M 


où TS est l’expression contravariante de (99) et où l’on a posé 


(206). SMS = np 0e E L—PÉP, + 00, )g* + 2N°%P — 20°? 


+ 2P*PÉ — 200. 
Dans cette expression les symboles Li, I, N5, Or’, sont définis 
par les relations suivantes : 


ais — — | [pp D = 
(207) RE as (Pal PR] + 1e [ni] 
(208) — gl L— C2 Pau 
(209) N° — ( Perf — Prest )( TE + conjuguée — : SSg** 


(210) Of — Ë Qbeee DE Ga de (= [a ]) + conjuguée — ; BB£g**, 


En effectuant les calculs indiqués dans (202) nous obtenons expli- 
citement : 


dr? ir eV— g 


5 (AC y = me : z 
(211) Fr Ts (u*r, + Ui'Qs + SP° — BQ° + conjuguée )- 


Enfin, en effectuant les calculs indiqués dans (203) et (204), nous 
obtenons les relations suivantes ainsi que les conjuguées : 


au* — sf dS Fr 
(212) lan = V— 8 g* [| TX 

aux dB 

REX Se LE CA SANS ET 
en [me [8]- 


En tdentifiant (!) les équations (178) et (205), (179) et (211), (184) 
avec les équations (212), (185) et (213), on obtient les équations qui 


(x) J'ai indiqué cette méthode d'identification pour trouver le principe correspondance. dans 
mes travaux antérieurs (voir ma note parue dans Bull, Ac. R. des Belg. CI. des Sciences, 5e série, 
t. XIII. Séance du 2 août 1927). 


» 
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expriment le principe de correspondance en mécanique ondulatoire 
Les équations (178) et (205) donnent 


(214) | Nutul -- PT — M°É 


où MxË est donné par (206) 
Les équations (179) et (211) donnent 


eV— £ = == = Re < À 
(215) | Co = — ae (ue, PAU IQ CSPRE= OS conjuguée) 


Enfin, les équations (170) et (211) ainsi que (184) et (212), donnent, 
lorsqu'on tient compte de (194) et (195), 


(216) Dé 797 


DES x g® 

On voit donc qu'il est possible d'exprimer le éenseur matériel, le 
courant électrique et les courants photoniques proprement dits en fonction 
des potentiels photoniques de Wluttaher. 


28. THÉORÈME DE L'IMPULSION ET DE L'ÉNERGIE ET THÉORÈME DE 
LA CONSERVATION DE L’'ÉLECTRICITÉ, EXPRIMÉS AU MOYEN DES POTEN- 
TIÉLS PHOTONIQUES. — Retournons au théorème de l'impulsion et 
de l'énergie (102) et remplaçons dans le tenseur T*f la partie massique 
NuuË + Pa par sa valeur Mt donnée en (206). On aura alors 


(217) (ur + Th). = 0 


qui exprime le théorème de l'impulsion et de l'énergie photonique. 
De même, retournons au théorème de la conservation de l’élec- 
tricité (91), et dans l’équation 


(218) 7 


remplaçons €{, par sa valeur (211). 


On aura le théorème de la conservation de l'électricité (exprimé au 
moyen des potentiels photoniques). 
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29. HOUATION DALEMBERTIENNE OU DE PROPAGATION DES POTEN- 
CIÉLS PHOTONIQUES. — En remplaçant dans (212) et (213), les U°Ÿ 
et og par leur valeur (188) et (189) on obtient les équations du second 


ordre : 
r ! oÙ niec V u r?m>?c" e? M D 
(219) g (P°). Te q De, a # h æ (P h de # h? (x ET mc? me }P: 
2Tiec un 
Fe HP, — H"Q, | — 0 
et 
dec .v! 2m?c? 2 . 
(220) ON) Se 66Q, pe TE æ"Q"), ne 1 D (z FE TGS 4,4) Q 


ainsi que des équations conjuguées. 


80. HQUATIONS ONDULATOIRES DE DIRAC. — Nous allons appliquer 
les équations générales ci-dessus au cas particulier d’un champ de 
Minkowski et nous choisirons les variables de telle façon que les ds? 
définissant celui-ci soit de la forme 


(221) ds? — dt — 3 (dx? + dy? + dr). 


Au lieu du potentiel vecteur ®, que nous avons utilisé ci-dessus, 
nous utiliserons avec Whittaker le potentiel ®’, lié au précédent par 
la relation 


(222) NE TER 7 


En utilisant les notations que nous avons définies précédemment, 
nous aurons 


I I I 
(223) P'=— = À, P=— e A Ps = — e A: 
WC 
Nous poserons de même 
db’ dd" 
22 X,,=— — + 
(224) D TRI ET 
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ce qui nous donne, en vertu de (175), 


TI dP D: mn FE 
(225) = —. ( SES = _ïye, 


Nous supposerons, dans ce paragraphe, qu’il n’y a ni polarisation 
électrique ni polarisation magnétique, de sorte que nous pourrons 
écrire 


3€ 
(226) X,, — — 2 H,,,. 


Remarquons que chez Whittaker, les indices y, », etc. varient de 
0 à 3 au lieu de 1 à 4. Pour plus d’uniformité nous utiliserons ici cette 
manière de procéder et choisirons 0, I, 2, 3 comme permutation fonda- 
mentale. 

Grâce à (221), nous pouvons déduire très facilement de (226) les 
composantes Hi ; à savoir 


(227) NE TH" 


Ecrivons, comme précédemment, 


X1 Er c? H,, = C sf X;: — c? dJG23 = — c 6 
I A ë T : 
(228) Re “= EH, en te Toys = — : 
T Be T 1 
nee XXE m5 Pia = — "5 4, 
Il en résulte, (227), que 
| KA H. X25 — 36, | 
(229) KP TH, XX — 06, 
XI =VCEH, KT 


Rappelons que H et # sont les forces électrique et magnétique. 
Posons encore, avec Whittaker 


(230) P,= P, Pi = — 
et | 


(231) Q=®% Q=—i9, M=— 9, Q=—:Q. 


TH. DE DONDER 


Nous aurons donc, grâce à (221) 


(232) P? = Je Pi = — Che P? = — CPS P—=— GE 
(25) LP QUMOE=S 0 nm QETc0,, = — Q.. 
Il nous reste enfin à nous occuper des U,. Nous poserons 
; | I 
Ur ss ut 
: I 
(234) US — c Sy ( U:: = ci by 
I | NO: 
US D c Se | ve — c 


ce qui donne immédiatement 


Ur cs, LOST A 
(235) Ua = — cs, ui — cb, 
US — — cs, LD 


Introduisons les notations (223) et (236) dans les équations (212), 
(213). Nous obtenons 


1110 ] 
_- [5] —[rot b] = mc P + [grad S] 
2rt 1 [os 
fiv = mere |] 


I [ob 271 
=. [3] + [rot s] = mc Q — [grad B] 


c L'oé 

Les crochets qui figurent dans ces équations ont la même signifi- 
cation que précédemment. 

Dans les équations ci-dessus, les symbôles P, Q désignent les vec- 
teurs ordinaires dont les composantes sont respectivement (Pz:, Py, P:), 
(Qr, Qv, Q:). Les scalaires $ et B et les vecteurs b(bz, by, ) et 
S(Sz, Sy, S:) sont définis par les relations (188) et (189) qui deviennent 
KE RE 


[div b] — 2e mc Q: — = (| 


27? 1[o L 
es mc b = — [rot P] — 2 | 2] — [grad Q:! 
27 LA 1 [oP 

TMS [rot Q] — : [5] — [grad P,] 


(237) | 
272 SE 1 [oP, 
h me S = [div P] + © [5e] 

2r1 net 1 [oQ, 

DE meB — [div Q] + = | | 


Are 
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À titre d'exemple, montrons ce que devient la première équation 
vectorielle (237) quand on effectue les opérations indiquées par les 
crochets ; nous obtenons 

I os 


(238) — de — rot b— grad $S + = — { i Vs—[A-0]—A. SE 


mc P. 


Les équations (236) ainsi explicitées peuvent être considérées comme 
les équations ondulatoires, les équations (237) donnant la valeur des 
vecteurs b et s, et des scalaires B et $S. Introduisons ces valeurs dans 
les équations (236). Nous obtenons : 


| DP 2 P]+ H-P,—[H.Q] — H.Q,} = 0 


2rie 


+ L(H:P) —(#.0)i — 


DQ Lee POI TO HP Ci 0 


DP, + 
(230) 


Lo (H-Q) + (4-P) | — 


où D est l’opérateur 


2? 2? SR NT on dre 
Gao Dessert [ 


=  . c? (V2 — A? — A3 — A?) 


Na +; + A, 


ù 
c ot + À . 


F 


Posons 
—P,+10Q;; V,—=P,+1P,; Vs——P,—1i0.; Y%,—=Q,—10, 
lo re D M 0e A 0h 
Grâce à ces notations les équations (230) sont équivalentes au sys- 
tème 
DY, + ci (286, + HOybs + 386.4, — Hubs + Hdi — Hiÿs) = 0 
Ds 2 (5004, — Dos — d06ba — His — Hu -+ Ha) = 0 
Ge) | 
DR = (96 Ya — Hoybs — 163% — Hs + 1H,% — H;4:) = 0 
Dh + _ Rs — Bb — 1H, — Eh, — Hibi + His) = 0 
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auquel il faut ajouter un système équivalent mais où l’on a remplacé 
les L par des w. Les équations (242) sont les équations du second 
ordre auxquelles satisfont les fonctions de Dirac d,, d,, d, d,. Si l’on 
pose 


a S,+1B; a——5,—15,, a——10,+$; a—0b,—1,0, 


CO TR RE Vi 


les équations (236) et (237) sont équivalentes à quatre systèmes de 
quatre équations ; le premier a la forme suivante 


PoYa + (Pa — Pa)hi + PaŸs = — mc a, 
] DoŸ2 + (21 + 1P2)Ÿs — Pay = — MC a 


40 | RS NE ee 
Por + (D, + 12), — Dao = mc a, 
où 
À 
(245) D=— Hi +IV;, pa; ete. 


Pour obtenir le 24 système, il suffit de permuter d et «x ; à savoir 


(246) À Por + (Pi — ÉPa)a, + Pros = — moy, 


Pour obtenir le 3e système, il suffit dans (244) de remplacer Ÿ par w 


et à par f, et enfin pour le 4m système, de permuter w et 6 dans le 
troisième. 


Il en résulte que si 4, de, ds, W est une solution des équations 
de Dirac, 


(247) au == Pa — y — LS M = I, 2, 3, 4 


constitue une solution de (236) et (237). 
Examinons ce que deviennent ici les composantes 


Ce = = Cle (œ = 0, I, 2, 3) 
du courant électrique (215). Nous poserons 
Aie Se . . 
(248) Gr Co = en = fs 
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En introduisant dans (215) les notations (223) à (235) , on obtient 
0 ——; (a + boat das dia +, B,—wBe +oB+o, Bi + conjuguée) 
(249) { = à (duos Voas + Vas + bas + &iB,+wBsoso +0 + conjuguée) 


nn nsc 


de sorte que pour une solution de la forme (247), on obtient 


p = — 2(db + db + WU + Wu) 


(250) je = 20 (bu us + dire ph) 


qui sont précisément les expressions trouvées par Darwin (:) confor- 
mément à la théorie de Dirac. 


(x) DARWIN. Proc. Roy. Soc. Série A., vol. 118, n° 780, p. 660, éq. (3. 3). 


— II5 — 


10. 


DE 
12. 
1e: 
14. 
I5. 
16. 
17. 
18. 
19. 
20. 
2%. 


22 
23. 
24 


25. 
26. 
CR 
28. 


29+ 
30. 


TABLE DES MATIÈRES 


La Gravifique einsteinienne 


Les équations fondamentales du champ gravifique................ Sete eee 
Léshidentités \gravitiques PA Rene RTE TO EU TE on 
Théorème de l'impulsion et de l'énergie............... 1... re 
Ondes setirayons sravinquese Mec-+-sue-meerre carie. ver. 
Champ gravifique massique...................... ao Eat Eire see 
Champ gravifique électromagnétique et massique............ HÉRE PHP O UT QUE 
Equations maxwelliennes........... et CR: see SAGE 
Forme lagrangienne et re canonique du théorème de l'impulsion et de l’énergie.... 


La Mécanique ondulatoire de de Broglie-Schrüdinger 


Mécanique relativiste des charges ponctuelles ...................... —- 
Equation relativiste de la mécanique ondulatoire............ Tape OC Ce 
Electrodynamique des corps en mouvement 
Equations gravifiques ............ ne ose re LL AR AO O0 nes os 
Equations électromagnétiques.......... AO EP en ee OR AUS 
RétOUR aa ONnS OAI TQUES ce ee eee Sosa ce 
Hystérèse électromagnétique ........4..... se en een Se ee sé ster.e 
quation de Conte ee. rc ee Te eee ss nero 
Pancipedtoudamentalidenlélectromagnétosinctont 2.7"... 
Calcul de K{(”») en fonction des déformations ........ RS none to 
Calcul dussecondmembrende (Cr) eee RERO a 
Tenseur de l’électrostriction ......, rues eue PE OL Ne OC vo Eu 
Hensions de adiatons een ee Peer ete re sas Me Eee SE 
Formule de Helmholtz-Lippmann .…........ 6... . s….. dette 
Thermodynamique relativiste 
Premier principe de la Thermodynamique......... RUE re te hnoe SE 
Deuxième principe de la Thermodynamique ............ SCT Re Suutere 


Thermodynamique des systèmes électromagnétiques doués d’hystérèse et animés d’un 


mouvement quelconque Sri On Fe Non nor 


Mécanique ondulatoire de Dirac généralisée 


Equations gravifiques et ho ORAN Ve dame SE UT Se 
Equations photoniques......... RAC Dre VA UT CR 5 ee 
Principe de correspondance pour es systèmes OndUIA tres enr e ASUS 


Théorème de l'impulsion et de l’énergie et Théorème de la conservation de lé électricité, 


exprimés au moyen des potentiels photoniques. . MT OC . 
Equation dalembertienne ou de propagation des potentiels photoniques.… . 
Equations ondulatoires de Dirac 


Sur quelques points de la théorie 


des probabilités 


PAR 


G. POLYA 


Les pages qui suivent reproduisent, abstraction faite de quelques 
altérations, les Leçons que j'ai eu l’honneur de faire à l’Institut 
Henri Poincaré en mars 1920. Je diviserai ici les matières en deux 
chapitres. Dans le premier, il s’agit d’une propriété caractéristique 
de la loi de Gauss et dans le second de quelques autres lois de fréquence. 
Ces dernières ressemblent à la loi de Gauss en ce qu’elles résultent, 
elles aussi , de la superposition d’un très grand nombre de petits effets 
fortuits, mais elles diffèrent de la loi de Gauss en ce que les effets 
fortuits dont elles résultent ne sont pas indépendants. Dans le premier 
chapitre, je donnerai les démonstrations en détail. Dans le second 
je m'efforcerai surtout de présenter les problèmes d’une manière 
intuitive et, quant aux démonstrations, je me bornerai quelquefois 
à une indication de la marche à suivre. 

Qu'il me soit permis d’exprimer ici tous mes remerciements à la 
direction de l’Institut Henri Poincaré et en particulier à M. Emile 
Borel de l’honneur de leur invitation qui m'a permis d'exposer mes 
remarques devant un auditoire d'élite. 
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I. — Sur une propriété caractéristique de la loi de Gauss 


I 


Envisageons des mesures pour lesquelles la probabilité d’une 
erreur ne dépend que de la valeur numérique de l’erreur. Désignons 
par w(x) la loi d’erreurs ; c’est-à-dire que la probabilité qu’une erreur 
tombe entre les limites 4 et b est donnée par 

b 


À o(x)x. 


a 


Supposons que # observations indépendantes d’une grandeur phy- 
sique dont la « vraie valeur » est / ont donné les valeurs 


li, b, le, sb be 


Les erreurs commises sont Z — /,, } — 1,, … 1 — 31 et la probabilité 
de les avoir rencontrées ensemble est proportionnelle à 


pU—h) gU—h)... el — 1h). 
En réalité on ne connaït pas / et nous devons trouver une valeur 
plausible de / en connaissant les valeurs /,, L, … HN fournies par l’obser- 
vation et en supposant connue la loi d'erreurs o(x). Il est naturel 


de prendre comme point de départ la distribution de probabilités 
dont la densité au point x est proportionnelle à 


(1) PU — hi) q(x — h) ++: px — ln). 

J'admets que nous étions d'accord jusqu'ici. Mais maintenant il 
s’agit de nous décider entre plusieurs chemins possibles. 

a) On peut choisir comme valeur plausible de la grandeur mesurée 
la valeur la plus probable, c'est-à-dire la valeur x — /’ qui rend (x) 
maximum. C’est le choix fait par GAUSS. 

b) On peut aussi choisir la valeur probable, c’est-à-dire la valeur /” 
déterminée par l'équation 

is 
| (x — 2") @(x — 1,) p(x — hs) . +. o(x — l,)dx = 0. 


u 


— L2 


PÉTTO 
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C’est ce que BERTRAND a fait remarquer (1) en critiquant le choix de 
GAUSS. 


c) On pourrait choisir aussi la valeur médiane, c’est-à-dire la valeur 
l" déterminée par 


le +00 
| EX — li) +. ox — 1,)dx — [ete — 1,) o(x — l,)dx. 
A té 


On pourrait encore choisir les abscisses d’autres points remar- 
quables de la fonction (r), mais arrêtons-nous ici. 

POINCARÉ, en critiquant la critique de BERTRAND, a remarqué que 
les deux chemins 4) et b) mènent au même but, sila loi d'erreurs o(x) 
est celle de GAUSS (?). C’est une propriété très commode de la loi de 
GAUSS : elle nous épargne une hésitation pénible.On est ainsi amené à 
se demander quelle est la loi d'erreurs la plus générale qui donne 

Lisa 
LL, lo, … Im étant quelconques ? 

Voilà la question à laquelle je tâcherai de répondre dans ce qui 

suit (?). 
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Nous demandons de la loi d'erreurs qu’elle nous fournisse le même 
résultat par le chemin a) et par le chemin b), quels que soient les nombres 
donnés Z, L, … Mn. C’est une chose bien compliquée, tâchons d’en 
démêler les postulats essentiels. 

Tout d’abord, la méthode a) doit être applicable. Nous posons donc 
comme premier postulat : 

(I) Etant donnés des nombres réels quelconques L, l, … ln, la fonction 
(x) atteint son maximum pour une seule valeur de x. 

Considérons maintenant #D mesures 


l,, RREN eS buse be cos Un, Un, in 


4 


où chacune des valeurs Z,, L, … ln est répétée p fois. En cherchant la 


(x) J. BERTRAND, Calcul des probabilités, 2e éd. (Paris, 1907). Voir p. 171-172. 

(2) H. Poincaré, Calcul des probabilités, 2° éd. (Paris, 1912). Voir p. 237-240 et aussi p. 176 
et les suivantes. 

(3) La réponse est formulée au n° vi. J'ai communiqué le résultat au Congrès international de 
Bologne, le 5 septembre 1928. Voir les Actes de ce Congrès. 
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valeur plausible pour cette combinaison de mesures, il faut d’après 
la méthode a) chercher le maximum du produit (1) élevé à la p-me puis- 
sance. La situation de ce maximum ne dépend pas de p. Donc le résultat 
fourni par la méthode b) doit être, lui aussi, indépendant de ? ; 
prenons cela comme second postulat : 

(II) Etant donnés des nombres réels quelconques L:, l», … ln, l'équation 


(2) fe — À) [e(x — h) q(x — h) :.. e(x — h)]Pdx = 0 


donne la même valeur pour L quel que soit l’entier positif p. 
J'ajoute encore les hypothèses supplémentaires que voici : 
(1') L'abcisse du maximum de o(x) divise l'axe des x en deux parties : 
dans chacune d'elles o(x) est monotone. 
(III) La fonction œ(x) possède une dérivée seconde continue, 
re 
o(x) > 0 et l'intégrale J xp(x)dx existe. 
=—00 
Le postulat que p(x) ne possède qu’un seul maximum absolu est 
contenu en (I), comme cas particulier ; (1') demande encore que o(x) 
ne possède qu’un seul maximum relatif. (III) demande en somme une 
certaine régularité de o(x). Ilest bien possible que les postulats (1') et 
(III) puissent être remplacés par d’autres moins restrictifs sans nuire 
au résultat final tandis que (I) et (II) me paraissent essentiels. 
Notre but est de chercher la fonction o(x) la plus générale satisfai- 
sant à (1), (1'), (II) et (II). 


III 
J'arrive à deux lemmes dont nous aurons besoin plus tard. 
Lemme À. — Supposons que f(x) et g(x) soient des fonctions continues, 
que g(Ëé) > 0 et que f(x) atteigne son maximum au point E et nulle part 


ailleurs, que f(x) Z 0 et f(x) tende vers o pour x > + et enfin que 
l'intégrale 


+00 
F tecaa1 sooax 
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soit convergente. Alors 


/ + co 


SEE 71 fac f{x)"dx — j(E). 


Ce qui est essentiel dans ce lemme est bien connu. Les hypothèses 
entraînent l'existence de trois nombres positifs e, à et f3, « < £, tels 
que 

I) SI 2e pour [x —i] PB, 
g(x) > 0 pour _ [x—#|<$, 
1 > JE) — ° pour [Xx—Ëêl<e 


On déduit le lemme A de l’inégalité double 


ta ER + 0 
eo" fl acadr — G@— 2 ( [+ [7 ) 1e6o 1e 
é—a — © +8 
+0 
= [ei fardx S 16 re frais ) 1 46) 
Lemme B. — Admettons que la fonchion f(x) est continue, que 


O< f(x) <'f(Ë) pour (x—EË) > 0 et que f(x) est croissante pour x < Ë 
ei décroissante pour x > £. Alors on peut conclure de 


(3) 1! XX PEr == 0 pour Del 29 


que f(x) est une fonction paire. 

Avant de commencer la démonstration observons que de la seule 
hypothèse (3) on ne peut pas conclure que f(— x) — f(x) comme 
le montre l'exemple de la fonction f(x) qui est égale à r dans les trois 
intervalles (— 8, — 5), (1, 4) et (5, 7) et égale à o partout ailleurs : 
elle satisfait à (3) sans être paire. 

L/hypothèse faite sur la forme de la courbe y — f(x) et l'existence 
de l'intégrale mentionnée dans (3) (D = 1 suffit ici) entraînent f(x) —0 


pour x—>» +, 


= PAL == 
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Si l’abscisse £ mentionnée dans le lemme B n’était pas 0, on pourrait 
déduire du lemme’ A que le premier membre de (3) a le signe de £ 
pour # suffisamment grand, ce qui contredirait l'hypothèse. Donc 


ma 


Ë — 0, le maximum de f(x) est /(0). 

Soient donnés les nombres y et , 
(4) OLE< y<}(0). 
Déterminons la fonction W(f) — W(f ; y, €) en posant 


O pour OZ=ÉZ<y—Ee 


Hs =} pour y <é< /(0) 


et en supposant que {V'({) est linéaire pour y — « < fé < y. D'après 
le théorème de WEIERSTRASS sur l’approximation des fonctions conti- 
nues on peut déterminer des polynomes ®,(é), &(#), … 


PBa(t) = Ani + Anot Se Anat° Sr <e PT 
tels qu’on ait uniformément pour o = # < f(0). 


(5) lim æ@,(t) — Y(i). 


n —> 
En utilisant l'hypothèse (3) pour p = 1, 2, … n on obtient 
+ co 


(6) J xf(x) Pn(f(x) )dx 


re na 
— f na f(X) + Anafl)? 2 2e + Gun (X)" dx = 0. 


Les polynomes fa(?) étant uniformément bornés pour 0 << f(0), 
l'intégrale (6) a une majorante de la forme 


+00 


D fixtio 


cn) 


où Z est une constante positive. On conclut ainsi de (5) et de (6) que 
+-00 


(8) J xf(x) W((x) 3 y, e)dx — 0. 


—@ 
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Cette dernière intégrale est encore majorée par (7) (avec / = 1) 
si on y fait tendre « vers o. Désignons par x, la plus petite et par x, 
la plus grande racine de l’équation 


f(x) = y. 


Ces racines existent d’après (4) et puisque la fonction continue 
f(x) tend vers o quand x —> + © ; on voit aussi que 


(9) OU 
D'après l'hypothèse faite sur la forme de la courbe y — j(x) 
(9 Zy pour A£<xr<£ 7% 
et f(x) < y si x est extérieur à l'intervalle (x,, x,). On obtient donc 
+0 X9 
Oo — lim ET Do, gr [rcredr at 
M à 
Tenant compte de (9), nous avons 
Ki = — X 


et puisque c’est vrai pour un y quelconque entre o et /(o) la fonction 
f(x) est paire. 


IV 


Le lemme B étant démontré, retournons à la recherche des fonctions 
p(x) satisfaisant aux postulats (1), (L'), (II) et (III). 

Observons qu’une substitution linéaire opérée sur x est sans in- 
fluence sur les postulats (1)-(IIT). En changeant x en x + a, où a est 
une constante facile à déterminer, on peut donc supposer que 


co 
(10) ie xo(x)dx = 0. 
D'après l'hypothèse fondamentale (II), on conclut de (10) que 
+00 
frere #10 pour Pi 12) 3, cer 
= 09 
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Ces équations, en nombre infini, entraînent en vertu du lemme B 
— dont les autres conditions sont remplies d’après (1') et (III) — que 
p(x) est paire 
(11) gx) = p(x). 

Soit donné un nombre c quelconque. Observons qu’en vertu de (11) 
la fonction 
(22) pt 6) px + oc) = px) pe) 

a: ue À (ne ue 


est paire. Elle n’atteint son maximum qu’une fois, d’après l'hypothèse 
fondamentale (I), elle doit donc l’atteindre au point x — 0. On a ainsi, 
en faisant usage encore une fois de (11) 


(13) gle— x) ec+a) <a} pour #>o. 
En posant 
C4 Sn USA, 
(x4) log &(x) — W(x) 


et en observant que x et c étant arbitraires x, et x, le sont aussi, nous 
pouvons exprimer (13) par 


a “HEC E) pour 424 
(15) exprime que la courbe y — W{(x) est partout concave vue d’en 
bas. 


Nous obtenons de (15) que, pour x, 2 %, 


L | Dre) + DV — 1) À PU rue 


VU=—I V—I V—=I 
ce qui exprime que la courbe 
(16) = W(X — hi) + W(x — de) + ee + W(x — ), 


elle aussi, est partout concave vue d’en bas. En particulier, nous 
obtenons ainsi, que cette courbe est croissante à gauche et décroissante 
à droite de son maximum unique. Observons que nous savons main- 
tenant davantage que nous n'avons postulé au commencement en 
posant (1). Nous savons maintenant assez de l’allure de la fonction (1) 


s 
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pour rendre applicable le lemme B à l'équation (2) valable pour 
P = 1, 2, 3, … Nous obtenons que la fonction (1) est une fonction 
paire de (x — ). 

En somme, nous avons obtenu la propriété suivante des fonctions 
o(x) satisfaisant aux postulats (1), (1'), (II) et (LIL) : Etant donnés 
des nombres réels quelconques L,, l, … l 11 existe un nombre 1, tel que 
D(x — À) p(x — 2) … (x — M) soit une fonction paire de (x — 1) et 
une fonction décroissante de (x — 1}. 

Si une loi d'erreurs o(x) jouit de cette propriété et Z,, L, … ln sont 
des mesures, elle fournit évidemment la même valeur avec les deux 
méthodes a) et b) expliquées au n° x : la courbe y — œ(x — /,) … 
(x — }h) est symétrique par rapport à l’axe vertical x — / qui contient 
aussi bien son maximum unique que son centre de gravité. Maintenant 
nous pouvons être assurés que les postulats (I), (II), (l') et (III) em- 
brassent toute la question du n° 1, ce qui n’était pas évident à priori. 

11 y a plus. Une loi d'erreurs w(x) jouissant de Îa dite propriété 
de symétrie fournira la même valeur aussi avec la méthode c) ou avec 
toute autre méthode raisonnable. 

Les déductions suivantes qui sont encore nécessaires pour déter- 
miner la forme des fonctions w(x) cherchées s’appuyeront sur la 
propriété de symétrie que nous venons d'établir. 


V 


Jusqu'ici nous n’avons pas fait usage de la première partie de l’hypo- 
thèse (III) concernant les dérivées de o(x). Il serait désirable de 
continuer le raisonnement sans faire appel à ces dérivées ; cela con- 
viendrait mieux au caractère approximatif de la théorie des proba- 
bilités qu’on ne devrait pas perdre des yeux en faisant des déductions, 
comme M. BOREL a si justement fait remarquer. Toutefois nous ferons 
usage des dérivées de (x) dans ce qui suit, parce que c’est plus 
commode. 

Occupons-nous de la fonction 


(x) — log p(x). 


En vertu des hypothèses réunies sous (III), F'(x) et W''(x) existent 
et sont continues. L'inégalité (15) entraîne, comme on voit facilement, 
que Y’’(x) est une fonction décroissante au sens stricte, le cas où W'(x) 
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est constante dans un intervalle est exclu comme le signe d’égalité 
dans (15). On conclut de (11), (14) que 


(7) View) 
et puis, W'(x) étant strictement monotone, 
(18) ME TO D OT LME TEE RO 


Soient donnés un nombre positif y et un entier "#1 supérieur à 1. 
Considérons la fonction de x 


PC 9) + ME (x) = A(x). 
En vertu de (17) et de (18) 
Ho) = W(—9)> 0, A(9) = (m — 1jV y) < 0. 
Donc il existe un » 
OSEO 

tel que Ar) — 0; donc 
(19) V(— y) + mT(r) = 0. 

C'est-à-dire que la dérivée de la fonction de x 
(20) g(X — y) (x + n)" 


s’annule pour x — 0. Mais puisque T''{x) est strictement monotone, 
la dérivée de (20) ne peut s’annuler qu’une fois. L'abscisse x — 0 
est donc pour la fonction (20) l’abscisse de symétrie dont l'existence 
a été établie au n° IV, la fonction (20) est paire. On a donc en faisant 
usage de (r1) 


(21) GX — Y) GX + n° = p— x — y) p(— x + n}" 
= px + y) p(x — n)", 


En éliminant m entre (19) et (21) et en utilisant (17) on obtient 


POS Pr) a Se 
Pat = 1) Sel 
ce qu'on peut écrire 
NACRE LA eee ’ 
(22) 2W/(y) P(x — ) a LA Ci n) a [w X + 4) — Viy— 4, P'ln) 
= D) — Va — y] 


Regardons y comme constante et faisons tendre l’entier #1 vers co 
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On conclut de (19) d’abord que (7) — 0, puis, W'(x) étant stricte- 
ment monotone, que » —> 0. En définitive (22) nous fournit pour 
nt —> 


(23) 20 (x) y) = [F(x + y) — Wix — y)]Y"(0). 
En différentiant (23) par rapport à x et à y nous obtenons 
(24) 2W°(x) W{y) = (x + y) + W(x — yjTE" (0). 


Observons que F'’(0) — o est exclu par (23) puisque W'{x) ne 
s’annule pas identiquement. En posant 


(25) ra = fU) 


nous obtenons finalement 


(26) 2/(%) {(y) = {x + y) + fx — y). 


C'est une des équations fonctionnelles classiques traitées par 
Caucxy (). Observons qu’en vertu de l'hypothèse (III) la fonction 
f(x) donnée par (25) est continue ; elle est paire, voir (17), et se réduit 
à z pour x — 0. Dans ces conditions la solution ne peut avoir d’après 
CAUCHY qu’une des trois formes suivantes : 


QT 
(27) ee = f(x) ee | cosh bx — : (er 2e Em) 
| cos x 


b étant une certaine constante positive. La troisième forme, cos bx 
est exclue par (17) et (18), la première n’est qu’un cas limite de la 
seconde, pour b —0. En observant que W''(0) doit être négative 
et W'(x) impaire, voir (17), (18), on obtient par deux intégrations 


(28) W(2) = 33 [1 — cosh (8x)] + log d, 
ce qui contient comme cas limite (b—0) 


(28*) W(x) = — _. + log 7; 


(a) A. L. CAucHY, Cours d'Analyse de l'Ecole Royale Polytechnique (Paris 1821). Voir p. 114- 
122, 
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b, c, d sont des constantes positives. En appliquant la transformation 
opérée au commencement du n° IV dans la direction inverse on trouve 
la forme définitive de o(x). 


VI 


Voici le résultat : 

Chaque fonction o(x) satisfaisant aux postulats (1), (1), (IT) et (III) 
est nécessairement d'une des deux formes suivantes 
(29) o(x) = decb"* [x — cosh b{x — a)] 

HO 
(29*) p(x) — de . 
b, c, d étant des constantes positives et a une constante réelle quelconque ; 
(29*) est le cas limite de (20) pour b —+ o. 

On peut ajouter qu'inversement 4, b,c étant donnés, b > 0, c> 0, 
a quelconque, on peut déterminer d de manière que o(x), donnée par (29) 
ou (29*) soit une loi d'erreur. Celte loi d'erreur fournira la même valeur 
pour la grandeur mesurée avec les trois méthodes a), b), c) expliquées 
au n° I, quelles que soient les mesures. 

En effet, o(x) étant positive, pour qu’elle puisse être admise comme 
loi d'erreur elle doit encore satisfaire à la condition 


+ 


IEC 27 


——00 


ce qui n’exige qu'un choix convenable de 4. 
Puis, Li, , … In étant donnés, l'équation en / 


(30) sinh b(Z — 1,) + sinh b(! — }) +... + sinh b(Z — 7) — 0 
aura une unique racine réelle, vu que la fonction 


i : REES e ? 


2 


sinh x — 


est, strictement monotone et varie de— 00 à + avec x. (On peut même 


dire que la racine Z sera comprise entre le plus petit et le plus grand 
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des nombres Z,, L, .… /n.) W(x) étant donnée par (28) — je prends, 
ce qui est sans importance, 4 — 0 dans (29) — on a, en vertu de (30) 


n 


Da W(x — 1) — n(cb? + log d) 
V—TI 
ñn 


— — cb? Ÿ° [cosh b{x — 1) cosh b(! — ,) + sinh b(x — J) sinh b(! — L,)] 


V—I 
ñn 


= — cb"? cosh (x — 1) Ÿ cosh b(! — 1). 
U—I 
C'est bien une fonction paire de (x — !) et une fonction décroissante 
de (x — 1}?. 
On peut vérifier encore plus facilement que la fonction (29*) possède 
la même propriété de symétrie ; dans ce cas-là, même 


e(a(x — 4)) e(a(x — h)) +: e(alx —h)) 


est une fonction paire de (x — J*) ; 4, 4, … an, li, L, … L étant quel- 
conques et /* étant choisi convenablement. (Cette propriété caractérise 
d’ailleurs complètement (29*) puisqu'elle ne convient pas à (29); 
cela a une certaine importance quand on envisage des mesures qui 
n’ont pas toutes la même précision.) 

En somme nous avons complètement déterminé les lois d’erreur 
qui mènent au même résultat avec les deux méthodes a) et b). Elles 
forment la famille de courbes (29) à 3 paramètres, contenant la courbe 
de GAUSS (29*) qui ne dépend que de 2 paramètres. En supposant (10) 
c’est-à-dire qu’il n’y ait pas d'erreur constante, on obtient 4 — 0, 
donc le nombre des paramètres se réduit d’une unité. 


VIT 


Je vais ajouter quelques remarques sur la loi d'erreurs cyclique, 
dont la notion est due à M. von Mises (5). Je n'ai pas à examiner 
ici si cette loi est applicable ou non au problème certainement pas- 
sionnant de la « Ganzzahligkeit » (problème des valeurs entières tirées 
de l'observation) que son auteur a en vue; je m'occuperai d’elle comme 


(5) Physikalische Zeitschrift (1919), p. 490-500. 
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d’une notion mathématique apte à éclaircir par analogie les questions 
correspondantes sur la loi d’erreurs ordinaire. 

Admettons qu’en mesurant des angles la probabilité d’une erreur 
ne dépend que de sa valeur numérique. (Les «angles » sont des valeurs 
réelles quelconques envisagées mod. 27.) Comme loi d'erreurs cyclique 
je désignerai une fonction o(x) telle que la probabilité qu'une erreur 
tombe entre les limites + et 5 est 

É 
f star (a LB< a+ 27). 


(o à 


On doit avoir 


(5 += emzo  femax 


Si n observations indépendantes d’un angle dont la vraie valeur 
est À ont donné les mesures À, A2, … Àn, la probabilité de les avoir 
rencontrées ensemble est proportionnelle à 


PA — hi) A — he) + PA — 2). 


Si on ne connait pas la vraie valeur de l’angle mesuré, seulement 
les mesures À, À, … Ar fournies par l'observation et la loi d'erreurs 
o(x), on peut se proposer de former une conjecture plausible sur la 
vraie valeur. Mais ici surgit une nouvelle difficulté, inhérente au 
nouveau problème. Il y a évidemment des cas où on ne peut pas 
former une conjecture raisonnable, comme par exemple lorsque ;, 
Aa, + An désignent les directions de.” rayons d’un cercle qui sont 
dirigés vers les sommets d’un polygone régulier de n côtés inscrit 
au cercle. Il faut éviter de demander trop. Bornons-nous au cas où 
les mesures couvrent moins qu'un demi-cercle, donc où on a 


(32) K=hSinrShe Mrs 


Nous n'exigeons des méthodes suivantes un résultat déterminé que dans 
le cas (32). 


ILest naturel d'envisager la distribution de probabilités sur le cercle 
unité dont la densité au point x est 


(33) QÉX — À) p(x — À) » + p(x — d), 


» 
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(Le rayon passant par le « point x » forme l'angle x avec une direction 
donnée). On peut proposer deux méthodes. 

a) On prend comme valeur plausible de l’angle cherché la valeur 
qui rend (33) maximum. C’est la méthode proposée par M. von Mises, 
analogue à celle de GAUSS. 

b) On cherche le centre de gravité d’une masse, répartie le long du 
cercle unité dont la densité au point x est (33) et on détermine l'angle 
cherché par le rayon qui passe par ce centre de gravité. Cela revient 
à prendre comme valeur plausible de l’angle cherché celle des deux 
valeurs À satisfaisant à l’équation 

2T 
[ sin (x — À) gx — À) pÜx — do) +: + pÜx — A)dx — 0 
© 
qui est située dans le plus petit des deux arcs déterminés par ?, et 
}n ; n'oublions pas que nous avons supposé (32). Cette méthode n’était 
pas proposée par M. von MISES mais est suggérée par ses considérations. 
Nous avons un problème analogue à celui du n° I si nous cherchons 
les Jois d'erreurs cycliques pour lesquelles les méthodes a) et b) mènent 
nécessairement au même résultat. Il sera suffisant d’esquisser la marche 
générale du raisonnement. 
Comme au n° II, nous pouvons détacher les 2 postulats suivants : 
(1). Quels que soient À4, 22, … An satisfaisant à (32), la fonction (33) 
n’a qu'un maximum absolu entre o et 2x. 
(II). Quels que soient À, À, … An, satisfaisant à (32), la solution de 
l'équation en À 
2T 
fs (x — à) [p(x — À) px — X,) 2: g(x — n)]Pdx = 0 
O 

est indépendante de D, p = I, 2, 3, … 

Pour pouvoir appliquer les méthodes des n°% III-V j'ai dû en 
outre faire les hypothèses supplémentaires suivantes : 

(I’’). La fonction (33) n’admet qu'un maximum yelatif entre o et 27, 
c’est-à-dire qu’elle est monotone entre son minimum et son maxi- 
mum, les deux étant univoquement déterminés mod. 2r. 

(II). Les dérivées p'(x) et o''(x) existent et sont continues : 
p(x) > 0. 
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Avec des modifications faciles, les raisonnements des n°5 ITI-V mènent 
à l’équation fonctionnelle (26) dont il faudra maintenant prendre la 
troisième solution, inutilisable pour le problème précédent (voir (27)). 
Plus précisément, vu le postulat (I) et la périodicité exigée par (31), 
il faudra la prendre avec b — 1. D’où on obtient la loi d'erreurs cyclique 


(34) p(x) = de£ cos x 


C’est la loi qui a été déduite par M. von Mises L. c. d’un postulat 
différent, analogue au postulat dont Gauss a déduit sa loi. II me 
paraît intéressant que cette loi cyclique (34) conserve avec la loi 
de GAUSS une analogie aussi au point de vue de la comparaison des 
méthodes a) et b). 

Beaucoup plus loin l’analogie n’ira pas, me semble-t-il. En cherchant 
une loi d'erreurs cyclique qui provient de la superposition d’une 
multitude de petites erreurs indépendantes, on doit arriver à la loi 
de GAUSS « enroulée autour du cercle » ; c’est la loi 


"+ 00 (x — 2nx)? 


—_— I LE 2C 
= e 


— —— 00 


représentée, ce qui est assez curieux, par une série théta. 


IL. — Sur quelques courbes de fréquence résultant de la superposition d'effets 
fortuits interdépendants 


£ 


La situation privilégiée de la courbe de Gauss dans la statistique 
mathématique s'explique par deux raisons, l’une expérimentale 
et l’autre théorique. L/expérience nous montre que, dans des domaines 
très divers, il y a des courbes statistiques qui, avec une approximation 
raisonnable, peuvent être assimilées à la courbe de Gauss. La théorie 
nous apprend que c’est la courbe de Gauss qui doit résulter de la 
superposition d’un très grand nombre de petits effets fortuits indé- 
pendants. La constatation théorique élucide le fait expérimental ; 
plusieurs grandeurs dont la statistique nous est familière, comme les 
erreurs d'observation, les dimensions des organes des êtres vivants, etc., 
subissent des modifications dues à l’action d’une multitude de petites 
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causes incontrôlables qu’on peut raisonnablement considérer comme 
« dues au hasard ». 

Mais dans tous les domaines d’application on rencontre aussi des 
courbes statistiques qui s’écartent considérablement de la courbe 
de GAUSS et une tâche très importante de la théorie est de chercher 
d’autres courbes théoriques, auxquelles les courbes expérimentales 
pourraient être raisonnablement assimilées. Ce problème paraît, par 
sa nature, indéterminé. On connaît les travaux de PEARSON, BRUNS, 
CHARLIER, KAPTEYN et d’autres qui en présentent des solutions 
très différentes, chaque solution généralisant la courbe de Gauss 
dans une autre direction. 

T1 me paraît qu’il ne faut pas considérer ce problème d’une manière 
formelle et chercher des familles de courbes quelconques qui se dé- 
duisent de la courbe de GAUSS par une opération analytique quelconque 
ou partagent avec la courbe de GAUSS une propriété analytique quel- 
conque. Ce qu’il faudrait élargir c’est la base théorique de la courbe 
de Gauss et c’est le chemin suivant qui semble se présenter le plus 
naturellement. Gardons la conception fondamentale que la courbe 
théorique cherchée résulte de la sommation d’un très grand nombre 
d'effets dûs au hasard, mais abandonnons l'hypothèse supplémentaire 
que ces effets sont indépendants. Le problème mathématique qui 
s'impose de cette façon peut être formulé ainsi : étudier les lois limites 
provenant de la somme d’un grand nombre de variables aléatoires 
interdépendantes. 

C’est ce problème vaste, trop vaste peut-être, que j'entreprendrai 
d’éclaircir par quelques exemples choisis avec soin. Ils nous feront 
entrevoir quelques conclusions générales et nous mèneront à quelques 
lois limites de fréquences, qui ne seront pas simplement des formules 
commodes pour l’interpolation, mais des formules susceptibles d’inter- 
prétation. 


II 


Considérons une suite de variables aléatoires (f) 
Xi, Xe, X3, °°" Xn, Xn+1, * 
et concevons-les comme réalisées par des expériences successives, 


(6) Les variables aléatoires seront désignées par des caractères gras. 
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x, faisant le commencement et x,}, venant après x,. Le cas général 
d’interdépendance consiste en ceci que Les chances de Xx,+, sont 
déterminées par l'ensemble des valeurs que les expériences précédentes 
ont attribuées à X;,, X3, .… Xy. 

Un cas particulier important est celui où les chances de x,+; 
ne dépendent que de la somme x, + x, + -.. + x,, deux combi- 
naisons X:, X,, … X, de même somme donnant les mêmes chances 
pour x,+,. Dans ce cas-ci je parlerai d'influence globale. 

Un autre cas particulier important est celui où les chances de 
X,+, ne dépendent que de xr, deux combinaisons x;, x,, . X, avec 
la même valeur de xr donnant les mêmes chances pour x,+,. Dans 
ce cas-là je parlerai de l'influence du prédécesseur. Ce cas a été intro- 
duit par A. MARKOFF qui en a fait une étude approfondie sous le nom 
d'épreuves liées en chaîne (?). 

On dit que les variables aléatoires x;, x,, … X,, X,41,.. sont 2dé- 
pendantes, si les chances de x,+, ne dépendent pas du tout du résultat 
des épreuves précédentes, deux combinaisons quelconques x;,, x, …. X, 
donnant les mêmes chances pour x,+1. 

L'interdépendance des variables aléatoires quelconques x,, xs, … 
X,, Xn+1, … Consiste en ceci que les chances de x1 sont modifiées 
si on fixe les valeurs d’autres variables x;, X4,... Xp, X4, .… qui peuvent 
en partie précéder x, et en partie la suivre. La distinction faite entre 
influence générale, influence globale, influence du prédécesseur et 
indépendance a été basée sur la considération des chances de x,+: 
quand toutes les variables précédentes x,, x,, … x, sont fixées et aucune 
des variables suivantes x,+:,, X,+9, … ne l’est. Il est encore important 
de considérer les chances de x: lorsqu’aucune variable n’est fixée, ni 
À, Ni X9, … Ni Xh+1, …, les chances de xn «avant la première épreuve » 
que je nommerai les chances à priori de x». (Une confusion avec 
la signification demi-philosophique du terme « probabilité à priori » 
me semble facile à éviter.) Si les variables aléatoires X,, x,, Xs, … sont 
indépendantes, les chances de x» après la fixation des valeurs d’autres 
variables quelconques X;, Xx, … Xp, Xq, … sont égales à ses chances 
à priori. 


(7) A. A. MARKOrF : 4) Wahrscheinlichkeitsrechnung (Leipzig-Berlin, 1912), p. 272-298 ; b) Sup- 
plément à la 3° édition russe de l'ouvrage précité ; traduction française parue séparément (St.-Pé- 
tersbourg, 1913), p. 44-66. 


L 
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Je dirai que la suite x,, x,, … x,, … est homogène si les chances à priori 
de toutes les variables sont les mêmes. 

Les variables x,, x,, .. X,, .. peuvent varier d’une manière continue 
ou discontinue, elles peuvent aussi être des vecteurs dans un espace 
d'un nombre quelconque de dimensions. 

Nous avons le cas le plus simple, lorsque chaque épreuve ne pré- 
sente qu'une alternative, c’est-à-dire lorsque x: n’est susceptible que 
de deux valeurs différentes, disons de o et de 1 (n—1,2,3,.). Dans 
ce cas je parlerai d'épreuves où de variables alternatives. 

Le cas le plus simple, le plus classique et aussi le plus important, 
est celui où les épreuves sont à la fois alternatives, homogènes et 
indépendantes. Chaque épreuve devant présenter une alternative, 
chacune la même alternative, avec des chances indépendantes du 
résultat des autres, l'épreuve en question peut être représentée par 
un tirage d’une urne qui ne contient que deux espèces de boules, 
disons rouges et noires, l’urne présentant à tous les tirages la même 
composition. Il s’agit donc de tirages réitérés d’une urne dont on main- 
tient la composition constante, en remettant chaque fois la boule 
tirée. La suite de tirages ne dépend que d’un paramètre, de la proba- 
bilité de tirer rouge de l’urne en question que nous désignerons par p. 
C’est le schéma d’urnes le plus classique à propos duquel Jacques 
BERNOULII a trouvé son théorème et DE MoivrE a rencontré la Loi 
de fréquences qu'on nomme généralement la loi de GAUSS. 

C’est le schéma d’urne que je veux modifier essentiellement, mais 
aussi simplement que possible. Je considérerai des épreuves alter- 
natives, homogènes mais non-indépendantes et je fixerai le degré 
d’interdépendance par un seul paramètre, quis’ajoutera au paramètre p, 
probabilité à priori de tirer rouge, la même pour toutes les épreuves. 
On pourrait ajouter ce second paramètre d’une infinité de manières, 
mais je ne considérerai que deux structures d’interdépendance, essen- 
tiellement différentes : je prendrai un exemple de l'influence globale 
et un autre de l'influence du prédécesseur. 


IIT 


Un phénomène familier où il s’agit d’une influence globale (avec 
des à peu près inévitables en statistique) est le progrès d’une épidémie. 
Nous jugeons les chances d'éviter l'épidémie surtout d’après le chiffre 
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total des cas de maladie ou de mort. Le fait élémentaire dans le progrès 
d’une épidémie est la production des germes ; celui qui tombe malade 
produit de nouveaux germes, en nombre beaucoup supérieur et aug- 
mente par cela les chances de son entourage de tomber malade. En 
réduisant ce fait à son expression la plus simple et en y ajoutant une 
certaine symétrie, propice au traitement mathématique, nous sommes 
amenés au schéma d’urnes suivant : 

Une urne contient originalement N boules, dont R sont rouges et 
S noires, R + S — N. Nous faisons de l’urne des tirages successifs 
en ajoutant à l’urne, après chaque tirage, à la place de la boule tirée 
1 + À boules de la même couleur. Si À est positif, le nombre des boules 
augmente après chaque tirage, chaque succès obtenu favorise les 
chances des succès à obtenir, chaque insuccès gâte encore les chances 
des épreuves suivantes, le succès ainsi que l’insuccès sont contagieux. 
Ajouter un nombre négatif signifie enlever ; donc si À est négatif, 
le nombre des boules diminue au cours des tirages, chaque succès 
obtenu diminue les chances de succès ultérieurs ; mais aussi les 
insuccès sont de la même nature et chaque coup tend à amener un 
revirement de fortune. Si À — o, les épreuves sont indépendantes ; 
c’est le cas classique de la boule remise. Le cas À — — r, celui de 
la boule non-remise est aussi classique. 

Associons aux tirages successifs les variables aléatoires x,, xo, X3, … 
le premier tirage correspondant à x,, le second à x,, etc. Soit Xn —=1 
si le tirage correspondant fait sortir une boule rouge et soit xs — 0 
s’il fait sortir une noire. 

Après n tirages l’urne contiendra un nombre total de N + »A 
boules, dont R + (x, + x, + :-.  x»)A rouges. La probabilité de 
tirer une boule rouge au (7 + 1)-me tirage est donc 

RS CE ee Exp )A 


(1) ENT PH (XI Xa + css + xn)d 
En mt À 


I+ nm 
Je pose ici comme dans ce qui suit 
R S Lt 
(2) Ro PNA T2 D PRERSS 


(x) est la probabilité de l'équation xnz1 — 1. Klle ne dépend que de 
la somme x; + x, + :: + x des variables aléatoires précédentes :; 
c'est bien le cas de l’influence globale. 

La probabilité (1) contient deux paramètres, » et d ; p est la pro- 


L1 
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babilité à priori de l’équation x, — 1, d exprime par son signe la 
direction et par sa valeur absolue l'intensité de l’interdépendance. 
On peut attribuer à 9 une valeur quelconque de l'intervalle 0 < p < 1. 
Etant donné o, si la série comprend » tirages, d doit satisfaire à la 
condition 


(3) PE = 1)0 — 0; 


le premier membre est le numérateur de la probabilité d’avoir Xn — 1 
APTÈS X, = Xe — +: — Xn-1 — I. La condition (3) est l’unique à laquelle 
à doit satisfaire pourvu que p < 0. 

La structure d’interdépendance que je viens d'expliquer et qui 
peut être caractérisée soit par le schéma d’urnes considéré, soit par 
la formule (1), sera nommée brièvement la structure de la « conta- 
gion » (S). 


IV 


Le cas le plus simple de l’influence du prédécesseur est représenté 
par le schéma d’urnes suivant. 

Pour faire une série de tirages, on dispose de trois urnes, U,, U,, U. 
chacune contenant des boules rouges et noires. De l’urne A, on ne 
fait qu’un tirage, le premier de la série. Si la boule sortie au n-me 
tirage était rouge, le (n + 1)-me tirage se fera de U,, si elle était 
noire, il se fera de We. On maintient constante la composition de 4, 
ainsi que celle de U en remettant chaque fois la boule puisée. Les 
probabilités sont données par le tableau suivant : 


urne KL TR PERS | 28 
probabilité de boule rouge CET 
probabilité de boule noire CHSCT 


Donc p, est la probabilité du « rouge après rouge », s, du « noir après 
rouge », etc. On a 


(4) (CP + Gi = I, PTT: Ps 0e = EL. 


(8) Voir F. EGGENBERGER et G. P6LYA : a) Zeitschrift f. angewandte Mathematik u. Mechanik, 
t. III (x923), p. 279-289 ; b) Comptes-rendus, t. 187 (1928), p. 870-872 ; F. EGGENBERGER, Thèse, 
Zurich, 1924 ; À. À. MARKOFF, Bulletin de l’Académie Imp. des Sciences, Pétrograd (1917) p. 177- 
186 (M. Serge BERNSTEIN a eu l’amabilité de me signaler cet ouvrage écrit en russe). Essen- 
tiellement la même structure de probabilité a été envisagée par Léon BRILLOUIN en connection 
avec un problème de la théorie des quanta, Annales de Physique, t. 7 (1927), p. 315-331. 
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Associons aux tirages successifs les variables aléatoires x;, x°, 43, . 
de la même manière qu’au n° III. La probabilité que la boule sortie 
au (n + 1)-me tirage soit rouge ou, ce qui est la même chose, la 


probabilité d’avoir Xn+1 — 1 est donnée par 
(5) Pa (fe Ps) en fe EX a 
où j'ai posé 
Po S 
GRR SN ; ‘ 
(6) É Vo Va & p le 


La probabilité (5) ne dépend que de xx ; c’est bien le cas de l'influence 
du prédécesseur. 

Si + est positif, la structure des tirages favorise la ressemblance 
de chaque épreuve à la précédente ; elle favorise la dissemblance 
si r est négatif. Dans le premier cas, la valeur absolue de 7 exprime 
la force de la « tradition », dans le second la force de L’ « opposition ». 

Un exemple aussi simple qu'ingénieux de cette structure d’inter- 
dépendance est due à MARKOFF lui-même (°). Il est familier à tout 
le monde : c’est un texte quelconque dans lequel les lettres, consonnes 
et voyelles, se suivent. Entendons par ps, ps, 5, 55 les quatre pro- 
babilités suivantes déterminées expérimentalement par dénombre- 
ment direct : « consonne après consonne », «consonne après voyelle », 
«voyelle après consonne », « voyelle après voyelle ». Nous aurons un 7 
négatif, donc le cas de fa dissemblance et le schéma d’urnes expliqué 
correspondra approximativement aux faits observables (1), On 
pourrait penser que le cas de 7 positif, celui de la ressemblance, 
se trouve réalisé dans une longue lignée de descendants successifs. 
Les dispositions sont transmises de père en fils ; les influences ances- 
trales de prédécesseurs très éloignés peuvent se faire sentir, mais 
seulement indirectement, par l'intermédiaire du prédécesseur immé- 
diat (M). Quoiqu'il en soit, simplement pour avoir une expression 
brève, je parlerai du cas discuté des épreuves liées en chaîne, carac- 
térisée par (5), comme de la structure de 1’ « hérédité ». 


(o) Voir L c. 7) b) p. 56-66. 

(10) Mais seulement approximativement : la (n — 1)-me lettre semble exercer sur la (n + 1)-me 
une influence directe non négligeable bien que moins forte que celle qu’elle exerce indirectement 
par l'intermédiaire de la #-me ; et il y a encore la (n — 2)-me lettre, etc. Voir les résultats 
expérimentaux de MARKOFF !, c. 

(11) Cette idée a été émise, sous une forme plus précise par Serge BERNSTEIN, Mathematische 
Annalen, t. 97 (1926), p. 1-59, voir p. 40-41. 
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V 


Je vais calculer quelques espérances mathématiques, utiles à l'étude 
de la structure de la « contagion ». Ce calcul sera des plus simples 
si nous nous appuyons sur le fait suivant : 

La structure de la «contagion » est homogène, c’est-à-dire que toutes les 
variables aléatoires X1, X2, .…. Xn, … ont les mêmes probabilités à priori. 
Outre cela, l’interdépendance entre deux variables quelconques xx et 
x est la même. 

Cela est plausible puisque la probabilité (x) est une fonction symé- 
trique des variables qui précèdent xr71, mais cela ne va pas de soi. 
Le résultat énoncé contient une propriété importante de la structure 
de la contagion qui sera démontrée sous une forme un peu plus géné- 
rale au n° VII et dont je ferai usage dès maintenant. 

Valeur probable de Xn. l'espérance mathématique de x»7 pour # 
quelconque est la même que celle de x,, à cause de l’homogénéité. 
x, ñne prend que deux valeurs, z et 0, avec la probabilité p et s res- 
pectivement ; on a donc () 


(7) Ex, — Exp" 00 ip: 
Ecart de xn. Le carré de l'écart moyen quadratique de x: est, lui 
aussi, indépendant de #. On a 
(8) Era — p} = Eu pŸ = pr — p}° + (0 — pb} = ps. 
Coefficient de corrélation de xx et x1. I1 suffit de calculer, d’après 
le principe avancé, le coefficient de corrélation de x, et de x,, celuide 
xx et x: étant le même. Commençons par calculer la probabilité d’avoir 
les deux équations simultanées 


X1 — I, Xe — I, 


c’est-à-dire que les deux premiers tirages donnent rouge. La proba- 
bilité d’avoir x, — 1 est p. Etant donné x, = 1, la probabilité d’avoir 


Ô MR £ Por 
x, — 1est, d’après (1), +. La probabilité cherchée est la probabilité 
composée 
Pb 
D * 
+0 


(12) E désigne l’espérance mathématique. 
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En calculant de la même manière encore trois probabilités analogues, 
on obtient la table de corrélation de x, et de x; 


Xi | 
Pa JE O 


z | PP) LM 
F0 I +OÔ 

où (je) (so + à) 
s De EN) I +Ô 


Observons en passant que cette table est symétrique par rapport à 
sa diagonale principale, donc x, dépend de la même manière de x, que 
x, de x,. En l’utilisant on trouve 


OR Ne et ee tee 


Le coefficient de corrélation entre xx et x est, en vertu de (8) et 


de (9), 
1 no Me M 
(xo) VElxx DE p)? Ex, — p} I+OÔ’ 


il a le signe de et il croît avec d, s’approchant de l'unité quand d 
est positif et très grand, c’est-à-dire lorsque la contagion est forte. 

Valeur probable de r — x, + x, +... x». —TJI,a variable aléatoire r 
est le nombre, sujet au hasard, des boules rouges sorties aux # premiers 
tirages. Sa valeur probable est 


(1x1) Er = Ex, + Ex: +... Ex, — np 
en vertu de (7). 
Ecart de r. — On obtient son carré, en utilisant (8) et (0) 
(22) Efr— np) = Efxi — p + x — p + ee + x, — p} 
= REX, — pŸ + nn — 1) Elxi — p) (x: — p) 
_ (n — 1) 
Coeficient de dispersion. — On obtient le carré de ce coefficient 
(qui sera comme d'usage, désigné par Q) en divisant E(r — np)? 
pat #p5, ce qui serait la valeur de (12) si x,, x2, … Xn étaient indé- 
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pendants. Donc 
RE ni 
(13) re ans 


On a Q > rou Q<1 selon qued > 0 oud <o. C'est-à-dire que 
la dispersion est hypernormale ou hyponormale selon que les chances 
augmentent où diminuent avec le succès. Cela se comprend. 


VI 


La structure de 1’ « hérédité », telle qu’elle a été définie au n° IV, 
n'est pas nécessairement homogène. Désignons par px la probabilité 
à priori (voir n° IT) d’avoir x: — 1 et par 5, celle d’avoir xx — 0. On a 
(x4) Pe + % = I. 


Cherchons la relation entre px et px+1. 

Que la boule tirée au (k + 1)-me tirage soit rouge peut arriver de 
deux manières : elle est sortie de U, ou de U, selon que la boule précé- 
dente a été rouge ou noire. I,;a probabilité totale est donc 


ou, utilisant (6) et (14) 


(x5) Pts = TOx Poe 
Déterminons p par 

(16) P = TP + Po 
On aura 


here re) 


et, en répétant le raisonnement, 


(7) Prtm — P = (px — p). 
En posant À — 1 dans (17) on voit qu’en supposant 
(x8) MAP 1 H=ITP 
on obtient p-— ps = ps —="""— p c'est-à-dire que la série des 


variables aléatoires x,, X2, X:, … est homogène. 
A partir d'ici, j'introduis la nouvelle hypothèse (18) comme partie 
essentielle de la définition de « l’hérédité » dont la structure devient 
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par cela homogène. Observons encore, que d’après la définition com- 
plétée, cette structure ne dépend que de deux paramètres, p et +. 
En effet en posant 


(19) CES 
on tire de (4), (6), (16) et (x9) 

(ee de Er 
(20) lo =? — GS — 5 + TO 


En vertu de (18), (19), (20), la composition des trois urnes U,,U,, U:, 
est exprimée par p etr. On peut choisir dans l'intervalle 0 << 7. 
Etant donné p, r doit être tel que 2e, 95, 5, 5, données par (20), soient 
comprises entre Oo et I. Si l’on suppose » £ 5 il est nécessaire (et aussi 
suffisant) que 


(21) Re 0 


Calculons encore la probabilité d’avoir Les deux équations simul- 

tanées 
Xx —= 1, Xyim = I 
c’est-à-dire que le k-me et le (k + m)-me tirage donnent rouge les deux. 
La probabilité d’avoir x4 — 1 est p. Etant donné xx — 1, on obtient 
la probabilité d’avoir xxym — 1 comme le pr+m de la formule (17) 
à condition d’y mettre px — 1. La probabilité cherchée est la proba- 
bilité composée 
pl + Ts). 

En calculant de la même manière encore trois probabilités ana- 

logues, on obtient la table de corrélation de xx et de Xkym : 


X} 
JR O 
Xxm 
T e(e ee To) s(o ET to) 
O p(o — 7") o(o «+ rio) 


Observons que cette table étant symétrique par rapport à sa diago- 
nale principale, la variable xx est liée de la même manière à Xezm 
que Xk4m l'est à xx. Observons encore que les probabilités qui y sont 
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contenues ne dépendent que de "», valeur absolue de la différence 
des indices À et À + 7". 

En résumé : 

La structure de l’« hérédité » est homogène, c'est-à-dire que toutes les 
variables aléatoires X;, X,,:::Xn,-.- ontles mêmes probabilités à priori. 
L'interdépendance entre deux variables quelconques xx et x1ne dépend 
que de leur distance |k — I]. 

Nous pouvons maintenant calculer tout comme au n° V, les espé- 
rances mathématiques les plus essentielles à l’étude de la structure 
de l’hérédité. 


Valeur probable de xn. À cause de l’homogénéité 
(22) EXxn = EX = p°I + 0.0 — p, 
Ecart de xn. Ie carré de cet écart est 
(23) Ex — p} = Ex, — p} — po. 


Coefficient de corrélation de xx el Xkim. La table de corrélation 
trouvée ci-dessus donne 


(24) E(xx An e) (Xp = e) _ par”, 
On obtient de (23) et (24) le coefficient de corrélation 


E(xx - = o) (Xx4m Mes p)? 0 
(25) VEG (X% = p}? Hixre - p)° k 


I1 est positif dans le cas de la ressemblance (r > 0); dans le cas de 
la dissemblance, il est négatif pour m — I et son signe alterne quand 
parcourt successivement les nombres 1, 2, 3,--.+ ; dans les deux cas sa 
valeur absolue diminue si |r| diminue ou si # augmente. Tout cela se 
comprend, si l’on envisage l’exemple des voyelles et des consonnes 
ou celui des descendants successifs, fils, petit-fils, arrière petit-fils, etc. 

Valeur probable der — x, + x, +:::+ xn. — C'est le nombre moyen 
des boules rouges sorties en n tirages. C’est Er — np. 

Ecart de r. — On obtient pour son carré, en utilisant (23) et (24) 


(26) ER he) Ex: — pr x, —p +. x —p} 
— NE(x — p)? +2(n — 71) E(xi — 0) (Xe —p) +: + 2E(x1 — 0) (xr —p) 
= ps + 2(N — I)oot + 2(n — 2)p07? +. + 2portT 1. 
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Coeficient de dispersion. On a 


Hein 2 LUE MNT role RE ent 
CR ns ae 1} + (n—2)? + (n—3}° + ] 
TAC RÉ 
dns | ame du 


On voit, encore mieux par la première des deux expressions (27), 
que la dispersion est hypernormale (Q > 1) s’il y a « tradition » 
(7 > 0) et hyponormale (Q < x) s’il y a « opposition » (r < 0). 


VII 


Retournons à la structure de la « contagion » pour démontrer la 
proposition avancée au début du n° V. 
Calculons pour commencer la probabilité gr d’avoir 


PARA ES Ne at Xi = pe er, =Ù0 


c’est à dire la probabilité que parmi n tirages successifs les 7 premiers 
donnent rouge et les s derniers noir, où 


VAI SEEN: 


gr est une probabilité composée dont on trouve les 7 premiers facteurs 
en mettant # — O, I, 2,... 7 — 1 dans (1) et en donnant les justes 
valeurs à X;,, X,, . ; on trouve les s derniers facteurs en soustravant 
la valeur fournie par (r) de l'unité. La probabilité cherchée est donc 


NI 


re p(o +ô) (p+28)--. (0 (r—1x1)d 5 (s+ 8)... (sa +(s—x)à) 
PARTS I(1 +0) (1420). .-(1+(7—1)0) (x +78) (x + (7 +18): - (1 +(n—1)0)" 


La probabilité que parmi les n — 7 + s variables aléatoires 


X1, X», X3;, dan: Xn 


\ 


r déterminées soient égales à 1 et les autres s à 0, est également gr, 
quel que soit l’ordre prescrit pour les x et les 0, parce que la probabilité 
en question est composée de # facteurs dont les dénominateurs sont 
les mêmes que de ceux qui composent gr et les numérateurs ne sont 
changés qu’en ordre. 
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Introduisons » indéterminées (variables ordinaires) z,, 2, …. 2% et 
calculons l'espérance mathématique 
(30) Ez, "2,2%. TO 


Le terme 2,2 2 Sera multiplié par la probabilité d’avoir les 
équations simultanées (28), donc par gr. Mais un terme quelconque 
Zi Zi … 2, produit de 7 facteurs différents, sera également multiplié 
par gr. La somme de tous les termes analogues est la 7-me fonction 
symétrique élémentaire 


2122007 loin C1 2 rec er ue her EL Zn 
Donc 
(31) Ez*t2,*2. r “207 EN us qiZi cu G:La TR GnZn Te P(21, 22, 23, °° Zn) 


OÙ P(2, Ze, … 2») est un polynome symétrique de # variables. 
La connaissance de £(2, 22, .… 21) suffit pour répondre à une question 
de probabilité quelconque concernant X;, x,, …. Xn. Par exemple la 


probabilité d’avoir xx — 1 est donnée par 
Ces (EE 2: æ)) 
é | d2% MES, 


celle d’avoir xx — I et «1 — I par 


d (21, ÉD Zn) 
NAT 
E rot ( 22,92] BE — 0: 


etc. Mais (2. 22, … 2n) est symétrique en %, 22, … 2n. Par conséquent 
toutes les probabilités qui dépendent de X;, X2, .… Xn Sont parfaitement 
symétriques par rapport à ces n variables aléatoires. Cela contient la 
proposition à démontrer et davantage ; par exemple l’interdépendance 
entre trois variables quelconques xx, x: et xn est la même, quels que 
soient les indices, etc. 

Les calculs faits nous seront utiles encore pour d’autres questions. 
Par exemple la probabilité pr,s d'obtenir 7 fois rouge et s fois noir sur 
n = 7 + s tirages dans un ordre quelconque ou d’avoir, ce qui est la 


même chose, r — x, + X, +:.:+ xn — 7, est évidemment 


> ny = fras\plort)-:(p+ (112) o[at8)- {0 (sr) 
Ge) pe (je ee ee 5 
En mettant 4 = 2% = :::— 2% —e dans (31) on obtient la fonction 
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caractéristique de la somme r — x, + X2 +:: .L x», c'est-à-dire la 
fonction 


" n 
Y—0o —0O 


I 

{ 
Rd te 
TS. as 
r(S)T (5) o 

La dernière expression de qn(z) s'obtient à partir de (32) à l’aide 


de l'expression classique de l'intégrale eulérienne B par des transfor- 
mations faciles que je n’explique pas ici. 


4 (4) Tres + (I — x)]J'ax. 


VIII 


Nous devons encore connaître la fonction caractéristique de 
r—=X, + x, +...+ x, pour la structure de l’ «hérédité ». Pour gagner 
du temps, je cite ici sans démonstration la formule suivante de 
MaARKOFF (1) : Si pr, désigne la probabilité de tirer rouge 7 fois et 


noir s fois sur » — 7 + s tirages ou, ce qui est la même chose, d’avoir 
Pi Fond 


(we) 09 
0 DV pme A A Ro 


— (o + ot)u — (so + pr)0 + ru 
ÿ—O $s—0 


La fonction caractéristique cherchée est 


ñn 


(35) Eer: —= > Dineres — 22 pue dw, 


— 2m J  wnri 


Y—=0o 


la ligne d'intégration dans le plan de la variable complexe w étant 
un cercle parcouru dans le sens positif dont le centre est w — 0 et 
le rayon aussi petit que, z étant donné, la fonction analytique F(e:w, w) 
de w soit régulière à l’intérieur et sur le cercle. Cette forme n’est pas 


très explicite, mais tout à fait suffisante pour les passages à la limite 
dont on aura besoin. 


(z3) Voir L c. 7) b), p. 49, formule (10). On peut déduire de cette formule le cas limite traité par 
MaRKOFF ; pour une métlode simple, voir l'indication qui sera donnée à la fin du ne IX. 
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IX 


Rapprochons maintenant les deux structures d’interdépendance, 
celle de la « contagion » et celle de l” « hérédité » que nous avons traitées 
séparément jusqu'ici. 

Considérons pour les deux structures la loi de la variable aléatoire r, 
donc la probabilité pr,n-r d'obtenir exactement 7 fois rouge sur # 
tirages, 7 — 0, 1, 2,---"n. Cette probabilité est déterminée par la fonc- 
tion génératrice respective,(33) ou (35), même déterminée explicitement 
par (32) pour la première des deux structures. En portant les n + I 
probabilités 


(36) Don Din Da,n_2 APT RE Dao 


comme ordonnées équidistantes dans un ‘système rectangulaire, 
nous obtenons une distribution (une « courbe ») de probabilités ; elle 
dépend de 3 paramètres dans les deux cas, de p, à, n, respectivement 
de 0, 7, n. Comme en statistique nous n’attribuons de valeur qu'aux 
observations portant sur un grand nombre d'individus, il est naturel 
de chercher les formes limites résultant de cette distribution, quand #, 
le nombre des tirages, tend vers ©. Il y a différentes formes limites, 
la courbe de Gauss n'étant que la plus connue, puisque p et à comme 
p et Tr peuvent varier différemment avec #. 

On peut classifier les formes limites de deux points de vue différents : 
d’abord d’après l’ordre de grandeur de p, c’est-à-dire d’après la fré- 
quence de l'événement considéré, puis d’après l’ordre de grandeur de d 
respectivement de 7, c’est-à-dire d’après le degré d’interdépendance 
entre les événements individuels. 

Concernant la fréquence, je considère deux cas. 

a) Je parle d'événements wswels si p garde une valeur invariable 
(entre o et I, limites exclues) lorsque # tend vers l'infini. Dans ce 
cas no, la valeur probable des tirages rouges, tend vers ce. 

b) Je parle d'événements ares si la valeur probable np garde une 
valeur positive invariable quand # — ©, Dans ce cas p tend vers zéro. 

Je distinguerai trois degrés d’interdépendance. C’est Q, le coefficient 
de dispersion qui sera la base essentielle de la classification. 

1) Les événements indépendants sont caractérisés par d = 0, respec- 
tivement r — 0. Dans ce cas le coefficient de dispersion Q est égal à 
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l'unité et, pour les deux structures envisagées, la réciproque est vraie : 
Q = 1 implique d — 0, respectivement 7 — 0 en vertu de (7 3) et de (27). 

2) Je parle de contagion faible si d tend vers zéro, de manière 
que xd garde une valeur invariable (différente de o et supérieure 
à — r) et de faible ressemblance sir garde une valeur invariable 
(ç 0, — 1 <Tr<1) quand # —+ ©. Dans ces deux cas d'inéerdépen- 
dance faible Q tend vers une limite finie positive et différente de 1. 

3) Je parle de contagion forte si à garde une valeur positive inva- 
_riable et de forte ressemblance si 7 tend vers 1 de manière que (1 —:) 


garde une valeur positive invariable quand n —> œ. Dans ces deux cas 
2 
d'interdépendance forte *, tend vers une limite positive finie. 


Les trois degrés d’interdépendance ont été distingués d’après 
l’allure de Q, brièvement 


HOMO Or Q fini. 
vn 
Nous aurions pu aussi les distinguer d’après l’allure du coefficient 
de corrélation entre x, et x», les deux tirages les plus éloignés ; c’est 


S 


[0] . . . 
7 + respectivement #1, Nous avons dans les trois cas respective- 
— J 


ment 
Ô ù è . 
— 0 D) = © < fini 
1 Ts cr 355 
el — TTL, O -n— 1 fini 
(en admettant n(1 — r) — }, on a æ#-1 , e-à), 


La manière qui est la plus commode au point de vue du calcul, de 
caractériser les cas I), 2), 3) ainsi que a) et b) est celle adoptée dans 
leur définition. C’est d’assigner à un des paramètres 9, d où + ou à 
une combinaison simple d’un paramètre avec 7 une valeur constante 
quand #0. Je réunis ces constantes, caractéristiques des cas res- 


pectifs, avec la notation dont je me servirai plus tard, dans la tabelle 
suivante 


a) p b)no = h 
LUO'ES 2) nà = à 3) ) 
Gi (0 T NI — 7) = À, 
Observons à propos de a) et b) qu'il aurait été inutile d'introduire 
un cas auquel L->1 puisque l'on peut supposer 9 Æ 5. Observons à 
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propos du cas 3) que Q — « implique Q? > 1 donc d > 0, r 0. Obser- 

vons encore qu'on pourrait introduire un quatrième degré d’inter- 

dépendance, celui de l’interdépendance « ultraforte » caractérisé par 
4) 8 = + co 


Di OT 


: RME ù : 
ou par ce que le coefficient de corrélation Ts respectivement #7? 


est exactement 1. Ce cas est instructif comme terme de comparaison 
mais sans intérêt en lui-même : Le premier excepté, les coups sont 
forcés. Chaque tirage à partir du second suit nécessairement son 
prédécesseur, de manière que la probabilité d’avoir r — 0 est o, celle 
d’avoir r — n est p, tandis que celle d’avoir une valeur intermédiaire 
quelconque est o. 

D'après ce qui a été dit, il sera clair, ce qu’il faut entendre par le 
cas 14) par exemple. C’est le cas où les conditions 1) et 4) sont réunies, 
le cas d'événements usuels indépendants ( est fixé, à — 0, respec- 
tivement 7 — 0), le cas classique de BERNOULLI menant à la loi nor- 
male de GAUSS. 

Le cas 3b) c’est-à-dire le cas d'événements rares à interdépendance 
forte est celui où 


h, À, à sont des constantes positives. Dans ce cas on trouve facilement 
la distribution limite en envisageant po,r, la probabilité d’avoir r — 0 
c’est-à-dire de ne tirer que des boules noires pendant les 7 premiers 
tirages. Pour la structure de la « contagion » on a d’après (32) 


__ os + ô)-..(s + (mn —xr)à) 
ee LPO EU ER) 
h\ | h h 
=(1-;) en nn) 
h { I g I 
HO cit re mn lat 


pour celle de l’« hérédité » on trouve 
Ponte on (rip) (2 pfr 7) 
h h)\ n—1 
=(—;)(-4) 7 
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Donc la distribution limite est la même pour les deux structures et 
peu intéressante comme distribution de probabilité : r — © a la proba- 
bilité r et toutes les autres valeurs de r la probabilitéo. Donc le résultat 
est le même que si tous les tirages donnaient forcément noir dès 
le commencement. 

Pour les autres cas la distribution limite est plus intéressante mais 
moins facile à calculer. Prenons par exemple le cas 3a), le cas d'évé- 
nements usuels à contagion forte, si nous nous bornons à la pre- 
mière des deux structures. Dans ce cas 2 et à sont des constantes, 


0 << p 1,9 > 0. Introduisons la nouvelle variable aléatoire x — s 


qui prend la valeur s avec la même probabilité que r la valeur r. Sa 


fonction caractéristique est donc, d'après (33) 


#14 


(I — x) ic fe + (Ii à ax. 


| © H 
—.. 
j 
oo 
Lu 
©? © 


En observant que 


3 ñ « 
PE Poe = Ë Rte en He er 
> n 20° 


pour #7 —® on trouve la limite de la fonction caractéristique. On en 
déduit (#) que la distribution des probabilités de la variable aléatoire 
= =. devient continue à la limite, la densité de probabilité au point 


X — x étant 

a 

FN Su 
(37) ————— %x (1 =)" 


C'est de cette manière, en partant des fonctions génératrices (33 
et (35) qu'on peut calculer les distributions limites pour les deux 


(14) Dans ce cas comme dans tous ceux qui figurent dans ces leçons, le théorème que j'ai donné 
sur la limite de la fonction caractéristique (Mathematische Zeitschrift, t. VIII (xo20), p. r7r-18r) 
et le théorème analogue que M. Paul LÉVY a donné un peu plus tard (Comptes-rendus, t, 175 
(1922) p. 854-856) sont également applicables, 
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structures dans tous les cas qu’on va énumérer au n° suivant. L/exemple 
que j'ai traité est un des plus simples ; en particulier le calcul de (35) 
exige un choix judicieux de la ligne d'intégration. Mais je n’ai pas le 
temps de m'’arrêter aux détails de ces calculs. 


X 


Les distributions de probabilités que nous allons considérer, résultent 
toutes de la loi de distribution de la variable aléatoire r. On se rappelle 
que r prend les valeurs 


0, Fe 2, h 
avec les probabilités respectives 


Done Pins Pa,n-2s FT Pn,0 


et que r n’est autre chose que le nombre des boules rouges sorties aux 
n premiers tirages, le schéma d’urne pouvant être celui de la « conta- 
gion » ou celui de 1’ « hérédité ». Quelquefois nous devrons introduire 
une nouvelle variable aléatoire x — ar + b, transformée linéaire der 
avec des coefficients 4, b qui dépendent de # ; x sera à la limite une 
variable aléatoire continue, c’est-à-dire qu’elle prendra toutes les valeurs 
d’un certain intervalle ; mais la loi de probabilité de x peut être 
discontinue., Quelquefois nous appellerons r la variable de la loi limite. 
Mais alors r pourra prendre des valeurs entières non-négatives quel- 
conques 


L. 2, «+. M, 
avec des probabilités qui seront désignées par 

Po Pis Pa; CT Pins 
de manière qu’on ait 


Dr —= lim Dr eme 
N—>0 


Nous allons considérer 9 différents passages à la limite. Ce sont les 
cas 14) 1b) 2a) 2b) 3a) expliqués au n° précédent, le cas f) ou cas de 
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transition qui sera expliqué plus loin et les cas 3b) 4a) et 4b), triviaux 
en eux-mêmes mais instructifs pour la comparaison (1). 

1a) Evénements usuels indépendants. p constante, 0 p<1,0=7T—0, 
les deux structures envisagées coïncident dans le cas particulier com- 
mun de l'indépendance. En posant 


rare 


8 D 
(8) as 


n2s 
x prendra à la limite toutes les valeurs de l'intervalle — co, + ©. 
La densité de la probabilité au point x — x est 


C’est la courbe normale ou de Gauss. Ce résultat classique, la « Loi 
des grands nombres » est dû à De MOIvRE. 

10) Evénements rares indépendants. La valeur probable ns =} 
est constante, d — r — o, c’est encore le cas d'indépendance, cas 
particulier commun aux deux structures. La probabilité Pr d’avoir 
r — m à la limite est donnée comme coefficient de la fonction généra- 
trice 


(39) Ezr — Do + Dir +: + pra +... 


eh eh hn 
! 


= eh + at DA SR ehtz—1), 


Ce résultat est encore classique, c’est la « Loi des petits nombres » 
due à Porsson. 

24) Evénements usuels à interdépendance faible. o constante, 
O<p< 1, nd — d constante, r constante, p + d = 0 à cause de (3) 
et les quatre expressions (20) sont positives. En introduisant x par (38), 
x varie à la limite de —-œ à + ©, La densité de probabilité au point 
X — x est 


(25) Les cas 1a) 1b) sont classiques, ainsi qu'un cas particulier de 2a) (celui de la boule non remise). 
Pour les autres cas limites non triviaux voir : a) Structure de la « contagion »: tous les cas 
chez MARKOrr L. 6, 8) et EGGENBERGER et P6LYA L. c. 8) a) et b). Les courbes limites continues 
34) et {) furent rencontrées déjà avant par PEARSON ; voir n° XI. 8) Structure de l’ « hérédité » : 
24) MARKOFF L 6. 7) a) b) ; 2b) A. AgrpLt, Thèse, Zurich 1924 ; 3a) et 1) ont été indiqués la pre- 
mière fois, à ma connaissance, dans ces leçons. 


» 
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selon qu'il s’agit de « contagion » ou d’« hérédité ». C’est la courbe de 
GAUSS, comme au cas Ia), mais la dispersion diffère de la normale, 
le signe de la différence étant celui de d ou de 7, selon la structure. 
2d) Evénements rares à interdépendance faible. np = h, m=detr 
sont des constantes positives, 7 1. (d < 0 est exclu par (3), r <0 
par (21).) La probabilité px d’avoir r — m à la limite est donnée par 


(40) EP ho + hais + cs + Dot He = (1—(2— 1)d). d 


s’il s'agit de « contagion » et par 


(41) Er = Po + Dir +. + pur He —e 
s’il s’agit d’« hérédité ». 

3a) Evénements usuels à interdépendance forte. p, d, n(1 — 7) = À 
sont des constantes positives, o °° p <[ 1. On introduit la variable 


(42) Je 


qui prendra à la limite toutes les valeurs de l'intervalle o = x < 1. 
Pour la structure de la contagion on a une distribution limite conti- 
nue ; nous avons calculé auparavant sa densité au point x — x, voir 


G7): 


Pour la structure de l’« hérédité » on a une distribution limite en 
partie continue, en partie discontinue. Les valeurs extrêmes 
Me 0 TEST 
ont des probabilités finies, respectivement 
pÀ —5À 


GER pe ; 


les valeurs intermédiaires n’ont que des probabilités infiniment petites, 
la densité de la probabilité étant une fonction continue possédant 
toutes les dérivées. On déduit ces propriétés d’une discussion appro- 
priée de la fonction caractéristique 
a es à 
1 sinh É | 


= (3—À) I 
po 2" À cn a) LA UE Papa A) Ses D ver 


A = V2? + 2(2p — 1)2 + 
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et, comme d'usage, 


2 cosh x = e + er, 2 Sinh x — et —"e7*. 
Je note encore ici le carré de l'écart de x (1). 


E(x — p}° = € (e 


3b) Evénements rares à interdépendance forte. np, à, n(1 — 7) = À 
sont des constantes positives. Le cas a été traité à la fin du n° IX. 
La distribution limite peut être conçue de deux manières. On peut 
prendre r comme variable ; alors à la limite la probabilité d’avoir 
r = 0est 1 et celle d’avoir r — m est o si m>0.Ou bien en posant 
x — ar, X variera d’une manière continue dans un intervalle fini ou 
infini, dont l'extrémité gauche est le point x — 0 ; ce point aura la 
probabilité 1, tous les autres points ensemble la probabilité o. 

4a) Evénements usuels à interdépendance ultraforte. p constante, 
Op L1,9 = + et r — 71, les coups sont forcés à partir du second ; 
les deux structures coïncident tout à fait dans ce cas particulier. 
X étant défini par (42), la probabilité d’avoir x — o est 9, celle d’avoir 
Z.— rest pet celle d'avoir 0 x Testo. 

4b) Evénements rares à interdépendance ultraforte. p = 0, d = +, 
r — 1, les coups sont forcés à partir du premier. Comme distribution 
limite c'est la même chose que 3b). 

l) Cas de transition. Dans ce cas 


sont des constantes positives ; @ est une fraction constante, O < O< : ie 


Ce cas est, aussi bien au point de vue dela fréquence qu’au point de vue 
de l’interdépendance, la transition entre les cas 34) et 2b) comme 
on le voit du tableau suivant qui montre les puissances de # auxquelles 
np et Q? sont proportionnels. 


ne : an tn: —® b}n° 

Q:  3)n tn  2)mt. 
(16) Pour une connexion possible de cette formule avec des questions de physique voir G. J. TAy- 
LOR, Diffusion by continuous movements, Proceedings of the London M. $. 2e série, t. XX (1921), 


P. 196-210, en particulier p. 198-2017. 
(17) On peut arriver à la même loi limite sous des conditions un peu moins restrictives. 


Et 
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Posons 


Pour les deux structures, x variera à la limite continuement de 
O à ©, 

Pour la structure de la:contagion, on a une distribution limite 
continue ; sa densité au point x est 


XD ee. 


Euh 

(a) 

Pour la structure de l’hérédité, on a une distribution limite en 

partie continue, en partie discontinue. La probabilité d’avoir x — o 

est finie, — e-4. Les valeurs positives de X n’ont qu’une probabilité 
infiniment petite, la densité de probabilité au point x étant 


er Va h(2iVax) —0—X. S Ce 2e 
i/x AS (n+1)!n! 
X 5 ; Ë 


où J,(x) désigne la fonction de BESSEL d'indice 1. 
Je termine cette énumération par une tabelle et un diagramme que 
voici. 


| 10 — 1h 
Cas Limites F2 63 | 
I4 — ©, 00 (e) 2 
24 — ©, © I 2 = _ 
2b — 2 — y'a IN 
34 OL 2 4 3a ————— 3h 
(30) Ca CRUE 
(44) (SE I 3 
(4b) + CEA AS | 
t O, co I 3 a 4 


Figure 1. 


La première colonne de la table contient le nom du cas limite. 
I1 faut se rappeler que 1, 2, 3, 4 y désignent des degrés croissants 
d’interdépendance, a les événements usuels, d les rares. Les cas tri- 


= 55 — 
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viaux sont entre parenthèses. La seconde colonne donne des limites 
de la variation de la variable aléatoire, si cette variation est continue ; 
un trait horizontal désigne les cas où cette variable ne prend que les 
valeurs 0, 1,2, 3, .… ; les cas 3b) et 4b) sont traités ici comme discontinus. 
La 3-me colonne donne p, le nombre des paramètres de la distribution 
sous sa forme considérée dans l’énumération, et la 4-me colonne, pour les 
distributions à variable continue, le nombre total © des paramètres 
qu'on peut y introduire en faisant usage d’une transformation linéaire 
arbitraire. 

Dans le diagramme toutes les distributions, aussi les discontinues, 
sont considérées comme admettant une transformation linéaire arbi- 
traire de la variable aléatoire. Alors la différence entre 1a) et 2a) 
disparaît ; c’est pour cela qu'ils sont réunis par un signe — (mis 
verticalement) de même que 4b) et 3b). La flèche allant de 2b) à rb) 
indique que 1b) est cas limite de 2b). En effet, en faisant 4-0 dans 
(40) ou en mettant r — o dans (41) on obtient les deux fois (39). Inver- 
sement, la distribution 2b) qui contient 2 paramètres ne peut pas être 
cas limite de 1b) qui n’en contient qu’un. Les autres flèches ont une 
signification analogue. Le cas 3b) est considéré ici comme continu. 
Il est à remarquer que la courbe de GAUSS, ra) ou 24), est cas limite 
de 1b) 2b) t) et 3a), donc de toutes les autres distributions non tri- 
viales. Qu'on observe aussi la situation de cas t), intermédiaire entre 
3a) et 2b) mais, d’une autre manière, aussi entre 3a) et 2a). 

C’est l’analogie des deux structures de la « contagion » et de « l’héré- 
dité » qui me paraît le plus remarquable en tout ceci. Elle est tellement 
complète qu'en construisant la table et le diagramme, il était super- 
fÎlu de séparer les structures. 


XI 


Les résultats énumérés au n° précédent ouvrent la voie à d’autres 
recherches. Je dois me borner à quelques indications. 

1. L/énumération des cas limites est-elle complète ? Je ne saurais 
donner une réponse définitive (15). J'esquisserai une réponse partielle 


(18) En tous cas il manque le cas où la variable aléatoire X étant déterminée par (42), X = po 
a la probabilité 1 à la limite. Ce cas limite, trivial en soi-même, peut être atteint d’une manière 
intéressante, par exemple lorsqu'il y a une «opposition extrême » dans la structure de l’ «hérédité » : 
I 1 (— 1} 
0=G—=- T = — 2 —  —- 

24 n 110 pe 


2? Q 
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en me servant d'une modification de la méthode par laquelle M. Karl 
PEARSON a obtenu primitivement les distributions qui proviennent 
ici, aux cas 3a) et {), de la structure de la « contagion ». Pour cette 
structure, on obtient de (32), en écrivant pour abréger p- à la place : 
de Prin-r, 


Prin — Pr __ (x —28)(r + 1) — (n +r)(p — 6) 

“5 Pres (78 + p)(n — r) ï 
\ $ Pr — Pr, __ (1 —28)r — (n + 1)(o — À) 
(44) PRE ER TR TE) \ 


Si nous considérons 7 comme variable, le dénominateur de la frac- 
tion dans le second membre est un polynome de second degré dont les 
racines sont réelles : 


r=—f, n pour (43), 
ç 
V—= 70, Moser pour (44). 


Le domaine de variabilité de la variable aléatoire 7 est compris entre 
0, première racine pour (44) et #, seconde racine pour (43). 

Une transformation linéaire x — 4r + b changera (43) et (44) en 
équations analogues avec des propriétés correspondantes des racines. 
Admettons que dans un cas limite où la variable aléatoire est x, il 
y ait une densité de probabilité y — œ(x) qui est une fonction déri- 
vable de x ; admettons en outre que la convergence vers le cas limite 
soit « suffisamment bonne ». Alors (43) et (44) coïncideront à la limite 


dans une même équation de la forme 


EAU do + X 
(45) y dx bd, + bix + box? 


où Le dénominateur bd, + b,x + b,x? n'a pas de racines 1maginaires. 
Observons que le degré du dénominateur peut s’abaisser à la limite 
si l’une ou l’autre des racines devient infinie. On voit que le domaine 
de variation de x est compris entre les deux racines, entre © (ou — co) 
et l'unique racine ou enfin entre — © el + © selon que le degré du 
dénominateur est 2, I, OU O. 

Faisant abstraction d’une transformation linéaire, on peut admettre 
que le dénominateur du second membre de (45) est x(1 — x) où x our. 


ANNALES DE L'INSTITUT H. POINCARÉ. IL 
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Ces trois cas figurent effectivement dans l’énumération du n° X. 
On obtient en partant de (43) ou de (44) 


If) + (3x) £ 
“dr to 5 AE au Cas 34) 


ù # 
Dennee au cas f) 
14 x 

D Ed au Cas 24). 


Les racines de b,, + b,x + b,x? étant supposées réelles, une discussion 
facile des intégrales y de l'équation (45) montre que, abstraction 
faite toujours d’une transformation linéaire, ces trois cas 34) ) 2a) 
sont les seuls où /'ydx converge entre les limites admissibles (racines 
et Ho comme spécifié ci-dessus). 

D'après cela, en tant qu’il s’agit des distributions limites conéinues 
de la structure de la « contagion », notre énumération est complète. 

2. Comme il a été dit, les considérations précédentes ne sont qu’une 
modification de la première déduction que M. PEARSON a donnée de 
ses courbes de fréquence bien connues provenant de (45) (1). 

Je n'ai pas à discuter ici, bien entendu, l’intérêt des applications 
que M. PEARSON et son école ont faites de ces courbes. Tout ce que 
je dirai se rapportera à la première introduction citée, entre laquelle 
et le présent traitement, il convient de signaler les différences suivantes: 


(a) Je considère ici (43) et (44) à la place du rapport AP — PI a :. £ 


considéré par M. PEARSON. 

(b) Ie point de départ est ici le schéma d’urnes de la « contagion », 
schéma plus général que celui de la boule non-remise dont s’est servi 
M. PEARSON. 

(c) Les cas considérés ici sont des cas limites ; pour chacun d’eux 
j'ai indiqué soigneusement de quelle manière il faut varier les para- 
mètres du schéma initial pour y arriver et quelle est la signification 
du passage à la limite pour la structure des probabilités. Les cas 
considérés par M. PEARSON proviennent de l'équation (45), bien 
sûr, mais ils ne proviennent pas pour cela d’une structure de proba- 
bilités, puisque M. PEARSON fait varier les paramètres d’une façon 


(19) Philosophical Transactions, t. 186 (1895), p. 343-414. 
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quelconque, sans se soucier, comme d’ailleurs il le dit explicitement lui- 
même, si ses formules restent ou ne restent pas susceptibles d’inter- 
prétation dans les termes du problème initial. 

La dernière remarque explique pourquoi nous n’avons pas obtenu 
ici toutes les courbes de M. PEARSON bien que nous soyons partis d’un 
schéma d’urnes plus général (2). 

3. Notre déduction donne à chacune des courbes limites obtenues 
une signification, la rattache à une structure d’interdépendance qu’on 
pourra, avec une approximation plus où moins grande, reconnaître 
au matériel statistique fourni par l'observation. 

Un exemple suffira pour montrer le genre d'application que j'ai 
en vue. Envisageons la statistique d’une épidémie rare, comme la 
petite vérole ou le choléra le sont de nos jours. Les individus observés 
sont les habitants d'un pays. On enregistre les cas mortels par suite 
de l’épidémie par mois. A la fin d’une certaine période on aura constaté 
M, mois sans cas mortels, M, avec exactement 1, … M; avec exacte- 
ment 7 cas mortels. Il s’agit de donner une expression théorique à 
laquelle la suite 


(46) PR OS ROSE ae 


pourra être assimilée. 
Observons d’abord que l'épidémie est are. Donc il est naturel 
d'examiner si la suite (46) est porportionnelle ou non à la suite des 


probabilités 
(47) Por Pas Pa Pre 


données par une des formules (30) (40) (41) relatives aux événements 
rares (cas Ib) et 2b)). 

Observons ensuite que les cas d’épidémie ne semblent pas être 
indépendants : un mois avec beaucoup de cas est un mauvais présage 
pour le suivant. Cela élimine à peu près 4 priori la formule (39) du cas 
1b) d'événements rares indépendants. 

Reste à décider entre les deux formules (40) et (41) du cas 2b) ; 
la première provient de la structure de la « contagion », la seconde de 
celle de 1 ’« hérédité ». Dans la première structure l'influence mutuelle 
mesurée par le coefficient de corrélation est la même entre deux indi- 


(20) Que les autres courbes de M. PEARSON ne sont pas cas limites du schéma d’urnes de la 
contagion, a été montré encore d’une autre manière, plus rigoureuse par M. EGGENBERGER. 
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vidus quelconques de la population, dans la seconde l'influence 
mutuelle décroît rapidement avec la distance croissante des individus 
(supposés rangés en s'rie linéaire). On dirait que la réalité est quelque 
part entre les deux structures. Si l'on envisage les schéma d’urnes, 
c’est la première structure qui semble être plus près de la réalité 
(voir le commencement du n° III). 

Les calculs numériques exécutés (2) montrent que (39) (indépen- 
dance) est décidément en désaccord avec les faits; (41) («hérédité ») 
donne une approximation meilleure et (40) (« contagion ») va encore 
mieux. S'il s’agit de la statistique de cas mortels par suite de certains 
accidents industriels, l'indépendance est encore en désaccord, et les 
deux structures de la « contagion » et de l’« hérédité » fournissent 
des approximations à peu près égales. 

4. Il y a des explications plausibles, il n'y en a pas d’obligatoires, 


au moins pas en statistique. Un exemple frappant de ce principe 
I 


est donné par la structure de la «contagion » au cas 3a) ; si p = & —d— 
la densité de probabilité correspondante, donnée par (37) devient 
égale à 1 pour 0 << x < 1, c’est-à-dire que la distribution de probabilités 
devient uniforme dans l'intervalle 0,1. Personne n’a trouvé nécessaire 
d'expliquer une répartition uniforme de probabilité par quelque chose 
d'autre, mais nous voyons qu’elle peut être expliquée par une structure 
compliquée d’interdépendance. Sans doute elle pourrait être expliquée 
encore d’une infinité d’autres manières (??). 

Voici un autre exemple. Si on rencontrait des séries d'événements 
alternatifs et homogènes montrant une dispersion anormale, mais 
réparties d’après la loi de GAUSS, la statistique correspondrait à notre 
cas 24). On pourrait l'expliquer tout aussi bien en supposant une 
Structure de « contagion » qu’une structure d’«hérédité », bien que ces 
structures soient différentes; mais elles mènent, dans le cas 24), à la 
même répartition des fréquences. 

5. Les deux structures envisagées, celle de la « contagion » et celle 
de l’« hérédité » sont différentes, mais leurs cas limites montrent un 
parallélisme remarquable, comme on à vu aux n° IX et X. Dans 


(21) Pour les calculs avec (39) et (40) voir EGGENBERGER L. c. 8). Les calculs concernant (41) 
ne sont pas encore publiés. 
) (22) Des remarques analogues s'appliquent au cas 5 = p = 5, n fini, plus intéressant encore 
à certains égards. 
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quels cas ce parallélisme va-t-il, malgré la différence des structures, 
jusqu’à la coïncidence des lois limites ? 

I faut écarter ici le cas de l'indépendance ainsi que celui de l’inter- 
dépendance ultraforte, donc les cas 


Id) 1b) 44) 4b) 
parce que dans ces cas particuliers communs les deux structures ne 
sont pas effectivement différentes. Elles sont effectivement différentes 
mais donnent la même loi limite dans les cas 


24) 3b) 
_et donnent des lois limites différentes dans les cas 


34) l) 2b) 

Laissons encore de côté le cas 30), trivial et pas trop différent de 4b). 
Alors le cas 2u) reste le seul cas où la même loi limite a pu s'imposer 
aux deux structures, malgré leur différence effective, et cette loi limite 
commune est la loi de GAUSS. 

Ce fait me semble symptomatique ; il nous fait entrevoir que la loi 
de GAUSS jouit d’une situation exceptionnelle non seulement dans la 
théorie des probabilités indépendantes, mais aussi dans celle des 
interdépendantes. En effet, en poursuivant les recherches de MARKOFF, 
M. Serge BERNSTEINN a réussi à montrer dans un travail trèsimportant(?) 
que c’est la loi de Gauss qui doit résulter de la superposition d’un 
grand nombre d'événements interdépendants, pourvu que l’interdé- 
pendance soit suffisamment « faible ». (Je n'ai énoncé qu’une des 
conditions et d’une manière assez vague.) 

Du résultat de M. BERNSTEIN on ne peut pas conclure qu’en dehors 
de la loi de Gauss il n’y a pas d’autres lois qui puissent s'imposer 
en même temps à des structures d’interdépendance très différentes. 
On ne peut pas conclure cela des exemples traités dans ce travail 
non plus, mais c’est ce qui est rendu plausible par tout ce qui a été 
dit. I1 serait important de préciser ces remarques et de trouver quelque 
propriété de la loi de Gauss qui caractérise sa place dans la théorie 
des probabilités interdépendantes comme elle a été caractérisée dans 
celle des indépendantes. Ce problème est encore bien obscur ; toutefois 
j'ai tenu à le signaler avant de terminer ces leçons. 


C3) Wen) 
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PAR 


M. Pauz LÉVY 


MESDAMES, MESSIEURS, 


Permettez-moi d’abord d'adresser à M. Emile BOREL l'expression 
de ma reconnaissance pour le grand honneur qu’il m'a fait en m'in- 
vitant à prendre la parole dans cette salle. Permettez-moi aussi de 
prévenir ceux d’entre vous qui ont entendu les remarquables confé- 
rences de M. POr.yA, que la présente conférence a un objet bien diffé- 
rent, et plus modeste ; ce ne sont pas des résultats nouveaux, maïs des 
résultats connus depuis longtemps, et pour la plupart classiques, 
que je vais avoir l’honneur d'exposer devant vous. 

1. CONSIDÉRATIONS GÉNÉRALES SUR LES ERREURS DANS LES MESURES 
PHYSIQUES. —- Le but que je me suis proposé est de vous faire com- 
prendre le théorème fondamental de la théorie des erreurs. D’après 
ce théorème, l'erreur accidentelle, commise dans la mesure d’une gran- 
duer physique, obéit à la lor de GAUSS. 

I1 faut d’abord préciser ce qu’on appelle erreur accidentelle. Le 
caractère essentiel d’une telle erreur est de dépendre du hasard ; 
si l’on recommence la mesure dans les mêmes conditions, l’erreur 
accidentelle variera ; elle sera tantôt positive, tantôt négative, et 
l’on peut espérer que ces différentes erreurs se compenseront si l’on 
prend la moyenne des résultats obtenus ; la précision est alors aug- 
mentée. L'erreur qui n’est pas accidentelle est dite systématique ; 
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on ne peut pas espérer la corriger de la même manière, et il n'y a pas 
de bonne mesure si l’on n'arrive pas à éviter ou corriger toutes les 
causes d'erreurs systématiques. 

La distinction entre ces deux sortes d'erreurs est souvent assez 
délicate. Contentons-nous de le montrer par un exemple, celui de 
l'erreur de lecture. Supposons pour fixer les idées que l’on ait à lire 
la position d’un trait de repère en face d’une règle graduée en milli- 
mètres, et qu’on se contente de lire le nombre entier de millimètres, 
sans chercher à apprécier la partie fractionnaire £ dunombre inconnu X. 
Si cette partie fractionnaire est inférieure à o, 4, pour fixer les idées, 
on lira le nombre entier # immédiatement inférieur à x; deo,6à1 
au contraire, on lira # + 1 ; entre 0, 4 et 0, 6 on hésitera ; on lira 
tantôt n, tantôt n +1, la probabilité de cette seconde valeur augmen- 


z-X 


tant avec £, suivant une formule que nous supposerons linéaire. Dans 
ces conditions, cela ne sert à rien de recommencer la lecture ; un 
même observateur ne pourrait guère oublier le résultat de la lecture 
précédente et relirait le même nombre, # où # + tr. Si l’on prend la 
moyenne des lectures faites par un grand nombre d’observateurs 
indépendants, on obtient une valeur x, fonction continue de X = # + £, 
et l'erreur commise x — X est une fonction continue de Ë, représentée 
par la figure ci-contre. C’est une erreur systématique ; la répétition 
des expériences ne permet pas de la corriger. 

Supposons au contraire que, pour mesurer la distance de deux traits 
de repère, on déplace à chaque lecture la régle graduée : £ varie d’une 
lecture à l’autre, les différentes valeurs entre o et 1 étant également 
probables. La moyenne des différentes erreurs x — X diffèrera alors 


L 
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peu de sa valeur théorique, qui est la valeur moyenne de la fonction 
x —- X de £ que nous venons de définir quand £ varie de o à I ; cette 
valeur est zéro ; il en sera ainsi à chaque extrémité de la règle. Dans 
ces conditions on peut arriver à une précision théoriquement indéfinie : 
l'erreur de lecture est devenue accidentelle ; on peut maintenant la 
corriger en recommençant la mesure et en prenant la moyenne des 
résultats obtenus. 

D'une manière générale l'erreur de lecture est systématique si la 
répétition de la mesure que l’on envisage implique la répétition exacte 
des mêmes lectures ; elle devient accidentelle dans le cas contraire. 

Bien entendu l’on ne peut pas toujours transformer l’erreur systé- 
maetique en erreur accidentelle. Cela ne serait d’ailleurs pas toujours 
désirable ; il importe en effet de s'attaquer directement à chaque cause 
d'erreur et de la corriger aussi exactement que possible. I1 ne doit 
rester, dans une bonne méthode de mesure, que des erreurstrop petites 
pour être corrigées ; chacune est négligeable. Mais leur réunion peut 
constituer une erreur appréciable. L'erreur accidentelle apparaît ainsi 
comme constituée par la réunion d'erreurs indépendantes et très petites. 

Ce résultat n’est pas encore assez précis pour constituer le point de 
départ d’une théorie mathématique. I1 faut préciser la relation entre 
ces erreurs élémentaires et leur résultante. L'hypothèse la plus simple 
consiste à supposer que ces erreurs élémentaires s'ajoutent ; c’est 
sur cette hypothèse que repose la théorie mathématique classique, 
qui conduit à cette conclusion que l’erreur accidentelle obéit à la loi 
de Gauss. Mais nous devons nous demander si elle est vérifiée en fait. 

Cela dépend naturellement des cas. Le calcul des probabilités nous 
habitue, dès ces débuts, à cette idée que certains cas semblent faits 
exprès pour donner l’occasion d'appliquer la théorie. C’est ainsi que 
a notion de probabilité a un sens très clair lorsqu'on parle du tirage 
d’une carte dans un jeu bien battu. I/étude de ces cas présente un 
grand intérêt, et aide à comprendre un peu ce qui se passe dans des 
cas plus complexes et dont la théorie ne peut être faite exactement ; 
les résultats obtenus dans les cas simples ne peuvent s'appliquer 
d’une manière exacte à ces cas complexes ; mais on ne doit pas être 
surpris de constater qu’ils s'appliquent d’une manière approchée. 
Ainsi on ne s'étonne pas de constater la constance de certaines fré- 
quences, même quand aucune théorie ne permet de les calculer. 

Dans la théorie des erreurs, un cas, où la théorie s’applique assez bien 
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est celui de l’erreur due à la réfraction atmosphérique dans les obser- 
vations astronomiques. Cette erreur comporte une partie systématique 
due à la densité croissante des couches atmosphériques traversées 
par le rayon lumineux. Cette erreur systématique corrigée, il reste 
une erreur accidentelle, due à l'influence des perturbations de l’atmos- 
phère ; elle apparaît bien comme la somme des déviations ainsi pro- 
duites aux différentes altitudes. 

Considérons maintenant le cas d’un télémètre, ayant pour objet la 
mesure d’une distance x, et la déduisant de la mesure d’un angle 6 
très petit et variant en raison inverse de x. La théorie s'applique assez 
bien à l'angle 8, et l’erreur d9 commise sur cet angle obéit à la loi de 
Gauss, dont nous rappellerons tout à l’heure la définition ; c’est une 
loi symétrique, donnant les mêmes probabilités à deux valeurs égales 
et opposées. Il en résulte pour l'erreur sur x une loi différente, évi- 
demment dissymétrique ; la théorie ne s’applique pas à l’erreur sur x. 
Mais si la mesure est précise, c’est-à-dire si l'erreur dx commise sur x 
(et non seulement chacune des erreurs partielles dont elle résulte) 
est très petite, 09 et dx sont proportionnels, et l’erreur sur x obéit 
à la loi de Gauss. 

D'une manière générale si une mesure est précise, c’est-à-dire si 
une variation de la quantité à mesurer x de même ordre de grandeur 
que l'erreur possible apparaît comme étant sans influence sur les con- 
ditions de la mesure, x et les différentes fonctions de x que l’on peut 
avoir à considérer varient proportionnellement dans ces limites, et 
il importe peu de savoir à laquelle de ces fonctions la théorie s'applique 
le mieux ; on peut l'appliquer à x et dire que l'erreur sur x obéit à 
la loi de Gauss. 

On peut exprimer la même idée sous une autre forme en disant que, 
les variables Ôx étant assez petites pour que l’on puisse négliger leurs 
carrés, toutes les formules deviennent linéaires. L'erreur totale ne 
peut donc être que la somme des erreurs partielles qui la constituent. 

Mais ce n’est pas là un raisonnement rigoureux ; même si les erreurs 
partielles sont très petites, d’autres formules que l'addition sont 
possibles. M. FRÉCHET a étudié le cas où l’on prendrait pour l'erreur 
résultante X la plus grande des erreurs partielles, prise en valeur 
absolue ; il a montré que dans cette hypothèse on peut développer 
une théorie analogue à la théorie classique ; mais la loi de Gauss 
est remplacée par une loi différente. 


» 
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Imaginons un examen dans lequel chaque candidat, ayant à effectuer 
un certain nombre de mesures, commettrait dans ces différentes me- 
sures des erreurs de valeurs absolues x, %2, .…, 4n. Supposons que l’exa- 
minateur dise : « Je ne crois pas au calcul des probabilités; peu m’im- 
portent la valeur moyenne ou la valeur quadratique moyenne de 
ces erreurs ; ce qui m'intéresse c’est la plus grande erreur que le can- 
didat risque de commettre. Je ne tiens donc compte que de la plus 
grande en valeur absolue des erreurs commises, et c’est d’après sa 
valeur X que je note le candidat. » Un tel examinateur aurait sans 
doute tort. Supposons quand même qu’il existe et cherchons à prévoir, 
avant l'examen, la note du candidat ; nous trouverons qu’elle dépend 
de Ja loi de M. FRÉCHET (du moins si # est assez grand). 

Mais un tel exemple est exceptionnel. Les fonctions que l’on peut 
rencontrer dans l’étude de la composition des erreurs sont le plus 
souvent des fonctions continues, à dérivées continues, et même 
développables en série de Tavror. $i les erreurs sont très petites, 
on peut négliger leurs carrés et leurs produits ; l’erreur totale est une 
fonction linéaire des erreurs partielles, et la théorie classique s'applique. 

Cette théorie, nous l’avons dit, nous conduira à ce résultat que 
l’erreur accidentelle obéit à la loi de GAUSs. Nous voyons comment 
il faut comprendre cet énoncé. Il repose sur certaines hypothèses, 
notamment sur l’hypothèse de l’additivité des erreurs partielles, et 
sur la petitesse aussi bien de chacune de ces erreurs partielles par 
rapport à leur somme que sur la petitesse de cette somme. Ces hypo- 
thèses semblent naturelles, et seront sans doute assez exactement 
vérifiées dans la plupart des applications. Mais il ne faut pas oublier 
que les résultats de la théorie ne sont applicables qu’à cette condition ; 
qu’il est nécessaire dans chaque application, soit de s'assurer que ces 
hypothèses sont bien vérifiées, soit de vérifier expérimentalement 
que la lo° de l’erreur est bien celle de GAUSS ; qu’en général il n’en sera 
ainsi que d’une manière approchée. Il ne faut pas oublier notamment 
que la théorie suppose que l’erreur accidentelle soit la somme d’une 
infinité d’erreurs indépendantes et infiniment petites ; or dans les 
cas les plus favorables il ne peut être question que d’un nombre très 
grand, mais fini, d'erreurs partielles. 


2. EXPOSÉ DE LA THÉORIE MATHÉMATIQUE. —- Rappelons d’abord 
la définition de la loi de Gauss. On d't qu’une var able £ obéit à la 


— 167 — 


PAUL LÉVY 


loi de Gauss réduite si la probabilité de l'inégalité 


(1) ae 
est, quels que soient + et G, égale à l'intégrale 
LEURS 
I PAR 7 
== € dE, 
(2) a). 


c'est-à-dire à l’aire hachurée de la figure ci-dessous, où l’on a tracé 
la courbe de GAUSS 


x? 
I 2 

0 

(3) Foie 
La valeur probable de £ est alors zéro, et sa valeur quadratique 
moyenne est 1 (propriétés caractéristiques des lois réduites). La loi 


4 


Fig. 2. 


de Gauss générale est celle à laquelle obéit la variable x — af, a étant 
une constante. 

Considérons alors une variable x, obéissant à une loi de probabilité 
quelconque, à cela près que sa valeur probable doit être nulle et que 
sa valeur quadratique moyenne a une valeur finie m. Si cette loi 
dépend d’un paramètre variable, nous dirons qu’elle tend vers 


x 


celle de GAUSS, ou d’une manière plus précise que la loi à laquelle 
obéit la variable réduite £ — à tend vers la loi de GAUSS réduite, si 


la probabilité de l'inégalité (1) tend vers la valeur (2) ; on voit aisément 
que cette convergence est nécessairement uniforme. 

La loi de GAUSS a une propriété importante, facile à démontrer : 
elle est séable. On entend par 1à que, si plusieurs variables indépen- 
dantes %,,%2, …., 4n, obéissent à la loi de GAUSS, il en est de même de 
leur somme 


(4) X = Xi + Notes + An. 


» 
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Mais ce qui est plus remarquable encore, c’est que, si les lois auxquelles 
obéissent x,, x, …, 4n sont quelconques, la loi résultante à laquelle 
obéit X ressemble plus à la loi de Gauss que les lois composantes, 
et tend pour # infini vers celle de Gauss, du moins si chacun des termes 
est très petit. En d’autres termes : la somme de n erreurs indépendantes 
et très petites, dont les valeurs probables sont nulles, obéit à une loi, 
tendant, pour n infini, vers celle de Gauss (1). 

C’est bien le résultat qu’il s’agit d'établir pour la théorie des erreurs. 
Ii faut d’ailleurs encore préciser le sens de l’expression « très petites » 
qui figure dans l’énoncé. Elle signifie que les erreurs sont très petites 
par rapport à la valeur quadratique moyenne M de l'erreur totale X ; 
rappelons que M est lié aux valeurs quadratiques moyennes ",, M, …, 
Mn par la formule 


(5) M° = mi + mi +. + me, 


Pour qu’on puisse dire que les erreurs sont très petites, il suffit alors 
évidemment qu'aucune d’entre elles ne puisse en aucun cas dépasser 
eM, © étant un nombre très petit, et tendant vers zéro pour # infini. 
Mais, la valeur quadratique moyenne d’une variable donnant assez 
bien idée de l’ordre de grandeur auquel il faut s'attendre pour cette 
variable, et les valeurs beaucoup plus grandes, si elles sont possibles, 
étant très peu probables, on peut se demander s’il ne suffirait pas que 
Ma, Ma, …, Mn, Soient inférieurs à M. Il n’en est pas ainsi ; mais cette 
hypothèse devient suffisante si l’on ajoute une condition supplémen- 
taire très peu restrictive, qu’on peut énoncer par exemple en disant 
qu'on commet une erreur relative négligeable sur la somme (5) en ne 
tenant compte que des valeurs de x,, %, …, x inférieures à «M. On 
peut donner aussi d’autres formes à cette condition supplémentaire, 
de sorte que le théorème fondamental est susceptible de recevoir 
plusieurs énoncés légèrement différents. Celui que nous venons d’indi- 
quer résulte des travaux de M. Lindeberg. 

Nous ne pouvons, dans cette conférence, donner la démonstration 
complète de ce théorème fondamental. Mais nous indiquerons briève- 
ment le principe de trois méthodes différentes auxquelles se rattachent 
les recherches relatives à ce sujet. 

(x) Sans la restriction que les x, aient leurs valeurs probables nulles, le résultat ‘subsisterait 
pour X, à une constante additive près ; mais si cette constante n’était pas nulle, elle devrait être 


rattachée à l'erreur systématique, et le résultat relatif à l’erreur accidentelle subsisterait sans 
modification. 
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La première méthode, celle de Bernoulli, repose sur l'étude du jeu 
de pile ou face. À chaque coup le gain x peut avoir les valeurs + I et 


DR: : 
— 1, chacune de ces valeurs ayant la probabilité; . Après n coups, 
les valeurs possibles du gain total X sont 


ñ, n — 2, REA, 


leurs probabilités respectives étant 
CO] I 
2 an RE ne 
La valeur quadratique moyenne de X, d’après la formule (5), étant 
Vn, il est naturel de considérer le gain réduit £ — — et, en désignant 


par dé = _ la différence de deux gains consécutifs, de mettre la 
probabilité de ce gain sous la forme fn(£)d£. Un calcul facile montre 
alors que /n(£) tend, uniformément, pour infini, vers la fonction de 


I 6? ’ er ARE pui E- 
GAUSS en A ; cela revient à dire que le gain réduit obéit, à la limite, 
T 


à la loi de Gauss réduite. 

On voit ainsi clairement la manière dont s’introduit cette courbe ; 
elle est dessinée en quelque sorte par la représentation graphique des 
coefficients du binôme. Seulement cela n’apparaîtrait pas si l’on mar- 


quait les points d’abscisses p et d’ordonnées C'(g — . + D) ; il faut 


en outre multiplier les abcisses par _ et les ordonnées par A: 

Ce résultat obtenu, il est facile d’abord d'en déduire une nouvelle 
démonstration du fait que la loi de Gauss est stable, ensuite de le 
généraliser de plusieurs manières, en supposant que la loi de probabilité 
de chacun des gains successifs 41, %o, …, An ne soit pas celle du jeu 
de pile ou face, et même qu'elle varie d’un coup à l’autre. Maïs il n’est 
pas facile d'arriver ainsi au degré de généralité nécessaire pour la 
théorie des erreurs, où l’on ne sait pour ainsi dire rien sur les lois des 
erreurs élémentaires, si ce n’est qu’elles sont petites ; une théorie 
introduisant d’autres hypothèses ne peut donc pas expliquer la grande 
généralité du fait que l’erreur accidentelle obéisse (au moins d’une 
manière approchée) à la loi de GAUSS. 


Passons donc à la seconde méthode, celle de M. Lindeberg (1920 
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et 1922). Chacune des erreurs partielles xn,(h — I, 2, …, n), étant 
représentée par nr, Mn étant sa valeur quadratique moyenne, de 


sorte que {x est l'erreur réduite correspondante, l'erreur totale réduite 
est 


| X mm ; Ma » Mn » 
(6) . Me RNA ni 


M. LINDEBERG considère d'autre part des variantes auxiliaires 
E, €, .…, en, Obéissant à la loi de GaAUSS réduite, et indépendantes les 
unes des autres. On voit aisément que la somme 


(7) L'an + Bee TEE 

obéit aussi à la loi de GAUSS réduite. Or on passe de la somme (7) à 
la somme (6) qu'il s’agit d'étudier en remplaçant successivement 
£,', par Z, %' par &, .…, £n' par 1. Une discussion des erreurs ainsi com- 
mises , relativement élémentaire quoiqu’elle nécessite quelques pré- 
cautions, montre que, sous les hypothèses rappelées tout à l’heure, 
non seulement chacune de ces erreurs, mais l’erreur totale commise 


en remplaçant = par = est très petite. À la limite, obéira bien à la 
loi de Gauss réduite. 

La troisième méthode dont nous voulons parler repose sur l'étude 
de la fonction caractéristique ; cette méthode, déjà considérée par 
CaucHY, et même incidemment par quelques prédécesseurs de 
CAUCHY, a été renouvelée par CHARLIER, et nous-même en avons 
développé systématiquement les conséquences. 

La fonction caractéristique, que nous désignerons par o(f), est la 
valeur probable de l’exponentielle imaginaire e””, x étant la variable 
obéissant à la loi étudiée, et £ une variable auxiliaire. Aïnsi, dans le 
cas d’une loi continue, où chaque intervalle dx a la probabilité f(x)dx, 
on a 


Pr + 00 
(8) e(0 = cf(x)dx 
et la formule de FOURIER 

| J + © Z 
Q) = [ea 
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nous permet de calculer f(x) connaissant (f). La fonction des proba- 
bilités totales 


se déduit alors aisément de la formule (9) par intégration sous le 
signe d'intégration. Il vient ainsi 


ÊTE ES —Ùx 
I I—e 
(10) F(x) — F(0) = = JL Nr o(t)dt. 
Or, à condition d'écrire cette formule sous la forme 
Sr nn 77 
PAR I —é 
(x1) F(x) — (0, — 2 der Je 3 «(é)dt, 


on a une formule applicable à toutes les lois de probabilité, continues 
ou non. Donc, dans tous les cas, une loi de probabilité est bien définie 
par la donnée de sa fonction caractéristique o(f) pour toutes les 
valeurs réelles de £. 
On voit d’ailleurs aisément que (0) — 1 ; que pour une loi réduite 
®'(0) = o et p''(0o) = — 1, de sorte que o(é) — r et W(f) — log o(t), 
2 


ont pour valeur principale — ; . Dans le cas de la loi de Gauss réduite, 


(4) = — £, et pour être sûr qu’une loi tend vers la loi de GAUSS réduite 


il suffit de montrer que W(f) + ; tend vers zéro, et cela uniformément 


dans tout intervalle fini. 


Considérons maintenant deux lois de probabilité indépendantes, 
de fonctions caractéristiques 


(12) ei(t) — V'opéitrs galt) = DATE 


0 


4, désignant la probabilité de la valeur x, pour la première loi, et 
Bx et y, définissant de même la seconde loi. La fonction caractéris- 
tique de la loi résultante, à laquelle obéit la somme x + y, est 


(13) g() = danpreitte + 3) À aff) qa(t). 


Donc, dans la composition de deux lois de probabilité (par addition 
des variables correspondantes, supposées indépendantes), les fonctions 
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caractéristiques se multiplient ; leurs logarithmes W'(f) s'ajoutent. 
Ce résultat n’est d’ailleurs pas spécial aux lois discontinues, comme la 


notation employée pour écrire les formules (12) et (13) pourrait le faire 
penser. 


Considérons alors # variables indépendantes x,, x, …, %n, obéissant 
à la même loi réduite ; nous poserons 


(x4) Log &(#) = 49 = — À [x + w(], 


o(é) étant la fonction caractéristique de cette loi ; on a w(o) — 0. Pour 
leur somme X, W({) est renrplacé par 


VE = 29) = — n À [1 + uw], 


et pour la somme réduite Lis le logarithme de la fonction caractéris- 


Vn 


(= 


2 
Ii tend uniformément vers— dans tout intervalle fini, de sorte qu’à 


tique est 


la limite Un obéit bien à la loi de GAUSS réduite. 


Nous venons de traiter le cas où les différents termes %, %2, .…., Xn 
obéissent à la même loi. Il suffirait de quelques lignes de plus pour 
traiter le cas général. On arrive ainsi à mettre en évidence que la 


/ . x ,e « . . PIE 2 
somme réduite ;, obéit à la loi de GAUSS, pourvu que soient vérifiées 


les trois conditions suivantes, toutes les trois essentielles : 

1° Chacune des erreurs x a une valeur quadratique moyenne finie " ; 
(la valeur probable de x est alors aussi finie et on peut, comme nous 
l’avons fait, la supposer nulle). 

2° Le plus grand des nombres m° est très petit par rapport à leur 
somme M? ; 

3° Les fonctions w(f), liées par la formule (14) aux fonctions carac- 
téristiques o(f) des lois auxquelles obéissent les variables réduites 
È— ï, sont toutes majorées, au moins dans un petit intervalle com- 
prenant la valeur £ — 0, par une fonction Q(i) nulle pout { — 0 et 
continue près de ce point. 
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Toutes ces conditions reviennent d’ailleurs à limiter le rôle des 


grandes valeurs des variables x, de sorte que la conclusion que TT: obéit 


à la loi de Gauss a bien la portée voulue au point de vue de la théorie 
des erteurs. 

Nous n’insisterons pas sur la troisième de ces conditions. La nécessité 
de la seconde est bien évidente : si l’une des erreurs partielles n’était 
pas très petite par rapport à l'erreur totale, il ne serait pas possible 
que la loi à laquelle obéit cette erreur fût sans influence sur l'erreur 
totale. Si les erreurs considérées sont en nombre infini, et se présentent 
comme les termes d’une série, il est alors essentiel que la série soit 


MÈ + Ms + + Ma 


divergente ; autrement chacun des termes serait une fraction 
appréciable de la somme. En cas de convergence de cette série, même 
si la série 

Mi + Ma His EM Eee 
est divergente, on ne doit pas s’attendre à trouver la loi de Gauss. 

À fortiori ne doit-on pas s’y attendre si les carrés des variables x 
n'ont pas leurs valeurs probables finies. Dans ces cas les grandes 
valeurs de ces variables ont des probabilités suffisantes pour influer 
d’une manière sensible sur la loi de probabilité de leur somme. Des 
circonstances très diverses sont alors possibles ;il peut arriver notam- 
ment qu’on obtienne à la limite une loi stable autre que celle de GAUSS ; 
il existe en effet une double infinité de lois stables. 

3. — CONCLUSIONS RELATIVES A LA COMPENSATION DES ERREURS. — 
Tels sont les caractères essentiels de la théorie mathématique qui 
permet de comprendre le rôle de la loi de GAUSS dans la théorie des 
erreurs. Terminons cet exposé par quelques remarques sur les méthodes 
employées pour la compensation des erreurs accidentelles ; on a beau- 
coup discuté la question de savoir si ces méthodes pouvaient se jus- 
tifier indépendamment de la loi de GAUss. Les conclusions de ces 
discussions sont les suivantes : 

10 S'il est vrai que l'erreur accidentelle obéisse à la loi de Gauss, 
ces méthodes sont les meilleures possibles ; toute autre méthode donne- 
rait une précision moins grande ; 

29 Si l'erreur n’obéit pas à la loi de GAUSS, mais n’en diffère pas trop 
ces méthodes conservent une valeur incontestable ; maïs on ne peut plus 
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affirmer qu’elles soient les meilleures possibles. Files restent bonnes 
tant que la loi de probabilité de l'erreur comporte une valeur quadra- 
tique moyenne, et permettent, au moins en théorie, d'augmenter 
indéfiniment la précision par la répétition des mesures ; 

30 Ces méthodes tombent au contraire complètement en défaut 
dans le cas des lois comportant une erreur quadratique moyenne 
infinie. De telles lois ne semblent guère pouvoir se rencontrer dans une 
bonne méthode de mesure. Mais s’il arrivait qu’une telle loi fût vérifiée, 
il importe de remarquer que dans ce cas l’on peut se tirer d’affaire 
par d’autres méthodes. Signalons notamment celle qui consiste à écarter 
dans des proportions déterminées les plus grandes et les plus petites 
valeurs trouvées, et à prendre la moyenne des mesures conservées. 
Cette méthode est applicable dans tous les cas et, même dans le cas 
où la méthode classique est meilleure, n’est pas beaucoup plus mau- 
vaise qu'elle. 


Sur quelques applications des équations intégrales 


PAR 
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1. — Je suis profondément ému par le grand honneur qui m'est 
fait. Je le suis d'autant plus que le grand nom auquel est consacré 
ce magnifique institut représente pour moi non seulement un beau 
chapitre de l’histoire de la science, maïs aussi une personnalité vivante : 
je me souviens parfaitement de Poincaré, de ses mouvements, de ses 
gestes, de sa voix et de l’impression profonde que produisait sur nous 
— auditeurs — sa pensée créatrice toujours active et toujours vivante, 
et de ia douleur que nous ressentîmes en apprenant un jour sa fin pré- 
maturée. 


2. — Je me souviens surtout du sentiment d'inquiétude scientifique 
dont étaient imprégnés ses derniers cours et les derniers mémoires 
qu’il publia. La crise actuelle de la science ou plutôt la crise de crois- 
sance de la science dont nous sommes tous témoins s’ébauchait déjà 
de ses jours. Il lui arrivait souvent de tenir des propos profondément 
sceptiques, parfois même pessimistes, mais il lui arrivait aussi souvent 
de prévoir le développement ultérieur de la science, et peut-être son 
scepticisme même provenait-il de ce qu’il comprenait parfaitement la 
fragilité des théories physiques, voyait clairement tous les coups de 
pouce nécessaires pour mettre une théorie en concordance avec l’expé- 
rience, et se rendait compte qu’on pouvait toujours sauver une théorie 
par un coup de pouce bien choisi et bien appliqué. 


V. A. KOSTITZIN 


3. — Dans mes conférences je m'occuperai précisément d’un coup 
de pouce donné à la mécanique et à la physique mathématique sous 
la forme de la physique héréditaire créée par l’éminent mathématicien 
qui fit récemment ici-même l'exposé magistral de la biologie mathé- 
matique. 

La physique héréditaire de M. Volterra est une tentative sinon de 
résoudre, du moins d’amoindrir les difficultés qui résultent de l’anti- 
nomie ancienne entre la physique ordinaire et la physique moléculaire, 
entre la mécanique des corps bruts et la mécanique statistique. Depuis 
l’époque de LAPLACE avec son idéal scientifique (qui fut aussi l'idéal 
scientifique de son temps), avec sa conception hardie d’un vaste ésprit 
capable de déduire le passé et l’avenir de l’univers de son présent, la 
science a subi beaucoup de changements, et bien de prémisses du 
raisonnement de LAPLACE se sont volatilisées. Le point matériel a 
ondoyé. On a perdu l’habitude de se servir de la notion de force sans 
réfléchir sur sa nature. Les équations différentielles ne sont plus ni la 
forme unique ni la meilleure forme pour exprimer les phénomènes, de 
même que les données initiales de la mécanique de LAGRANGE ne sont 
plus du tout suffisantes pour déduire le passé et l’avenir d’un système 
ou de l’univers entier. Et la conception mécaniste du commencement 
du xIx® siècle ne nous apparaît actuellement que comme une extrapo- 
lation hasardeuse et injustifiée des lois souvent encore douteuses du 
monde macroscopique dans le monde microscopique. 


4. — Dans ces conditions la tentative de totaliser les traces des 
actions subies par un système, la tentative d’escompter l’histoire du 
système étudié n’est pas en elle-même plus paradoxale que la totali- 
sation des influences extérieures sur un élément. Ici, comme 1à, on 
opère avec des notions d’origine statistique, on fait des hypothèses 
souvent grossières sur la nature des actions intra-moléculaires, on 
construit des modèles mécaniques en utilisant pour ce travail délicat 
des matériaux hétéroclites. Dans les deux cas, en raisonnant on se 
base sans s’en rendre compte sur des prémisses plutôt intuitives. 
Il n’y a là-dedans rien qui puisse nous choquer : il est très {naturel 
d'utiliser dans notre raisonnement les éléments immédiats de notre 
expérience. Progressivement ces images grossières sont remplacées 
par des notions plus fines, plus élaborées. Des lois plus exactes se 
dégagent peu peu à des hypothèses approximatives. C’est une évolu- 
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tion par laquelle la physique moléculaire a déjà passée, et ses succès 
mêmes permettent à la physique héréditaire d'évoluer plus rapidement. 

Serait-ce une méthode transitoire qui s'explique par la faiblesse de 
notre technique ? Je ne le crois pas. Il ne faut pas oublier que les choses 
temporaires sont les plus durables. La physique héréditaire n’est pas 
un stade dans le développement de la science -— c’est une méthode 
qui fera toujours partie de notre outillage scientifique par suite de la 
nature même de notre connaissance. Pour étudier un système, nous 
serons toujours obligés de 1e séparer du reste de l’univers et de chercher 
ensuite un équivalent aux actions du reste de l’univers. Nous serons 
toujours obligés de substituer aux interactions des particules de la 
matière des moyennes observables et de chercher ensuite l'équivalent 
des actions négligées. Ces équivalents peuvent se présenter sous une 
double forme : 1° les forces fictives et les conditions limites introduites 
dans jes équations, 29 l’histoire du système. Certes, c’est un artifice, 
mais c’est un artifice naturel et inévitable. 


5. — En mettant en équations n'importe quel problème de la 
philosophie naturelle on cherche toujours la forme la plus simple et 
on choisit presque toujours la forme linéaire, celle de HooKkE. Tout en 
comprenant les raisons de ce choix je veux signaler quand même que 
pareille simplification va presque toujours à l'encontre des intérêts de 
la recherche scientifique et que le temps est venu de nous débarrasser 
peu à peu de cette habitude. On peut en effet citer de nombreux cas 
où la forme linéaire ne donne pas et ne peut donner même la pre- 
mière approximation, la fonction n'étant pas développable en série 
de TAYLOR. On peut citer d’autres cas où l’on perd en première approxi- 
mation une propriété essentielle du problème comme l'interaction 
des individus dans la biologie dynamique ou l'interaction des particules 
dans un grand nombre des problèmes physiques. Il sera sans doute 
nécessaire dans l’avenir de passer à l'étude des équations non linéaires, 
et alors à côté des difficultés analytiques on aura des difficultés philo- 
sophiques. Pour en donner l’idée je me permets d'indiquer rapidement 
un exemple intéressant. 

Prenons une équation intégrale non linéaire 


g(x) = f(x) +? Î El t)Ef+(#)]dé, 
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K{(x, t) étant pour fixer les idées un noyau bien connu 


K(x, t) = x(1 — 6) x Lt 
K(x, t) = Wir — x) Ph 


et supposons les fonctions F(+) et f(x) développables en série de TAYLOR. 
Il est très naturel d'arrêter ces développements aux premiers termes. 
Ecrivons donc 


I 
Q(X) = 5 + X + à [ K{(x, tg(t) + mp°(8) dt. 
Le] 


Cette équation intégrale se résout par des fonctions élliptiques de 
JACoB1, et pour chaque paire de valeurs À, 14 on a une multiplicité et 
parfois une infinité de solutions alors que l’équation linéaire corres- 
pondante se conduit suivant les règles de FREDHOLM, et d’autre part 
en choisissant convenablement la fonction F(c) on peut assurer l’uni- 
cité de solution de l’équation initiale. Nous allons voir que la multi- 
plicité de solutions se présente aussi pour quelques classes d'équations 
linéaires, et une question très importante se poserait alors sur le 
rapport de ces solutions multiples avec la réalité. 


6. — Après ces préliminaires nous allons nous occuper d'équations 
intégrales linéaires de la physique héréditaire. Il en existe trois formes 
X 
G) 4) = 0 22 fl ec — net 
| e 
(2) p(x) = f(x) + à [ ÉRN tye(t)dt 
X—hN 
(3) o(x) = f(x) + À 1 gx — tjo(t)dt. 


Dans le premier cas on suppose que l’histoire commence à l'instant f, ; 
dans le second cas on suppose que l’action héréditaire appréciable 
est bornée à un petit intervalle de temps, et dans le troisième cas 
on considère toute l’histoire du système. Nous allons étudier surtout 
l'équation (2). 

Il faut remarquer tout d’abord que sa forme nous permet de la 
rapprocher des équations aux différences finies. On peut donc 
s'attendre à ce qu’elle en ait toutes les particularités et à ce qu’on 
puisse y appliquer les mêmes procédés de résolution. D'autre part 
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c’est une équation intégrale linéaire et on y peut appliquer les méthodes 
servant généralement à leur résolution et tout d’abord la méthode des 
approximations successives. Nous verrons par la suite que les deux 
méthodes se complètent mutuellement. 

Désignons l'opération 


de g(x — t)p(t)dt 
X—h 


par le symbole Gs et écrivons l'équation (2) sous la forme 
p = f + Gy. 
En y appliquant # fois l’opération G on obtient ; 


HT 
p—=f+ Ÿ NC} + Ge 
Rk=—=1 


ou, ce qui revient au même 
X 
X—n 


et = do + Da [7 ete — aie +20 [7 ete — arte 
2; X—) e 


Il est facile de voir que les fonctions Gn(y) sont déterminées par 
des relations recurrentes 


&n(y) — de g(y — U)gnni(u)du  (n—1}h<Ly<Lnh 
LOS de ; g(y — w)gni(w)du hR<y<L(n— 1h 
0) = [lang (ide o<y<x 

£n(y) — 0, en dehors de ces limites. Remarquons que chaque 


nouvelle itération fait reculer la limite inférieure de l’intégration ; 
donc, si nous voulons obtenir une solution par cette méthode nous 
devons supposer que l’histoire nous est connue depuis — ©. 

Passons rapidement en revue les résultats de l'application de 
l’opération G aux diverses fonctions élémentaires : 


a) JON =" ; Cr an 
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en désignant par km l'intégrale 


Gta) = ko (x — au + né UEske — Bi?) 
à) FE) =; mi Grf(x) — 6" ÿ(a), 


en désignant par y(x) l'intégrale 
‘h 
yla) = [ pales de 
/ O 
B 
o je = [votes 
A 
B 
Gf(x) =f V(Pe ln (at)dt. 
A 
On déduit facilement des formules précédentes 


3 nh 
ju Ln(x — 2)dz — [ AE ER Te 
e/ X n 0 


°nh 


16 gr(a)e "dr = (a). 


Ecrivons maintenant la formule 


n 
p—f+ Der + AHGnHE, 


ÉT 
sous la forme 
x de a% 
pÜx) = f(x) + [ ajde V'X-to,(x — 2) + H Re 
; en 2 FER Rss (NX —2)o(2)d2. 
La série 
T(Â; x — 2) — > M1 g(x — 2) 
R=1 
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vérifie l'équation intégrale 


TO; x — 2) — g(x — 2) +3 [ gx — u)TQ ; u — 2)du[oL<x—2<h] 


x 
M; #— 9 = f à 
X— 


On voit donc que pouro<x—z<h 
; TO; x — 2) = RA; x — 2), 
en désignant par R(À ; x —- z) le noyau résolvant ordinaire de VoL.- 


TERRA. 
On peut maintenant établir le caractère de convergence de la série 


elx) = fn + DC»). 
n=—=I 
Supposons par exemple que g(z) n’est jamais négatif et que la fonction 
f(x) est bornée 
RCRESUE 


On obtient alors la série majorante 


CO 
| n EE M 
M+ ŸYMIAl# = TE 
N—I 
On voit que la série est convergente tant que | À | &, < 1 et diverge 
dans le cas contraire. La constante | À | k, joue un rôle important, 
et nous l’appeilerons mesure d’hérédité. 
Supposons ensuite que f(x) croit comme e°%* : f(x) sv e%*; on a en 
ce cas 
Grf(x) 2 e7" (a), 
et la série majorante a la forme 
ax 
+ VA fre yn(a) — € 
D a RAIRIEE 
Ces exemples montrent que dans un nombre de cas Ia série 


X 


G) er) = ft) + D rGrf(x) = (9 + sh (A; x — 2)f{2)dz 


—— C 
R==T 
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peut nous donner une solution de l’équation (2) mais il est très facile 
de s’assurer que la formule (3) n’en donne pas la solution complète, 
et qu’il faut examiner l’équation homogène. 

Nous avons indiqué l’analogie qui existe entre les équations aux 
différences finies et l'équation (2). Or, précisément de la même façon 
qui permet d'obtenir l'équation caractéristique 


1—eh—0 
de l'équation fonctionnelle des fonctions périodiques 
fx + h) — f(x) = 0 


et de former le système complet des fonctions trigonométriques, 
on peut construire l'équation caractéristique 


h 
(4) 1e) à [| ete vds 


de l’équation intégrale homogène 


X 
(5) (x) — » [ dure t)p(#)dt. 
Pa 
L’équation (4) a en général une infinité de solutions an ce qui nous 
permet d'utiliser le système de fonctions 
ea, ns 
pour développer sous certaines conditions en série une fonction arbi- 
traire dans l'intervalle de mesure }. 
Supposons maintenant que l’on ait 


À Lo, g(2) = 0. 


Dans ce cas l'équation (4) a une solution réelle unique «,. Cette solution 
est positive si la mesure d’hérédité Àk, > 1 et négative dans le cas 
contraire. Quant aux solutions complexes, on les trouve en résolvant 


le système 
, h 
| 1= À ji g(z)e cos nzdz 
[e) 


vh 
| OX | g(a)e* cos n2d£. 
v Oo 


(6) 


SUR QUELQUES APPLICATIONS DES ÉQUATIONS INTÉGRALES 


Remarquons que pour g(z) décroissant la partie réelle de & = £ + in 
est toujours négative. D'autre part dans ce cas il n’y a jamais de solu- 
tions purement imaginaires. Donc si la mesure d’hérédité est inférieure 
à l’unité la solution de l’équation (2) est toujours stable. Or, à vrai 
dire, nous n’avons guère de raisons pour supposer dans tous les cas 
g(z) décroissant, et alors il faut imposer des conditions supplémentaires 
pour assurer la stabilité de la solution de l'équation (2). 

On peut maintenant obtenir une solution complète de l'équation (2) 


D #0 = + DOG + set + D Aer LA enr 
en ajoutant à la série (3) le développement 
AE + D Autre + Aer, 


Les coefficients A: sont soumis à une seule condition — la convergence 
de la série. On voit également que la fonction f (x) dans ce cas ne 
peut être choisie arbitrairement, tout comme dans le troisième théo- 


x 


rème de FREDHOLM, à cette différence près qu'ici le caractère de 
croissance de f (x) tient place de la condition d’orthogonalité 


1 ex(x)f(x)dx = 0. 


Supposons en effet que f (x) croisse comme Aezx. Il est facile de 
voir que le coefficient « doit satisfaire à la condition 
Y(4) < Yo) 
d’où résulte l'inégalité 
CI CTS 
Si cette condition n’est pas vérifiée, la série (3) diverge. 


Nous terminerons ce paragraphe par un exemple. Soit g(z) — ex, 
Le système caractéristique prend la forme 


E— — 16 ET cos n# 
n = he ET Sin n4 
où 
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Il est facile de voir que les courbes (8) ont un nombre infini de 
points d’intersections. 


7. — Existe-t-il quelque chose d’analogue au premier théorème de 
FRrEepHoLM ? Le spectre de l'équation (2) est continu, l’intervalle 
de l'intégration est infini ; donc, il faut chercher la solution de (2) 
sous forme d’une intégrale de FOURIER. Supposons donc f(x) repré- 
sentable sous la forme 


a) & À *e HE [ F À frpéiale — De 
27 © : 


= ‘0 — 


Nous obtenons en itérant 
5 LR Se s , 


ce qui nous permer d'écrire une représentation formelle de o(x) 


+ © ?+ 0 
I da jalx — Ë 
AIS de I — fe) | HOeste — Sd 


— Cds; 
ou 

D D ip ee 

(x) — (62) Von Ja — ‘er f(E)eialx — S)dE, 

Or justement nous avons vu que pour g(z) décroissant l'équation 
caractéristique n’a pas de solutions purement imaginaires. Donc le 
dénominateur 1 —}y(ix) ne peut s’annuler que pour « = 0 ; On voit 
qu'ici aussi la mesure d’hérédité joue un rôle essentiel. Il faut remar- 
quer en outre que pour l'application de l'intégrale de FOURIER il 
est nécessaire de connaître la fonction f(x) non seulement dans le 
passé, mais aussi dans l'avenir. 


8. — Nous avons donné la solution L(x) de l'équation homogène (5) 
sous forme d’une série. Or, il est possible de lui donner la forme d’une 
intégrale en supposant L(x) connue dans un intervalle de mesure 4, 
par exemple dans l'intervalle (0, 4). En désignant par R(1 ; x — w) 
le noyau résolvant ordinaire de VoLrERRA et en supposant 4 <x<2h 
on obtient la solution de l'équation intégrale 


x h 
p(x) = À J g(x — to(t)dt + À J + g(x — t)V(t)at. 
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sous la forme, 
h 
= [HO 224086 
O 


La fonction H, est définie par les relations 


rh +3 
HER 2) À | RAA; x — u\g{u — 2)du (O<L3<x%x — h) 
/h 


TX 
Hi; 4,2) — gix — 2) +2 | RAA; x — u)g(u — 2)du (X—h<2<h). 


Il est facile de voir qu’on a en général pour nh < x < (n + 1)h 


?h 
o{x) = À | H,4 ; x2)V(2)dz, 


PA) 


HA étant définie par la relation recurrente 
h « 
Hal; x,2) = f HA 3 9 + (0 — 2) JE; # — (n = 1)h, Yéy. 
0 


On peut donc calculer la fonction p(x) de proche en proche. La fonction 
p(x) ainsi calculée est en général discontinue pour x — k ; au contraire, 
la continuité est assurée pour x — nh (n > 1). La possibilité de calculer 
o(x) en partant de ses valeurs dans l’intervalle (0, k) peut servir de 
preuve indirecte de la fermeture du système de fonctions [ex]. 


9. — Il est facile de.s’assurer que cette particularité n’est pas 
spéciale aux équations (2) et qu’on la retrouve dans les équations 
de FREÉDHOLM et de VOLTERRA aux limites infinies. 

I1 est toujours possible en effet de représenter la fonction g(x) 
sous forme d’une série de fonctions d'HERMITE s’il s’agit d’équation 
de FREDHOLM, ou d’autres fonctions analogues, s’il s’agit d’équation 
de VoLIERRA. En arrêtant ces développements à un terme quelconque 
on peut obtenir une fonction continue suffisamment approchée de 
g(x). Cette fonction utilisée comme noyau de FREDHOLM ou de VOL- 
TERRA jouira des mêmes propriétés que ‘g(x). Un petit exemple suffira : 
soit 
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noyau de l'équation de FREDHOLM 


rie 2 ' 
px) — À | en — 2 p(z)dz. 


&7 =—:@ 


L'équation caractéristique à la forme 


+ © Dh, 
12 | Eire. 
— oo 


En posant « — Ë + 1n on obtient 


E2 — 1? 


ën s. 
4 touts 
e cos > VE 
ue ë. 
CSI > —= © 
Il en résulte 
En — his 0 
62 — n° ñn 
ÉORE * : om 
Er MR 


On voit donc que pour À < 0 il faut choisir # impair alors que pour 
150 n doit être pair. Désignons dans les deux cas la quantité positive 
ee par e4, Nous aurons 

T 

Le n°? — 4A 
LE 4AE? — arr? — 0 
E — 2A + 2VA? + rm. 


On trouve donc pour chaque À un système de fonctions [e**] tout 
à fait analogue à celles déjà considérées dans le cas de l'équation (5). 


10. — Nous avions admis jusqu’à présent l'hypothèse de l’absence 
du moment initial. Si l’on suppose avoir fait table rase de tout ce qui 
a précédé un moment donné {,, on peut obtenir une solution et une 
solution unique de l'équation intégrale 


x 


(x) = j(x) + À J gx — bof x>h+t 


x—h 
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laquelle pour 4, < x < h + 1, doit nécessairement prendre la forme 
d’une équation ordinaire de VOLTERRA 


a) a = 1) 2. [l'atx — toit 


Remarquons pourtant qu’on doit être bien sûr que l’histoire com- 
mence au moment 4, ; on peut toujours se demander s’il n’y a pas eu 
d’influences cachées de l’histoire précédente, d’influences pouvant 

| affecter la solution cherchée. En supposant que tel n’est pas le cas, 
| on trouve 


#4 
#2) = 16) + L R(;,2 — ujf{ujdu (Gb <x<b + h) 
À 
et on procède ensuite par les prolongements successifs. On a pour 
b+h<:< th + 2h 


gd 


c 1 o+h 
ax) = f(x) + À J gs — tjaléjdt + } À g(x — tje(t)at. 


L — 


Or 
o+h ty+h 1 
: EN : giz — tje(fdt — [ : gx — tdKflt) + à J R(à; { — u){(u)du] 
| 04h 
L _ J P,(2 ; x, 2)f(2)dz, 
1 Lo 


en désignant par P,(;; x,z) la fonction 
to+h 
P,0:4,27 = gx — 2) +7 4h gx — DR ;t— 2jdt (x—h<z<i,+h) 
Z 
h415 
P,(À: #,2) =) FL gx — DR( ; À — zdi (by 2 < % —- h). 
4h 
Cela nous donne 
o+h 

ex) = 1) +2 É P,0 ; 3,2)f\)dz + 

, to+h 

+ } 16 RO; 3 — 2ddAf(t) + À SL P, (2; 4, u)f(u)du]. 

5 Lo+h to 


On peut continuer ainsi indéfiniment en obtenant toujours une 
solution unique. 
rs 189 = 
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11. — Voyons maintenant la théorie des équations intégro-diffé- 
rentielles ; commençons notre étude par l’équation linéaire homogène 


(0) gl) + Brelx) = 1Ge(x), 
qui ne diffère de l'équation ordinaire des mouvements harmoniques 


que par la présence du terme héréditaire. L’équation caractéristique 
a la forme 


h 
(10) a? + P? — Ay(a) — À [ BCE Fdz. 
« oO 
En supposant À > 0, g(z) > 0 on voit que Ày(+) est une fonction dé- 
croissante de + co à O ; d’autre part sa croissance pour & —>— © 
est beaucoup plus rapide que celle de &? + f?. Les courbes 


(D) ES vol 

(IT) Vi Ay(a) 

ont toujours un point d’intersection. À mesure que À diminue ce 
point se déplace vers la gauche, et pour À — O0 on à 4, — — ©. D'autre 


part à, est positif pour 18, > BF? et négatif dans le cas contraire ; il 
est nul pour Àk, — (?. Posons ensuite «x — Ë + in. On trouve 


| 2 + BE — x [ gene cos n242 
[M e) 


(x1) \ h 
DÈT = — À 1 gene sin nzdz. 
[e) 


On voit sans peine que pour g(z) décroissant £ < 0 et qu’il n’y a pas 
de solutions purement imaginaires. Donc, en ce qui concerne les solu- 
tions complexes, la stabilité de la solution complète est assurée. Mais 
le terme A,eo% peut présenter des difficultés, surtout quand {2 est 
petit, car alors la condition d’instabilité 


M > 8 


peut facilement être remplie. 
Parmi les solutions de l’équation (10) il y a tout d’abord deux 
développables en série de puissance À 


Le] 


a — D Yrhr, 
n—0o 
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On trouve sans peine les valeurs des premiers coefficients », : 
Y = 1 


. le 
— de V0 dz 
Y1 2% Je g(a)e 


2 k 
= — era d [ zg(2)e T®dz 
 O 


Ces solutions que nous désignons par w et » correspondent aux 
solutions de l'équation caractéristique 2? + {2 — o sans hérédité. 
Mais il y a d’autres solutions en nombre infini on — Ën + 2m et 
pouvant servir à former un système complet de fonctions 


apX CAE 
0 CE ere 


e : cu 


On peut montrer que %n —> — © pour À —+ 0. Ces fonctions peuvent être 
utilisées pour construire une fonction 


V2) = Aer + D + Aer 


vérifiant l’équation (9). 


12. — Ia résolution de l’équation avec le second membre 
(72) gx) + Po(x) — 1G4(x) = (x) 
dans le cas où le moment initial #, est connu et où l’on fait table rase 
de l’histoire antérieure ne présente pas de difficultés et n’entraîne pas 
d’ambiguité. On peut facilement par un procédé dû à M. LaLEsco 
par exemple, transformer l'équation 
X 


(zzv) (x) + P'o(x) = f(x) + X j g(x — 2) par ( Lx LE + h) 


lo 
en équation intégrale ordinaire de VOLTERRA, on peut ensuite prolonger 
la solution obtenue dans les intervalles suivants. Multiplions en effet 


PU — x) 
F 


l'équation (12 bis) par sin dx et intégrons par rapport à 


x de , à t ; nous obtenons 


p() = p(éo) cos B(É — &) + pb) 
a RE dax fus (4 ETES dx 
l to 


sin u 4 ee 


—- 
0 
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La résolution de cette équation se fait sans difficultés. On écrit 
ensuite pour l'intervalle suivant 4, + h << x < 4 + 2h 


%0+h 
g(x — à) g(édt + À [ glx — ty e(b)dt. 
h 


to+h e/ X— 


[x 


gx) + Be(x) = f(x) + à 


On substitue dans la dernière intégrale la solution obtenue pour 
bh < x Lt, + h et on procède de la même manière que dans l'intervalle 
précédent. En désignant par gn(x) la valeur de w(x) dans l'intervalle 
b+nh<x<i + (n +1)h, on aura à résoudre dans chaque inter- 
valle l'équation de VOLTERRA 

En(X) = Pn_i(l) + nh) cos B(Xx — À, — nh) + 
sin — 1, —nh 


+ v'n1féo + nh) 


: , #2 to+nh 
LP qe ir | Et uen alidiit + 
1 


tot+nh É 

x % in Olx — f\olt — 
+à [7 mad [SRE à, 

to+nh u Î 


On voit que la solution dépendra de deux constantes o(4), o'(4) 
et que le problème se résout complètement. Mais il reste toujours un 
doute sur notre droit de négliger l’histoire antérieure au moment #. En 
effet, pour l'instabilité de la solution deux conditions sont suffisantes : 

10 l'existence de traces insignifiantes de l’ancien état du système ; 

29 la mesure d’hérédité dépassant £?, ce qui est bien possible quand 
F? est petit. 

On a alors à, > 0, et le terme séculaire A,e%* peut toujours 
apparaître à côté de la solution régulière de l'équation (12). Or, jus- 
tement ces conditions sont facilement réalisables. Si l’on étudie par 
exemple l’hystérésis magnétique on prend soin de faire passer l’échan- 
tillon étudié par un chauffage suffisant pour faire disparaître les traces 
de l’état magnétique antérieur. Et malgré cela peut on être sûr qu’il 
n'en reste pas de traces insignifiantes ? Est-on sûr que la mesure 
d’hérédité est toujours suffisamment petite ? 


18. — Supposons maintenant qu'on cherche une solution de (12) 
valable pour tout l’espace et que l’on aie 


ï +00 +00 s - 
= [Ua ['pée-te 
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Muiltiplions l'équation (12) par eux) et intégrons par rapport à 


x de — à + ©. Nous obtiendrons ainsi 


+00 +00 
J esltu-2)@"(x)dx + P? jl el) (x)dx — 
ne) —X 
Lo +00 
— } | éurzldx F & g(x — to(t)di al ex x) f(x)dx. 
© — D X—h —0 


Admettons que 


On a 


+ RS 
e etslux)@ Lee NY — — — ? f enrltu-alo(x)dx. 


0 — co 


D'autre part 


—+ 00 de 00 th 
. eux) dx : g(x er t)e(t)dt =} eo | este) p(x — t)dx = 
—c0 x—h —00 l 
+00 
— (24) de ptet di, 


On a donc 


Lo ARE UC (u—t) 4 
fx o(é JE u—i) _ + 
La _ pP—2— (4) 
et on obtient sans peine 


+00 d +-00 ; 
(13) o(x) — _ 12 Dr — À _ Xy(62) fa f(t)err 9 di, 


00 0 


Pour avoir la solution complète de (12) il faut joindre à l'expression 
(13) la solution de l'équation homogène, la solution qui comporte. 
comme nous l’avons vu une fonction et deux constantes arbitraires 
Cette solution a pourtant un inconvénient déjà signalé : elle suppose 
la connaissance de l’avenir alors qu’en réalité c’est sur le passé seul 
qu’il nous est permis de nous baser. 


14. — Pour éviter cet inconvénient multiplions l'équation (12) 


sin Ou — : 
par = 5, et intégrons par rapport à x de w — k à u. On obtient 
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ainsi l'équation 


| o(u) — œ(u — h) cos 8h + g'(u — k) _ Ph 

(14) Pn. ; 
| “re [ sin “ — x) f{x)dx + À d S(u — z)p(2)dz 
\ u—h os 


en désignant par $ la fonction 


h . 
Sup [MED nd (u— 2h <Lz<u— h) 


ù 


Su—2)— |" M  oe (u—h< 2<u). 
[e] 

I1 paraît au premier abord que pour résoudre l'équation (14) il 
nous faut connaître la fonction o(x) dans un intervalle de mesure 24. 
Or on voit sans peine qu’il suffit de connaître o(x) dans un intervalle 
de mesure k. En effet, supposons h << 4 << 2h et o(x) connu dans l’in- 
tervalle (0, k). Transcrivons l'équation (14) sous forme 


o(u) = g(u — h) cos £h + &'(u — h) EE + 


—- [2 {(x)dx + à [°stu — 2)o(2)dz + 


rh Po 
+ À [ S(u — z2)o(2)dz + J S(u — 2)p(2)4z. 
«/9 u—2h 
Or il est très facile de déterminer o(x) pour — h < x < o quand 
il est connu dans l'intervalle (0, k). Il suffit de différentier l'équation 


[e) x 


g(x — z2)e(2)dz + x | g(x — z)o(z)dz. 


[e) 


#2) + Bet) = 16) +2 [ 

V/ X—h 

Si gW(x) est la première dérivée ne s’annulant pas pour x = #, 

on obtient après k différentiations une équation intégrale de VOLTERRA 

et en la résolvant on trouve o(x) pour — k > x < 0. Le reste s'effectue 

sans difficulté. Dans le cas où la différentiation n’est pas possible 

ou bien g(4) — o quel que soit À on peut utiliser le procédé du 
n°12: 
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15. — On peut traiter de la même manière les équations intégro- 
différentielles de n'importe quel ordre. Soit 


G3) Date) = je +2 [ste — 9 D 


R=O R—0 


une équation intégro-différentielle aux coefficients constants. En 
la multipliant par eWw-*Jx et en intégrant par rapport à x de 
— ©œ à + © on obtient 


+ 00 
E %k ‘7 eut) (x)dx — J gvltu-2)f(x)dx + 


re) 


+00 ce 
in DE A emnrn te 


X— 


+ 00 PH 00 
j gite) M) (x)dx — (iv) l evtu-xo(x)dx 
— 


en supposant naturellement que 
QE œ) = g(Eoœ) — = NE) = 0. 
De même 


+0 x 
J enraidn | gt — 2et(ade = 


— 0 e X—h 


"+ 00 nz+h 
— ] eW)(z)dz J g(x — 2e dx — 


LES 3 


+ co 


] 


h 
el (z)eutu—s)dz f gx) vedx — 
ss 


dy —2 


A +00 
— (iy)y(dy) J cute (z)dz. 


On obtient donc après la substitution et l’inversion 


T + co 11, | 
(16) «eZ CE iy)B(6y) le à 
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en désignant par A(é) et B(é) les polynômes 
0 


A) = > ayt* 


k—0o 


DONS 


k=0o 
L'équation caractéristique a la forme 
A(a) = XB(a)y(2). 


On peut diviser ses solutions en deux groupes : 1° # solutions x 
tendant pour À —+ o vers les racines de l’équation A(x) = 0 ; 20 solu- 
tions «x tendant vers — © pour À — 0. Les premières correspondent 
aux solutions ordinaires des équations différentielles, les secondes 
permettent de construire un système complet de fonctions ez»*, 
Pour avoir la solution complète de (15) il faut donc adjoindre à 
l'expression (16) une somme de termes complémentaires 


ñ _ un 
NB dE Y Carr + One 
kR=—1 


k=0o 


Je n’insiste pas sur les conditions de validité de ces formules. 


16. Ces résultats se prêtent à une généralisation. En effet, il n’y a 
pas de raisons pour supposer que la fonction d’hérédité est la même 
pour toutes les dérivées de o(x). Dans les problèmes de magnétisme 
de fer par exemple on peut lier le terme héréditaire d’hystérésis à 
la première dérivée de o(x) et le terme héréditaire de trainage à la 
fonction + elle-même ; on aura alors l'équation 


ait) = (x) Le ïr a. — 2)p(z)dz + eh Us — 2)v'(2)dz. 


Il est donc tout à fait naturel d'écrire l'équation intégro-différen- 
tielle sous forme 


x 
(17) ae) = 100) + D | Eux — 2er 
DE ae X— Ax 
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En la multipliant par ex et en intégrant par rapport à 
x de — © à + o nous obtenons 


+00 
ÿ d EN eva) (x)dx = f due) f(x)dx + 


Rk=o 


Le 


m 


+00 x 
+ D je : euta)gx x ù ga(x — 2)9%)(2)dz. 


k=0o ei 
Nous avons vu déjà que 


+o net 
si évtu-2)e 0 (x)dx = a [| eut-s)e(x)dx. 
— 


— 
D'autre part 
+00 x ) 
[ and [7 gute— apnée = puisnti [7 etoevidr 
e x—hy . CH se) 


OR 


en désignant par yx(iy) la fonction 


he 
va(2y) = JL gu(x)e dx. 
[e) 


On obtient donc après substitution et inversion 


+00 
(18) at = ZE [ De fre eve—2){(ndz 


en désignant par D(iy) la fonction 


D(iy) — A(iy) — 2 va(2y) (y). 


É=0 


Quant à l’équation caractéristique, elle a la forme 


D(a) = A(a)— D na) = 0, 


et sa discussion donnerait des résultats analogues à ceux du numéro 
précédent. 
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17. Nous allons ébaucher l'étude des systèmes d'équations intégro- 
différentielles. Soit 


(x9) gi (x) + ÿ aupr(x) = f(x) + > fs ga(x — t)x(t)dt 


R=—1 


un système de ce genre. Supposons qu’on puisse résoudre le système 
homogène correspondant par des fonctions 


eo) = DA 


Cela nous donne immédiatement 


ñn 


n 
Agus + Daudu; = D Aurult) 
k=1 CE 


en désignant 


k 
ee 1 gutt)e Wat. 


L'élimination de coefficients Ax; ; nous donne l’équation caracté- 
ristique 
+ au — Yu(n), do — Yu), *-- 
dpi — Yu), LU? + Gp — Yaa(u), ++ | = 0. 


clé lois sols else rares ST Se Sas nu se a au ia he 


A(b) — 


L'étude de cette équation est assez pénible ; on peut pourtant 
s'attendre à ce qu'il y ait en général deux groupes de solutions : 
19 2n solutions correspondant aux solutions ordinaires de l'équation 
caractéristique dans le cas d’absence d’hérédité ; 2° solutions per- 
mettant de construire un système complet de fonctions eun*, 

Pour obtenir une solution du système (19) on peut appliquer la 
méthode de FoURIER. Multiplions (19) par e*{"-* 4x et intégrons par 
rapport à x de —o à ©. Nous obtenons 


+00 
M gi(x)eist vegas + Ya ” ni g(rert-dx = 


sf" x)eïx| u—X) fx + Y oh uaax | Lix(X SU t'px(t)dt 
Me ER 


or mé 


er 
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ou 
—-00 
CN xxx + D Loir — vix(22)] pre adx — 
— 00 
R=T 
co 
= f(x)ertua) gx, 
— 


Désignons les mineurs du déterminant A par le symbole Aw. On a 


+ co 
D Axg(12) J.. fax)est dx 


qg—1 
« iz(u—2%) —1— = . 
J or(x)e di= de 


—— 


Il en résulte 


+ +00 
u(x) = > Ex DE iz) NE fabertudt, 


Nous n’insisterons pas sur les conditions de validité de ces formules 
ainsi que sur la formation des solutions complètes du système (19). 


18. — Etudions maintenant le mouvement d’un point matériel 
dans les conditions héréditaires et en absence de forces. En désignant 
par xz(t) les coordonnées du point P et par gi et lx les coefficients 
d’hérédité on peut écrire les équations du mouvement 


ñn Î n { 
(20) x d jL ; ga(é — u)xi(u)du + D # 2 La(t — u)xx (u)du 
RE Ai Ê=I 


et silon a 4, <<, + h ce système prend la forme 


(204) xt) LS 1e LE Eu u)xi(u)du + D JL: PA ( (é — U | a0 (u)du. 


R=T Er" 


Imposons une condition supplémentaire : supposons que le mou- 
vement soit rectiligne ; on a par conséquent 


Li gr (to) té) — (6)] 
Xi(Ë) = Xifto) + 9(%) À 
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La substitution dans l'équation (20 bis) nous donne 


x! (to)e"(é) fe à: #7! (bo)e(£o) 1 en 
HET 2 E70 w/(o) hf Lin(é u)du = 


nn PATA) Î _ 4) | é 
+2 J) J, p(u)gin(é == u)du + > À : T6 ) oa(é u)du, 


Or les valeurs initiales sont tout à fait arbitraires. On obtient donc 


t 
(21) J eut — du = 0 (RE T, 2,00) 
Lo 
(22) Je sega — md + flat u\lafé — u)du = 0 GR) 
lo 


d 
(23) @"(£) =). o(U) gré = u)du f o'(u)lx(é — u)du. 


Il résulte des équations (21) que tous les coefficients d’hérédité 
gi sont nuls, et alors les équations (22) montrent que 


Li — 10 (2 == R). 
On voit facilement que 
L,, ER = 


et l'équation (23) prend la forme 
(24) p(Ë) — is D'(u)U(E — u)du (lo LÉ Lt + hi. 
0 
On voit de même que 
(24?) @"(E) — fs _. uit — u)du É> tb + h). 
L’équation (24) se transforme facilement : on a 
Q(E) = #)| à — pr (É— t — ut(du | + ( — bo)" (bo) + 


ri 1—u 
+ J o(u)du: Je l(z)dz. 
lo Oo 
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Cette équation de VOLTERRA se résout sans difficultés. En désignant 
par R(f — 2) le noyau résolvant de /{z), on a 


N A nt 11) 
200 = 8h) + 6) } 1 — he + Dar Rare 
© 
Ces formules montrent que l’on doit prendre la fonction /(2) néga- 
tive. Soit pour fixer les idées /(z) — — 1e, 
En ce cas R(z) — — e—(\+up, 


+ BE 4) Ar — ATH EH)] | 

P(E) = po) + P'(b0) NE Tu CME 
et si l’on admet que l’action héréditaire n’est pas bornée on voit que 
le mouvement va en ralentissant et que la vitesse se tend asympto- 


: PAL ! ia 

tiquement vers la limite o'(5) SEE PTE 

On arrive facilement aux mêmes conclusions dans le cas général. 
L'équation caractéristique de (24655) 


h 
a — à f Uu)e “du 
+, O 

a une racine &« — O0, une racine négative et une infinité de racines 
complexes. Les parties réelles de ces racines sont négatives si l’on a 


lu) < o 


k 
— [ UKu)du < 7. 
«/ 0 


19. — En terminant cette étude je veux revenir sur la question 
de l’infinité de solutions et leurs rapports avec la réalité. On est parfois 
tenté de dire que le problème réel ne peut avoir qu’une seule solution 
et que si le travestissement mathématique d’un problème physique 
en donne plusieurs c’est que ce travestissement a été mal choisi, 
On ne peut nier l'exactitude de cette réflexion dans certains cas mais 
il ne faut pourtant pas oublier les paroles de HEGEL : « tout ce qui 
est réel est raisonnable, et tout ce qui est raisonnable est réel ». Les 
conceptions de notre esprit ne sont pas des êtres imaginaires : elles 
ont parfois plus de réalité que maintes choses réelles. Notre appareil 
analytique, ainsi que toute la science, est le produit supérieur de l’évo- 


LOIR 
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lution biologique et de la vie sociale de l'humanité ; il est mieux 
adapté et plus intimement lié à la réalité des choses et aux besoins de 
notre conquête de la nature que l’on ne le croit ordinairement. Et s’il 
nous fournit quelque chose en plus de ce qu’on lui demande, il faut 
toujours attentivement chercher une réalité correspondante à ce 
surplus, une réalité qu’on trouve toujours à force de recherches et 
d’interprétations. Et nous avons toutes les raisons de croire que l’infi- 
nité de solutions n’est pas un défaut de ces équations mais plutôt 
une aubaine. Les expressions du genre 


po X+= 
gtx — tp(#)dé, [ |, gtx — he(bdt 
© "x—h VX 

se rencontrent dans plusieurs problèmes de la philosophie naturelle. 
On les trouve par exemple dans les problèmes apparemment techniques 
de la théorie des erreurs, d’interpolation, de la théorie des appareils 
enregistreurs, de la théorie de recherche de périodes, en un mot 
partout où il faut dégager les faits essentiels ou les lois empiriques des 
ensembles statistiques. Un fait observé est toujours une moyenne, 
en tout cas un résultat des opérations statistiques que l’on effectue 
parfois consciemment, mais dans la plupart de cas sans comprendre 
leur nature statistique. Or, la même moyenne observable peut corres- 
pondre à une infinité des distributions de faits. N'oublions pas d’autre 
part que les équations intégrales aux limites infinies sont très em- 
ployées dans la mécanique statistique pour la recherche des fonctions 
de distribution, et dans la mécanique ondulatoire. Or précisément 
dans le premier genre de problèmes la multiplicité de solutions peut 
correspondre à la multiplicité de lois de distribution possibles. Remar- 
quons que l'étude des fluides montra toute l'importance de la turbu- 
lence, mais l'introduction de la turbulence dans les équations de 
l’hydrodynamique est encore une opération statistique nécessitant 
À th 

n “dt, etc. Il n’y a donc rien d’éton- 

2 
nant à ce qu'une équation intégrale ait une infinité de solutions. 
La fonction arbitraire qui entre dans la composition de la solution 
complète est selon moi une circonstance favorable correspondant 
à la réalité des choses et nous permettant une certaine liberté de 


l'usage des intégrales # — . J 
t 
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mouvements. C’est peut être un défaut au point de vue des axioma- 
tistes extrêmes ; la conception statistique de la nature comporte 
heureusement toujours de l’imprévu et de l’imprévoyable pour une 
machine logique si perfectionnée qu’elle soit. On peut bien axiomatiser 
une théorie physique mais les axiomes utilisés ne seront que des 
postulats ou des conventions et encore faut-il toujours introduire, 
en les établissant, des postulats cachés comme on introduit sans 
nombre en établissant par exemple l'équation de BoLTZMANN. On peut 
construire un appareil logique mais ce sera un appareil d'usage très 
limité. Je ne veux pas dire par cela que les machines logiques ne soient 
pas utiles ; je veux dire seulement que dans la conquête et l'étude 
de la nature il y aura toujours place à l’activité créatrice humaine. 
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